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Abstract 



The physics of interacting Fermi systems is extremely sensitive to the energy scale. 
Of particular interest is the low energy regime where correlation induced collective 
behavior emerges. The theory of interacting Fermi systems is confronted with the oc- 
curence of very different phenomena along a continuum of scales calling for methods 
capable of computing physical observables as a function of energy scale. 

In this thesis, we perform a comprehensive renormalization group analysis of two- 
and three-dimensional Fermi systems at low and zero temperature. We examine sys- 
tems with spontaneous symmetry-breaking and quantum critical behavior by deriving 
and solving flow equations within the functional renormalization group framework. 

We extend the Hertz-Millis theory of quantum phase transitions in itinerant fermion 
systems to phases with discrete and continuous symmetry-breaking, and to quantum 
critical points where the zero temperature theory is associated with a non-Gaussian 
fixed point. The order parameter is implemented by a bosonic Hubbard-Stratonovich 
field, which -for continuous symmetry-breaking- splits into two components cor- 
responding to longitudinal and transversal Goldstone fluctuations. We compute the 
finite temperature phase boundary near the quantum critical point explicitly including 
non-Gaussian fluctuations. 

We then set up a coupled fermion-boson renormalization group theory that cap- 
tures the mutual interplay of gapless fermions with massless order parameter fluc- 
tuations when approaching a quantum critical point. As a first application, we com- 
pute the complete set of quantum critical exponents at the semimetal-to-superfluid 
quantum phase transition of attractively interacting Dirac fermions in two dimen- 
sions. Both, the order parameter propagator and the fermion propagator become non- 
analytic functions of momenta destroying the Fermi liquid behavior. 

We finally compute the effects of quantum fluctuations in the superfluid ground 
state of an attractively interacting Fermi system, employing the attractive Hubbard 
model as a prototype. The flow equations capture the influence of longitudinal 
and Goldstone order parameter fluctuations on non-universal quantities such as the 
fermionic gap and the fermion-boson vertex, as well as the exact universal infrared 
asymptotics present in every fermionic superfluid. 
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1 



Introduction 



This thesis is concerned with low- and zero-temperature properties of correlated 
fermion systems in two and three dimensions. Interaction effects accompanied by 
quantum and thermal fluctuations play an important role and produce fascinating phe- 
nomena such as, for example, new scaling laws in the vicinity of quantum critical 
points and exotic forms of superfluidity in cold atomic gases. 



1.1 Experiments 

Before presenting the outline of this thesis, we refer to two recent hallmark experi- 
ments which are related to computations performed in the present work. 

1.1.1 Phase boundary close to a quantum critical point 

First, we discuss the recent measurement of the finite temperature phase boundary 
of the antiferromagnetically ordered phase in BaCuSiiOg performed in Stanford (Se- 
bastian 2006). In this compound, the line of finite temperature phase transitions re- 
mains continuous down to the lowest temperatures and ends at a quantum critical 
point (QCP). In the vicinity of the QCP, the shape of the phase boundary becomes 
universal and follows the power-law: ~ (6 - dc)*^, where 6 is the non-thermal con- 
trol parameter -in this case the external magnetic field. The so-called shift exponent 
i/r is independent of microscopic details and depends only on dimensionality (d) and 
the additional dimensionality incurred from quantum fluctuations (z): 

a + z-2 

as first derived by Millis (1993). The relevant excitations of the material measured by 
the Stanford group (Sebastian 2006) are spin-dimers which can be described by an 
interacting Rose gas undergoing Bose-Einstein condensation at the critical tempera- 
ture. The original spin degrees of freedom exhibit antiferromagnetic order when the 
(collective) spin-dimer gas condenses. In such cases, the dynamical exponent is z = 2. 



2 1 Introduction 




Fig. 1.1. Phase boundary of the antiferromagnetic phase in BaCuSi20e close to the 
QCP (Sebastian 2006). (a) Values of the shift exponent (here v) obtained from fitting 
experimental points on the phase boundary. The data approach v = 2/3 in the in- 
termediate regime, and there is a distinct crossover toward v = 1 before the QCP is 
reached, (b): Estimates of the location of the QCP, Hd, obtained along with v during 
the fit. (c): Best fits to the phase boundary in the intermediate and low temperature 
regime on a logarithmic scale. 



In Fig. ll.ll (c). data from torque magnetometry for the phase boundary is exhibited, 
clearly in agreement with Eq. (|l.ll) . The area coloured in blue is the paramagnetic 
region and the yellow area corresponds to the antiferromagnetically ordered phase. 
Upon lowering the temperature when approaching the QCP, there is a crossover from 
three-dimensional scaling to two-dimensional scaling. This has been argued to be a 
consequence of effective decoupling of the 2J-layers of the material at very low tem- 
peratures and is specific to the geometrically frustrated lattice structure of BaCuSiiOg. 

In chapters |5] and [6l we compute zero- and finite-temperature properties of cor- 
related systems close to a QCP. The interplay of thermal and quantum fluctuations 
makes this an interesting subject for theoretical studies. Among other things, we in- 
troduce a new way to compute the shift exponent coming from the symmetry-broken 
phase with and without Goldstone modes by the use of modem renormalization group 
equations which continuously connect the quantum fluctuation dominated regime di- 
rectly at the QCP to the regime dominated by non-Gaussian classical fluctuations 
further up the Tc-line. 
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1.1.2 Phase diagram of an attractive two-component Fermi gas 

As the second experiment, we make reference to the seminal measurement of the 
phase diagram for a spin-polarized Fermi gas performed by the Ketterle group at MIT 
(Shin 2008). In this experiment, a different amount of spin-up and spin-down com- 
ponents of ^Li are loaded in an optical trap. By applying an external magnetic field 
close to a Feshbach resonance the population imbalanced gas is tuned to the reso- 
nantly (attractively) interacting regime, forms bosonic spin- singlets, and condenses to 
the superfluid state upon lowering the temperature. The non-thermal control parame- 
ter here is the population imbalance or polarization a = («-[•- n^) / (n-f -I- nj^) , where 
T and i refer to the two spin components with densities |. Increasing cr deprives the 
attractively interacting fermions of their respective opposite-spin pairing partner thus 
inhibiting the superfluid pairing beyond a critical polarization crdT). 

As shown in the left plot of Fig. II. 2[ an increasing polarization indeed suppresses 
the critical temperature and destabilizes the superfluid phase at low temperatures. In 
the inner region of the trap, a superfluid core forms and the spin polarization shows a 




Spin polarization Radial position. 



Fig. 1.2. Left: The cr - T -phase diagram for a homogeneous spin-polarized Fermi 
gas with resonant interactions (Shin 2008). The critical polarizations cr^ (black solid 
circles and square) and cr^ (grey solid circles) are displayed along the local T/Tf^, 
with Tff the Fermi temperature of the majority component in the spin-up state, at 
the phase boundary. The yellow area ( cr^ < cr < cTc) represents a region with phase 
separation until cr^ = cTc at the tricritial point. The triangle with ref. 8 marks Tc for 
the unpolarized Fermi gas at unitarity from Burovski et. al (2006) and the square with 
ref. 10 marks the position of the first order quantum phase transition from Lobo et. al. 
(2006). Right: At low T, the sample has a three-layer structure (Shin 2008): the core 
region (0 < r < Rc) of a fully paired superfluid with n^ = n^; the intermediate region 
(Rc <r< Ri) of a partially polarized normal gas; and the outer region (Ri < r < R^) 
of a fully polarized normal gas. 
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discontinuity at the boundary of the superfluid core r = R^, a signature of the phase 
separation of a superfluid and a normal gas (Bedaque 2003). The critical polarization 
(Tc = limr->«+cr(r) represents the minimum spin polarization for a stable normal gas; 
o-s = lim^_>K-cr(r) represents the maximum spin polarization for a stable superfluid 
gas as illustrated in on the right plot in Fig. 11.21 At o-cq ~ 0.36, in the bottom right 
corner of the left plot in Fig. 11.21 the Clogston limit of superfluidity which is a first 
order quantum phase transition was experimentally verified. 

One way to achieve a continuous quantum phase transition in attractive fermion 
systems involves cold atomic gases loaded in a honeycomb optical lattice (see Zhao 
2006 and references therein) and possibly graphene (Uchoa 2007, Castro Neto 2008). 
The attractive Hubbard model on the honeycomb lattice exhibits, at half-filling, a 
quantum critical point between a semimetal with massless Dirac fermions and an s- 
wave superconductor. 

Fermionic systems displaying superfluidity form a cornerstone of this thesis. In 
chapter |7l we investigate attractively interacting Dirac fermions in the vicinity of a 
continuous semimetal-to-superfluid quantum phase transition and in chapter [H we 
consider the superfluid ground state of the attractive Hubbard model at quarter-filling. 

1.2 Thesis outline 

This thesis is structured in two parts. In the first part, the methodological framework 
is presented. The second part contains various applications of these methods and con- 
cepts. In detail, the chapters in Part U have the following contents: 

• In chapter |2l we introduce the fundamental concepts necessary for an understand- 
ing of the applications in the second part. First, the standard model of interacting 
Fermi systems. Landau's Fermi liquid theory, and instabilities thereof are phrased 
in the language of the renormalization group (RG). The notion of spontaneous 
symmetry-breaking with a focus on Goldstone bosons is introduced. These mass- 
less bosonic excitations arise when a continuous symmetry is broken, as is the 
case for example in magnetic or superfluid systems. Then, we -after having pre- 
sented basic features of phase transitions in general- expose the conventional RG 
approach to quantum critical systems referred to as the Hertz-Millis theory. In 
transiting to chapter [2] we recapitulate the underlying ideas of the Wilsonian RG 
approach to quantum field theory. 

• In chapter [31 the functional integral approach to many -particle systems containing 
fermionic and bosonic fields is laid out. At the heart of the chapter is the deriva- 
tion of the functional RG framework in its one-particle irreducible representation 
geared toward systems with spontaneously broken symmetry. 
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The salient points of Part IJ are summarized in chapter |4l The applications presented 
in Part HI] are organized as follows. 

• In chapter [51 we extend the Hertz-Millis theory for quantum criticality to phases 
with broken discrete symmetry. Using the Hertz action, we compute the RG flow of 
the bosonic propagator and the effective potential with flow equations derived from 
the functional RG framework. Difl^erent dynamical exponents are distinguished 
and the shape of the phase boundary at finite temperature is computed for each. 
Large parts of this chapter have been published previously in Phys. Rev. B 77, 
195120 (2008) and the results for the Quantum Ising model are published in Phys. 
Rev. B 80, 085108 (2009). 

• In chapter[6l we analyze quantum critical systems with Goldstone modes. We em- 
ploy an exclusively bosonic effective action with two propagators: one associated 
with the transverse (Goldstone) dispersion and the other stemming from longitudi- 
nal fluctuations. The interactions between longitudinal and Goldstone fluctuations 
lead to strong renormalizations of the longitudinal propagator away from and at 
the phase boundary where both modes become degenerate. We compute RG flows 
to determine the phase boundary at finite temperatures in three dimensions. 

• In chapter |71 the main inconsistency of the conventional Hertz-Millis approach 
-that massless fermions are integrated out in one sweep- is cured by explicitly in- 
cluding fermions into a coupled RG flow for (Dirac) fermions and order parameter 
fluctuations. This way, we assess the mutual interplay of both types of fluctuations 
for the semimetal-to-superfluid quantum phase transition in two dimensions. At 
the QCP, we find that both the order parameter propagator and the single-particle 
propagator are non-analytic functions of frequency and momenta signalling the 
breakdown of the Fermi liquid. An improved version of this chapter is published 
in Phys. Rev. B 81, 125103 (2010). 

• In chapter [8l the low energy behavior of the attractive Hubbard model at quarter- 
filling as a prototype for systems with a superfluid ground state is analyzed. Various 
non-universal quantities such as the fermionic gap are computed. At the same time, 
the universal infrared behavior correctly emerges from the flow of the coupled 
fermion-boson action in agreement with the exact behavior of an interacting Bose 
gas. This chapter has been published in Phys. Rev. B 78, 014522 (2008). 

Part Hn] contains the conclusions of this thesis, criticism, and a description of future 
research projects. The numerical procedure developed to solve the flow equations is 
explained in Appendix lAl 
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Theoretical Framework 



2 



Underlying concepts 



2.1 Fermi liquid instabilities 

This section investigates conditions under which an interacting Fermi system might 
not be described well anymore by its original degrees of freedom, namely fermionic 
(quasi-) particles. To lay a basis, we first describe the standard model of interacting 
fermions: the Landau Fermi liquid theory (for references and review see Nozieres 
1964). This theory is formulated entirely in terms of fermionic quasi-particles -bare 
fermions dressed by interaction effects, hence the prefix quasi. The validity of the 
Landau Fermi liquid theory rests upon fulfillment of the following conditions: 

• the existence of a well-defined Fermi surface, that is, a (d- l)-dimensional hy- 
persurface in momentum space, to the vicinity of which the low-energy fermionic 

excitations are restricted to, 

• the existence of quasi-particles with a one-to-one correspondence to the non- 
interacting particles by adiabatically turning on the interaction. 

In the functional integral terminology of the renormalization group, the Fermi liquid 
is said to be described by the fixed point or scale-invariant action: 

i'<A ^ r-<A 

S [iff, lf\ = ijfk {iko - Vkp^r) 'l^k-T; fkp,k{:h-q/2lJfk'+q/2ll^k'-q/2(//k+q/2 , 

Jk ^ Jk,k',q 

(2.1) 

where jy. = T J^^^ / comprises (Matsubara) frequency summation and momen- 
tum integration restricted to a shell of width 2A around the Fermi surface, if/, ij/ are 
renormalized Grassmann fields rescaled by the fermionic Z^p-factor, k,. is the momen- 
tum deviation in radial direction from the Fermi surface, is the Fermi velocity 
projected on the Fermi surface, and f^^^^, the Landau function obtained as the (non- 
commutative) limit lim^j_>o limq^o with q/^o ^ of the two-particle vertex (for a 
review see Metzner 1998). Here and in the following, we employ the four-vector no- 
tation k = (ko, k). Considering only processes with zero-momentum transfer by taking 
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the limit ^ — > 0, the Hartree mean-field theory of this action leads to the energy func- 
tional phenomenologically postulated by Landau 

6E [6n\ = [ Vkp^r^^k + ^ / /kF,k;,5nkF<^«^ , (2.2) 

Jk ^ Jk,k' 

with (5nk the quasi-particle distribution function. In addition to constituting the stan- 
dard model of interacting fermions, the merit of Landau's theory lies in its predictive 
power for systems with spherical Fermi surfaces, the epitome being liquid ^He. Upon 
expanding the Landau function /kj.,k;.. in Legendre polynomials, one can -with knowl- 
edge of only the first few Legendre coeflicients- compute a larger number of physical 
observables such as, e.g., the specific heat or compressibility (Negele 1987). 

An important trademark of the Fermi liquid is the renormalized single-particle 
propagator determined by the quadratic part of the fixed point action as 

Gf{k) = ^ , (2.3) 

iko - Vk^kr 

which describes stable quasi-particle excitations with a velocity Vk^. Note that the 
interaction-induced reduction in quasi-particle weight by Zkp has been absorbed into 
renormalized field variables. Standard higher-order corrections to the fermion self- 
energy of Landau's theory are expected to be quadratic in frequency and temperature 
(Chubukov 2003). Fresh research by Millis, Chubukov and others, however, revealed 
-already for the generic weakly and locally interacting Fermi gas- non-analytic cor- 
rections to the Fermi liquid behavior of the susceptibilities and the specific heat (Chi- 
tov 2001, Chubukov 2003, Chubukov 2005). 

Indeed, much of the more recent research in correlated fermion systems has been 
devoted to systems where the Fermi liquid paradigm breaks down and other theoreti- 
cal descriptions must be invoked (Stewart 2001). Often, tendencies to destabilize the 
Fermi liquid are already visible in the bare bosonic response functions of the particle- 
particle (pp) channel: 

npp(qo,q) = - J^Gfo(k + q)Gfo(-k)oc .,;^^^(>«^ (2.4) 

and the particle-hole (ph) channel, 

Ilphiqo, <D = - J Gfo(k + q) G/o (k) oc -^si^c^^^p^ , (2.5) 

where Gjl = iko - ^k is the bare propagator with ^k the single-particle energy relative 
to the Fermi surface. The magnitudes of i7ph(<?o, q) and n^{qQ, q) in conjunction with 
the sign and momentum structure of the interaction may point at the low-energy fate 
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of the system. If, in a lattice structure for example, the interaction is local and repul- 
sive and max{77pp(^05 q)^ ^ph(^o> Q)} occurs at /Zph(0, Q = 2kf ), the system may form 
a charge- or spin-density wave eventually leading to charge or magnetic order. On 
the other hand, strong particle-hole fluctuations in some cases enhance pairing phe- 
nomena in the particle-particle channel via the Kohn-Luttinger eff"ect. Therefore, the 
interplay of both channels has to be assessed to make a statement about the true low- 
energy fate of the system. Moreover, possible infrared singularities in one or some- 
times even multiple scattering channels invalidate perturbative approaches. For certain 
geometries, such as a nested Fermi surface or when the Fermi surface intersects the 
van-Hove points, 77ph exhibits a logarithmic divergence for small energies. ZTpp also 
diverges logarithmically as a function of energy for an inflection-symmetric Fermi 
surface in any dimension (for a review see Abrikosov 1975). 

A situation of competing Fermi liquid instabilities where both, ZTph and ZTpp are 
important, occurs for example in the repulsive Hubbard model and can be ana- 
lyzed within a fermionic functional RG treatment which takes into account both, 
the particle-particle and particle-hole channel (and their mutual interplay) on equal 
footing (Zanchi 2000, Halboth 2000, Honerkamp 2001). In these RG flows, it is the 
summed up non-universal contributions to the flow that determine which of the sev- 
eral universal fixed points (Fermi liquid, antiferromagnetic, superconducting, etc.) the 
system scales towards. 

A drawback of the hitherto employed fermionic functional RG is its failure to pen- 
etrate phases with broken-symmetry as the vertex functions grow large at a critical 
scale thus invalidating the truncation devised for weakly to intermediately coupled 
systems. If one attempts to go beyond this scale, the vertices eventually diverge and 
the flow has to be stopped. To continue the flow into phases with broken symmetry 
is clearly desirable because many systems have symmetry-broken ground states and 
many experiments are able to verify or falsify quantitative predictions from calcula- 
tions in the symmetry -broken phase such as, for instance, the magnitude and symme- 
try properties of the superconducting gap in the cuprates. 

At present, two remedial approaches have been developed. In the first approach, a 
small initial anomalous self-energy for one or multiple channels is added to the bare 
action and as such offered to the flow (Salmhofer 2004, Gersch 2005, 2006, 2008). If 
now the symmetry-breaking occurs in one of the offered channels, this initial anoma- 
lous self-energy regularizes the vertices around the critical scale and the flow can be 
continued into the symmetry-broken phase. This method requires no decoupling of 
the interaction. 

In the second approach, the fermionic interaction is decoupled with a suitable 
Hubbard-Stratonovich field and one subsequently deals with a coupled fermion-boson 
theory (Baier 2004, Strack 2008). This approach yields easy access to the phenom- 
ena associated with spontaneous symmetry -breaking alluded to in section [2^ and we 
employ it in chapter [S] of this thesis in the context of fermionic superfluids. 
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2.2 Spontaneous symmetry-breaking 

Spontaneous symmetry-breaking occurs if the energetically lowest state of a system 
described by a quantum or statistical field theory exhibits different symmetry proper- 
ties than the original Lagrangian. The importance of spontaneous symmetry-breaking 
for physics can hardly be overestimated as has been recognized by the Nobel commit- 
tee in the last 50 years with a series of Nobel prizes. At the time of writing, Nambu has 
been awarded the Nobel prize "for the discovery of the mechanism of spontaneously 
broken symmetry in subatomic physics" (Nambu 2008). Quite appropriately, we will 
in parts employ Nambu's formalism in chapter[8]of this thesis. 

In elementary particle physics as well as condensed matter physics, a stunning 
range of observed phenomena finds an explanation with the help of symmetry con- 
siderations (Negele 1987, Peskin 1996). Perhaps most fundamentally, the Anderson- 
Higgs mechanism (Negele 1987) bridges the gap between quantum excitations of the 
vacuum in high-energy physics and statistical many-particle systems at low tempera- 
tures. In the former context, this 
scenario for dynamical mass generation 
of fundamental particles such as gauge 
bosons: by (spontaneous) condensation 
of the Higgs field the order parameter 
thereof generates a mass term for par- 
ticles covariantly coupled to the Higgs 
field. In the latter context of statistical 
physics, it was Anderson who first pro- 
posed this mechanism in 1958 to explain 
the Meissner effect in superconductors 
(Negele 1987): the photon (gauge bo- 
son) mediating the electromagnetic in- 
teraction becomes massive from the su- 
perconducting condensate (Higgs field). 

Goldstone bosons emerge when spontaneously breaking a continuous symmetry. 
In the low-energy sector of massless Quantum Chromodynamics, for example, the 
formation of mesonic two-quark bound states spontaneously breaks the chiral sym- 
metry of the original Lagrangian. As a result, the emergent pions can be understood 
as the Goldstone bosons of massless QCD (Peskin 1996). Many striking properties of 
low-temperature condensed matter systems such as magnets, superfluid Helium, and 
superconducting materials are attributed to Goldstone bosons. In superfluid systems, 
Goldstone bosons are an indicator of long-range phase coherence. In this case, the 
order parameter contains two degrees of freedom, one associated with the amplitude 
and the other with the phase of the order parameter. Above the critical temperature, 
the Hamiltonian is invariant under the rotational symmetry group 0(2). Below the 



mechanism is the generally accepted 




Fig. 2.1. Potential for spontaneous break- 
ing of continuous 0(N)-symmetry, drawn 
for the case N = 2. Fluctuations along the 
trough in the potential correspond to the 
Goldstone excitations (Peskin 1996). 
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critical temperature, the potential attains the form shown in Fig. 12. H and the energeti- 
cally lowest state is degenerate with the respect to phase transformations (i.e. location 
in the trough). Generalized to an A'^-component order parameter, the original 0(N)- 
symmetry is hidden leaving only the subgroup 0(N - 1), which rotates the Goldstone 
modes among themselves (Peskin 1996). Note that the form of the potential must re- 
spect the full 0(A^)- symmetry even below the critical temperature, but the field evalu- 
ated at its minimum breaks and reduces the symmetry. Theoretically, Goldstone modes 
are difficult to cope with in perturbation theory as the massless propagators cause se- 
vere infrared singularities in low enough dimensions necessitating a renormalization 
group (RG) treatment. 

In this thesis, the physics of Goldstone modes plays a central role. In chapter[6l we 
consider a quantum critical point adjacent to an ordered phase with Goldstone bosons. 
In magnetic systems, for instance, the original 0(3) symmetry is broken down to 0(2) 
and the Goldstone bosons are collective, gapless fluctuations in the direction of the 
electron spins with respect to which the ground state is degenerate. For antiferromag- 
netically ordered spins, these excitations form spin waves. 

In chapter m we compute the RG flow of various observables in fermionic su- 
perfluids. Here, the Goldstone bosons correspond to fluctuations in the phase of the 
superfluid order parameter and the original 0(2) symmetry is broken spontaneously. 

2.3 Quantum criticality 

In certain cases, the departure from the Fermi liquid occurs in the vicinity of a quan- 
tum critical point (QCP) where spontaneous symmetry-breaking occurs at zero tem- 
perature. This quantum criticality induced breakdown of the Fermi liquid has attracted 
enormous attention in recent years (Stewart 2001, Loehneysen 2007). This section 
recollects the most popular approach to quantum critical fermion systems, the Hertz- 
Millis theory (Hertz 1976, Millis 1993) and its shortcomings, after having recapitu- 
lated similarities and differences between classical and quantum phase transitions. 

The amount of thermal agitation in a multi-particle system can severly alter its 
macroscopic properties although it is described by the same particles, the same inter- 
actions, and the same Hamiltonian (Goldenfeld 1992). The point separating the two 
distinct -for example para- versus ferromagnetic- phases by tuning the temperature 
to Tc, is referred to as the critical point if the correlation length among the constituent 
particles diverges. At T^, thermal fluctuations are so strong that the system becomes 
self-similar on all scales, i.e., it becomes scale-invariant. The free energy develops 
non-analyticities which can be subsumed under a set of critical exponents two of 
which are independent related to the existence of two relevant scaling fields, namely 
temperature and the external symmetry -breaking field (Goldenfeld 1992). It now turns 
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out that -rather remarkably- lavish reservoirs of physical systems can be described 
by the same set of critical exponents distinguished only by their most basic symmetry 
properties and spatial dimensionality into universality classes. 

In recent times, one has come to en- 
visage systems where continuous phase 
transitions occur at zero temperature. 
Here, a non-thermal control parameter 
drives the phase transition and zero- 
point quantum fluctuations trigger quan- 
tum critical scaling of various quanti- 
ties often described by diff"erent expo- 
nents than their classical counterparts 
(Sachdev 1999, Loehneysen 2007). Such 
a quantum critical point lies between 
neighboring ground states -the generic 
situation is exhibited in Fig. 12. 2[ A pe- 
culiarity of quantum phase transitions 
is their increased eff'ective dimensional- 
ity incurred by quantum fluctuations in 
the direction of imaginary time (Hertz 
1976). In the functional integral formal- 
ism, the frequency term in the quadratic part of the action, + q^) 4> for 
z= 1,2, and (|i2|/|q|^"^ -i- q^) for z > 3, with (p being the bosonic order parameter 
field, indicates that characteristic frequencies scale as i3 ~ |q|^. When approaching the 
transition at T = 0, both the order parameter correlation length ^ and correlation time 

diverge as function of control parameter: 




Fig. 2.2. Generic phase diagram in the 
vicinity of a continuous quantum phase 
transition. 



(2.6) 



At finite temperature, however, the {d + z)-dimensional quantum system has finite 
length in the time direction, Lt ~ I /T. In Fig. 12.21 we observe that the existence of a 
QCP shapes large portions of the phase diagram. The black line of finite-temperature 
classical phase transitions separates the thermally disordered phase from the ordered 
phase and terminates at the QCP. Its shape is determined by the shift-exponent as al- 
luded to in chapter [T] A completely different regime is the high-temperature regime 
above the QCP where ^ Lt, the critical singularity is cut off by the finite tem- 
perature, and the boundaries are the crossover lines T \6 - 6^^^. Strikingly, in this 
quantum critical regime one can measure power-law dependencies of physical observ- 
ables up to rather high temperature with exponents not assuming the expected Fermi 
liquid values (Loehneysen 2007). 

For a computation of these non-Fermi liquid power-laws in itinerant fermion sys- 
tems, Millis (1993) has set up RG equations capitalizing on the fact that the zero- 
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temperature fixed point for the order parameter is Gaussian as long as J + z > 4. The 
procedure is the following. The microscopic four-fermion interaction is decoupled 
with a Hubbard-Stratonovich field corresponding to the incipient order. The fermions 
are integrated out from the functional integral in one sweep, but leaving their mark 
on the subsequent, usually polynomial, expansion in powers of the order parameter 
field. One then integrates out order parameter fluctuations on a Gaussian level (Mil- 
lis 1993). Already this results in a rich phase diagram with some materials fulfilling 
the Hertz-Millis exponents but others not (Stewart 2001, Sebastian 2006, Loehneysen 
2007). For certain cases, such as the SU(2)-symmetric antiferromagnetic channel, this 
procedure is ill-defined as the coefficients are singular functions of frequency and mo- 
mentum (Chubukov 2003, Belitz 2005, Rech 2006, Loehneysen 2007). Experimental 
discrepancies and the questionable procedure to integrate out gapless fermions before 
considering order parameter fluctuations suggest that there is room for improving the 
Hertz-Millis theory. 

In this thesis, we go beyond the Hertz-Millis approach in three ways. First, we 
show how to account for discrete symmetry -breaking and non-Gaussian classical fluc- 
tuations by use of truncated functional RG equations in chapter [5l Second, possible 
eff'ects of Goldstone modes and their fluctuations on the quantum phase transition 
when breaking a continuous symmetry are analyzed in chapter |6l Third, we present 
a coupled fermion-boson RG that is able to assess the interplay of gapless fermions 
with massless order parameter fluctuations in chapter |7] 

2.4 The renormalization group 

The perturbative expansion of most quantum field theories is plagued by unbounded 
expressions. Either the integrands are singular for small momenta and the phase space 
volume in low dimensions does not shrink fast enough to tame the infrared diver- 
gence or the phase space volume blows up for large momenta leading to ultraviolet 
divergences in higher-dimensional continuum theories. These divergences are either 
of physical origin signalling for instance the proximity to a phase transition with truly 
infinite observables, or of technical origin as result of an inefficient organization of 
the perturbation theory with all physical observables being in fact finite. 

In the 1950's and 60's, the obnoxious presence of these singularities in calculations 
proved a pertinacious problem to solve and severely hindered progress in statistical 
and quantum field theory. Then, in the beginning of the 1970's, Wilson -at the top of 
a scientist iceberg- reformulated the functional integral in terms of difi'erential equa- 
tions along the continuous flow parameter A. Typically, but not always, A is associated 
with the characteristic kinetic energy scale implemented by a cutoff" for the propaga- 
tor. This cutoff carves out the singularity from the spectrum and excludes modes with 
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lenergyl < A from the functional integral. This way, the singularity is approached peu- 
a-peu only at the end of the functional integration when A ^ and the renormalized 
theory is delivered. The solution of these so-called RG or flow equations enabled, 
among countless other groundbreaking discoveries, a controlled computation of criti- 
cal exponents for bosonic field theories (for a review see Wilson 1974). The extension 
of these ideas to non-Abelian gauge theories involving fermions and bosons led to the 
discovery of asymptotic freedom (Politzer 1973, Gross 1973). 

Purely fermionic systems in one dimension were considered with RG methods 
even before Wilson's Nobel prize in 1982 by Solyom (1979). Originally taken up 
by mathematical physicists (see references in Salmhofer 1999), the RG was subse- 
quently applied to interacting Fermi systems also in higher dimensions (for a review 
see Shankar 1994, Metzner 1998, and Salmhofer 1999). The existence of a Fermi sur- 
face greatly enriches the analysis as the fermion propagator becomes unbounded on 
a (d - l)-dimensional hypersurface to be compared with the much simpler situation 
for bosons where the propagator becomes large only at one point (the origin) of mo- 
mentum space. Rigorous proofs involving phase space arguments and the curvature of 
the Fermi-surface have established the RG as a mathematically well defined starting 
point for practical computations. Therein, the low energy projection of the vertices 
onto the Fermi surface involves flowing functions to parametrize transversal momen- 
tum dependences parallel to the Fermi surface which remain relevant until the very 
end of the flow (Zanchi 2000, Halboth 2000, Honerkamp 2001). For this reason, the 
Fermi surface is partitioned into patches and mapped onto a grid in momentum space. 
The vertex function is then parametrized by several thousand variables and the flow 
equations are solved numerically. 

Flowing functions can be dealt with conveniently within modem functional renor- 
malization group formulations (Polchinski 1984, Wetterich 1993, Salmhofer 2001, 
for a review see Berges 2002, Metzner 2005). Three formally exact flow equations 
have been devised: (i) the Polchinski (1984) scheme for the generating functional of 
connected amputated vertex functions ("V) utilized by Zanchi (2000), (ii) the Wick- 
ordered scheme (Salmhofer 1999) where the vertex functions are the expansion coeffi- 
cients when expanding "V in Wick-ordered polynomials of the source-field utilized by 
Halboth (2000), and (iii) the one-particle irreducible (IPI) scheme (Wetterich 1993, 
Salmhofer 2001) for the generating functional of the IPI-vertex functions utilized by 
Honerkamp (2001) and in a variety of contexts in Berges (2002). 

In the next chapter, we thoroughly present the most popular, the IPI-scheme of the 
functional RG. A particular strength of this scheme is the natural inclusion of self- 
energy corrections as all internal lines are fully dressed. Further, coupled fermion- 
boson theories and spontaneous symmetry-breaking can be treated easily within this 
scheme. We will apply this scheme for the computations in Part|lI]of this thesis. 
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In this chapter, we introduce the main computational tools employed in this thesis. 
First, we set up a unifying framework for fermionic and bosonic models by defining 
superfields whose fermionic and bosonic components are distinguished by a statis- 
tics index. Second, we define the one-particle irreducible generating functional, the 
effective action. Third, the functional flow equation for the effective action is derived. 
Finally, we explain how to devise truncations with non-zero expectation values of 
field variables to account for spontaneous symmetry -breaking. 

3.1 Functional integral for quantum many-particle systems 

We consider interacting Fermi systems described by the bare action, 

ro[l^,(A] = - / iJfKG~fQiK)l^K + / V{^, Qlf/K^I^K'+Qil^K'il^K+Q , (3.1) 

JK Jk,K',Q 

where ijf and i/r are Grassmann- valued fields, the index (K = {ko, k, cr) and its corre- 
sponding integration f^^ comprises frequency, momentum and internal states such as, 
for example, the spin projection, and GJq{K) = iko-^k is the bare inverse Green func- 
tion. The Fermi surface is defined as the (d - l)-dimensional manifold = 0. V^^, q 
denotes an arbitrary and in general momentum- and spin-dependent many-body inter- 
action of the lowest non-trivial (quartic) orderQ. In many circumstances, correlations 

^ Models of the form Eq. ( 13. iP have attracted attention for decades: in one- 
dimensional systems, various phases such as the Luttinger liquid (Solyom 1979, 
Voit 1994) occur; in the context of gauge theories, in vacuum and with more elab- 
orate internal degrees of freedom, Eq. A3.1\) may describe the Thirring or Gross- 
Neveu Model (Hands 1993); on the lattice in two and higher dimensions, Eq. 
(1X7]) includes the attractive or repulsive (depending on the sign of the in that case 
momentum-independent interaction) Hubbard model often employed in contempo- 
rary condensed matter physics (Micnas 1990, Baeriswyl 1995). 
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induce fermion-pairing phenomena and one finds the description in terms of the order 
parameter useful. The simplest bosonic action is the (^"^ -model 

^ol^] = - (f>*QGl^(Q)(f)Q + VK,K',QfK-QfK'+Q^K'^K, (3.2) 

jQ Jk,K',Q 

where (p is here a complex-valued field but may in general contain components 
referred to as 0{N)-mo&Q\s (Amit 1995). G^q(^) = - (^q -I- 6t>q -I- m^) is the bosonic 
Green's function with being the dispersion, usually but not always quadratic in 
momenta, and m\ = ^"^ the mass or inverse correlation length. The frequency term is 
also not always quadratic; for the non-relativistic Bose gas one has a complex linear 
term ~ iqo and for certain quantum-critical systems Landau-damping entails |^ol or 
even |^ol/q (Hertz 1976, Millis 1993). 

Finally, models with both, fermionic and bosonic fields, may be the subject of 
analysis. In addition to the kinetic terms with G^q(2) and G~j:q{K), to lowest order we 
then have to consider the fermion-boson vertex 

ro[<^,^,(A]= / gK,Q {fQil/-K+Q^K + 4>Q^K+Q^-K) ■ (3.3) 

Jk,q 

3.1.1 Superfield formulation for fermionic and bosonic models 

To develop a unifying framework, we combine fermionic and bosonic fields in a su- 
perfield S, where fermions and bosons are distinguished by a statistics index s = b,f, 
that is, 

Sb = ^, Sf='P. (3.4) 

To account for matrix propagators necessary to describe phases with broken symmetry 
such as, for example, superfluidity, we here use the fermionic Nambu fields 

n= ij^'l^ , %={h^,ilf-ki) (3.5) 
and bosonic Nambu fields 

The fermionic and bosonic matrix propagators are then given by G f{k) = -{W^ and 
Gbiq) = -{'I*q^q) , respectively. The superpropagator G(^) = -(SqSq) is diagonal in 
the statistics index. 

Note that the specific choice of the fermionic and bosonic Nambu fields depends 
on the physical situation under investigation and can be adjusted straightforwardly. 
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3.1.2 One-particle irreducible generating functional 

Functional integration of the bare action containing fermionic and bosonic fields, 
ro[<S,<S], and subsequent functional differentiation with respect to source fields S' 
and S' linearly coupled to S and S yields the connected m-particle Green functions 
(Negele 1987) 



Gm {K[, K'^; Ki, Km) - -(SK[-SK'^SK„,---'SKt)c 

gm Qm 



, (3.7) 



where (...)( denotes the connected average and the generating functional is obtained 
from the logarithm of the partition function 



(3.8) 



where the bracket (., .) is a shorthand notation for the inner product for superfields. 
The Legendre transform of G[S', S'] is the effective action 

nS, S] = £G[S', S'] = G[S', S'] - (S',S) - (S, S') , (3.9) 

where the conjugated field variables are 

r[tS,<S] generates the one-particle-irreducible (IPI) Green functions. Topologically, 
all diagrams are generated where cutting one line, irrespective if fermionic or bosonic, 
does not separate the diagram into two disconnected parts (Negele 1987, Amit 2005). 



3.2 Flow equations 

In this section, we derive the exact flow equation (Wetterich 1993, Salmhofer 2001, 
Berges 2002, Metzner 2005, Enss 2005) for the effective action, that is, the IPI gen- 
erating functional, adapted to superfields. 

3.2.1 Exact flow equation 

The exact flow equation is a reformulation of the functional integral as a functional 
differential equation. It describes the evolution of the effective action as a function of 
a flow parameter A, usually a cutoff. The cutoff can be implemented by adding the 
regulator term quadratic in the fields 
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to the bare action. The purpose of the cutoff is twofold. First, it regularizes fermionic 
and bosonic infrared divergences when dealing with, for example, a massless Gold- 
stone boson or zero-temperature fermions whose propagator becomes unbounded in 
the vicinity of the Fermi surface. Second, R'^ incorporates the continuous flow param- 
eter along which the functional integration is performed: as a function of decreasing 
cutoff", interpolates smoothly between the bare action Fq for A = oo and the full 
eff"ective action F, recovered in the limit A ^ 0. 

To regularize and interpolate, the matrix elements of R-^ have to fulfill the follow- 
ing properties: 

m \im^^^ R\Q) = ^ . (3.12) 

Integrating e"^""'^'' in the presence of source fields coupling linearly to S and S 
yields the cutoff"-dependent generating functional for connected Green functions 

G\S',S'] = -log j z)^z)^e-^o[^'^]-«'^[^.^]+(^''^)+(^'^') . (3.13) 

The cutoff-dependent effective action F^ is defined as 

F^[S,S] = £G'^[S',S']-'R'^[S,S] , (3.14) 

where £,G'^ is the Legendre transform of G'^ as specified without a cutoff in Eq. (13.91 ). 

Executing a scale-derivative on both sides of Eq. (|3.14l) and using Eq. (13.111) . we 
obtain 

d^F^[S,S] = d^G^[S',S'] - iS,R^S) , (3.15) 
where R'^ = ^^R'^ and we further compute 

'dG^[S',S'] ^^5G^[^',^']\ „ .^d^G^[S',S'] 




dS' J dS'dS' 
{F^{S, S]+n''{S,S]y ' 

dSds ' ' ^ • 
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where the supertrace Str traces over all indices with a plus sign for bosons and a minus 
sign for fermions. For the first term in the last line, we have used Eq. (|3.10l) and for the 
second term in the last line we have used the relation d^G/dS'dS' = {d^r/dSdS) 
augmented with the cutoff term from Eq. (13.141) . Denoting 

T^'^'[S,S] = ^'['^ff\ (3.17) 

and assembling Eqs. (I3.16[|3TT5] ). we obtain the exact flow equation for the effective 
action ^ 

-^r\S,S] = Str— — ^ . (3.18) 

Note that the definitions of vary slightly in the literature. In particular, F^ is fre- 
quently defined as the Legendre transform of without subtracting the regulator 
term 'R^, which leads to a simple additional term in the flow equation. 

To expand the functional flow equation (13.181) in powers of the fields, we write the 
Hessian of F^ as 

r(2^^[.S,.S] = -(G^)'' + f(2^^[.S,.S] , (3.19) 

where t^^^^[S, S] comes from terms which are at least quadratic in the fields. Defining 
= [(G^)"' - R"^]"^ and expanding in powers of F^^^^ yields 

—F^ = -Str (R'^G^) - Str [G/ff ^^^^ + f ^^^^G^f ^^^^^ + ...)] , (3.20) 

where 

G'k^ = G^R^G^ . (3.21) 

For a sharp cutoff, Eq. (13.211) denotes the so-called single-scale propagator (Salmhofer 
2001, Metzner 2005) as it then has support only for energies = A due to a ^-function. 
In general, G^'^ is smoothened and may thus have support on more energies, typically 
in an energy window around energies « A. 

Expanding both sides of Eq. (13.201) in powers of the fields and comparing coef- 
ficients yields an (infinite) hierarchy of flow equations for the vertex functions. The 
zero-point vertex, y^")'^, see also the first line of Fig. [3711 equals the effective potential 
with the analytic flow equation from the first term in Eq. (13.201) : 

— C/^ = -Str(R^G^) . (3.22) 
dA 



Flow equations for thermodynamic properties of the system such as the specific heat 
or the entropy follow from derivatives with respect to the temperature of Eq. (13.221) . 
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3.2.2 Spontaneous symmetry-breaking 

The low energy regimes of the models described in Eqs. (13.11 13.21) are often domi- 
nated by phases with broken symmetry, that is, states where interaction effects spon- 
taneously break one, or multiple, of the symmetries present in the original bare action. 
The physical properties in such a phase are vastly different from the phase where the 
original symmetry is preserved. 

In terms of our superfields, the symmetry-breaking is associated with a non- 
vanishing expectation value of one of the bosonic components of S, referred to as 
the order parameter. While in the symmetric phase the flow equations are expanded 
around = 0, in the symmetry-broken regime we expand around a generally non-zero 
value of the q = component of the bosonic field, that is, we expand r^[ifr, + 0] 
in powers of if/, ijf, 0, (p*, where may be non-zero (see also Delamotte 2004, Schiitz 
2006). 

We now illustrate how this expansion works. For this purpose, we restrict ourselves 
to the case of a real- valued scalar field and consider only local terms of the effective 
action. The ansatz for the bosonic components of the left-hand-side of Eq. (13.201) is 

ncf>]= ff2^y'''i^-^y ' (3.23) 

where the integration symbol denotes the standard position space integration. Per- 
forming a scale-derivative on both sides of Eq. (13.231 ) yields 

^n\ ' {n-l)\ 

Comparing powers of - a'^) between Eq. (13.241 ) and the right-hand-side of Eq. 
(13.201) yields the following prescription for the flow equation for the n-point vertex 

y(n)A _ rA^^ 1-loop IPI diagrams generated by with n external legs) 

_ (3 25) 

The expansion point will be chosen such that the bosonic 1 -point function y*^^^^ = 
for all A: 



I \ I \ 

\ I \ I 

Y _ Y 

y^'""^ ' (3.26) 

In this way tadpole contributions are absorbed into the flow of a. The vanishing of 
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Fi^. Truncated hierarchy of 1 PI vertex functions up to quartic order in the fields 
for phases with spontaneously broken symmetry. The cross denotes the flow of the 
order parameter d^a^. Usually, the presence of vertices with an odd number of legs is 
peculiar to phases with broken symmetry. It is important to note that for the general 
case of mixed Fermi-Bose theories, many more diagrams with normal and anomalous 
propagators as well as normal and anomalous vertices appear, see Chapter^ 
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the bosonic one-point function is equivalent to the condition that the minimum of 
the effective potential coincides with its first field derivative for all A: U'{(p\\^=a = 
(Berges 2002). 

In a general polynomial truncation of Eq. (l3.23l) up to quartic order in the bosonic 
fields, the hierarchy for the n-point vertices can be represented in terms of Feynman 
diagrams, see Fig. (13.11) . 

One aspect of the practicability of the IPI-formulation of the functional RG 
lies in its intimate kinship to conventional perturbation theory (Negele 1987). Eq. 
(13.251) contains precisely the diagrammatic corrections that would show up in a IPI- 
expansion of perturbation theory. Technically, a computation of the flow equations 
therefore involves the determination of prefactors and signs of the various Feynman 
diagrams. One way to achieve this is to compare coefficients in Eq. (13.201) . However, 
in symmetry-broken phases, the Nambu structure of the fields leads to at least 4x4- 
matrices for and V^^^, thus complicating the trace operation on matrix products 
with several of these matrices, and necessitating systematic use of software packages, 
for example Mathematica. 

In this thesis, we determine the prefactors and signs by use of Wick's theorem 
and Feynman rules directly (see chapter 2.3, Negele 1987; chapter 4.4, Peskin 1995; 
chapter 6, Weinberg 2005). 
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In chapter [21 some scenarios going beyond the conventional Fermi liquid description 
for interacting Fermi systems were discussed. After recapitulating the Fermi liquid 
fixed point in the functional integral formalism in section I2.1[ prominent competing 
and conspiring ordering tendencies of two-dimensional lattice fermions were intro- 
duced. We have argued that the fermionic functional RG is a useful tool to treat these 
cases but for computations within symmetry-broken phases a description in terms of 
the collective bosonic degrees of freedom seems more practicable. In section [Z2l the 
mechanism of spontaneous symmetry-breaking is introduced in detail. We zoomed in 
on the phenomenology of Goldstone bosons which are emitted in interacting Fermi 
systems with broken continuous symmetry for example with magnetic or superfluid 
order. Critical points that separate phases with different symmetry properties were 
brought into play in section [231 The standard phase diagram in the vicinity of quan- 
tum critical point was exhibited in Fig. 12.21 and the conventional Hertz-Millis ap- 
proach to theoretically describe quantum critical behavior was outlined briefly. Es- 
sential renormalization group ideas and references are summarized in section [241 

In chapter [31 the functional RG framework was presented as a promising tool to 
transform the before-mentioned theoretical concepts into real computations for sys- 
tems containing both fermionic and bosonic degrees of freedom. We defined the effec- 
tive action, that is, the functional which generates the one-particle irreducible vertex 
functions for superfields in section I3.1.2[ The exact flow equation for the effective 
action is derived in section 13.2.11 The central result of this chapter is contained in 
subsection |3.2.2t the prescription how to include scale-dependent expectation values 
into the hierarchy of vertex functions to account for spontaneous symmetry-breaking 
is pictorially shown in Fig. 13.11 

In the following Part [III we merge the theoretical concepts with the functional 
RG framework and apply it in four cases. The chapters are sequenced in order of 
increasing truncation complexity, consecutively build on and extend the previous one, 
but may yet be read independently. 
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Applications 
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Hertz-Millis theory with discrete symmetry-breaking 



5.1 Introduction 

Quantum phase transitions in itinerant electron systems continue to ignite consider- 
able interest (Sachdev 1999, Belitz 2005, Loehneysen 2007). In many physical situa- 
tions, a line of finite temperature second order phase transitions in the phase diagram 
terminates at a quantum critical point at T = 0. In such cases quantum fluctuations 
influence the system also at finite temperatures altering physical quantities such as the 
shape of the phase boundary. Consequently, in a complete description of the system at 
finite r, quantum and thermal fluctuations have to be accounted for simultaneously. 

The conventional renormalization group (RG) approach to quantum criticality in 
itinerant electron systems (Hertz 1976, Millis 1993), see also section [231 relies on 
the assumption that it is sensible to integrate out fermionic degrees of freedom from 
the functional integral representation of the partition function and then to expand the 
resulting effective action in powers of the order parameter alone. This approach has 
been questioned for magnetic phase transitions associated with spontaneous breaking 
of continuous spin rotation invariance, since integrating the fermions leads to singular 
interactions of the order parameter field (Belitz 2005, Loehneysen 2007). 

In this chapter, we focus on quantum phase transitions to phases with broken 
discrete symmetry. We analyze quantum and classical fluctuations in the symmetry- 
broken phase with an Ising-like order parameter near a quantum critical point. Our 
calculations are based on a set of coupled flow equations obtained by approximating 
the exact flow equations of the one-particle irreducible version of the functional RG. 
Quantum and classical (thermal) fluctuations are treated on equal footing. The func- 
tional RG has been applied extensively to classical critical phenomena (Herges 2002), 
where it provides a unified description of 0(A^)- symmetric scalar models, including 
two-dimensional systems. The classical Ising universality class has been analyzed in 
Refs. 9-11. In our approach, we can compute the RG flow in any region of the phase 
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diagram, including the region governed by non-Gaussian critical fluctuations. This al- 
lows comparison between the true transition line and the Ginzburg line, which so far 
has been used as an estimate of the former (Loehneysen 2007, Millis 1993). Specif- 
ically, we capture the strong-coupling behavior emergent in the vicinity of the tran- 
sition line as well as the correct classical fixed point for Tc, including the anomalous 
dimension of the order parameter field. 

Our results may be applied to Pomeranchuk transitions (Metzner 2003, Dell' Anna 
2006, Wolfle 2007) which spontaneously break the discrete point group lattice sym- 
metry, and to the Quantum Ising model (Sachdev 1999). An example of a fermionic 
model displaying a genuine quantum critical point already on mean-field level is the 
so-called f-model, where forward- scattering processes tend to deform the Fermi sur- 
face leading to Pomeranchuk transitions (Yamase 2005, Dell' Anna 2006). The control 
parameter here is the density/chemical potential. 

In Section [5]2] we introduce Hertz's (1976) action and the effective action for the 
Quantum Ising model, adapted to the symmetry-broken phase, which serves as a start- 
ing point for the subsequent analysis. In Section [53] we describe the functional RG 
method and its application in the present context, and subsequently derive the RG 
flow equations. In Section l54l we present a solution of the theory in the case T = 0. 
Section l53] contains numerical results for the finite T phase diagram in the region with 
broken symmetry. Difl"erent cases are discussed, distinguished by the dimensionality 
d and the dynamical exponent z- For Pomeranchuk transitions z = 3 and for the Quan- 
tum Ising model z = 1 . In particular, we compare the T^. line with the Ginzburg line, 
thus providing an estimate of the critical region size. In Section 15.61 we summarize 
and discuss the results. 

5.2 Bosonic action 

The starting point of the standard RG approach to quantum critical phenomena in itin- 
erant electron systems is the Hertz action (Hertz 1976, Millis 1993). It can be derived 
from a microscopic Hamiltonian by applying a Hubbard-Stratonovich transformation 
to the path-integral representation of the partition function and subsequently integrat- 
ing out the fermionic degrees of freedom. The resulting action is then expanded in 
powers of the order parameter field, usually to quartic order. 

The validity of this expansion is dubious in several physically interesting cases, 
in particular for magnetic transitions with SU(2)-symmetry, since the integration over 
gapless fermionic modes can lead to singular eff"ective interactions of the order pa- 
rameter field, which may invalidate the conventional power counting (Belitz 2005, 
Loehneysen 2006, Abanov 2004). Such complications do probably not aff"ect tran- 
sitions in the charge channel and magnetic transitions with Ising symmetry. There 
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are several indications that singularities cancel in that case, namely the cancellation 
of singularities in effective interactions upon symmetrization of fermion loops (Neu- 
mayr 1998, Kopper 2001), and the cancellation of non-analyticities in susceptibilities 
(Rech 2006). 

We therefore rely on the usual expansion of the action to quartic order in the or- 
der parameter field and also on the conventional parametrization of momentum and 
frequency dependences, which leads to the Hertz action (Hertz 1976, Millis 1993) 

Sum = '^yJ^4>p(M^ + p')4>-p + urn . (5.1) 



2^ J i2nY'"'y\p\^-^ 

Here (f> is the scalar order parameter field and (pp with p = (p, a;„) its momentum rep- 
resentation; CO,, = 2nnT with integer n denotes the (bosonic) Matsubara frequencies. 
Momentum and energy units are chosen such that the prefactors in front of and 
are equal to unity. The Hertz action is regularized in the ultraviolet by restricting 
momenta to |p| < Aq. The value of the dynamical exponent is restricted to z > 2. The 
case z = 3 for Pomeranchuk transitions is of our main interest. Formally, the results 
obtained with Su[(f>] as our starting point may be applied to systems with arbitrary 
z>2.Q 

Other interesting cases include spin systems such as the Quantum Ising model for 
which z = 1 and the associated continuum action reads (Sachdev 1999) 

where the frequency term enters with the same power as the momentum term. 

In the symmetric phase the potential is minimal at = 0, and is usually 
parametrized by a positive quadratic and a positive quartic term (Hertz 1976, Millis 
1993). Since we approach the quantum critical point from the symmetry-broken re- 
gion of the phase diagram, we assume a potential U[^] with a minimum at a non-zero 
order parameter (pQ-. 



(5.3) 

where (p and cp' are functions of x and r with cp = cpo + (p'. The parameter 

^ A comprehensive discussion of the origin of the co„ dependence of the action is 
given by Millis ( 1993). Eq. d5.iD is valid in the limit j^pj «; 1 which is relevant here 
since the dominant fluctuations occur only in this regime. 



l/T 



dT I d'^x 
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5= — = ^, (5.4) 

controls the distance from criticality. Approaching the phase boundary in the (5, T)- 
plane from the symmetry-broken phase gives rise to the three-point vertex V3 u 5, 
which generates an anomalous dimension of the order parameter field already at one- 
loop level. 

Although formally correct as a result of integrating out the fermions, the Hertz 
action (|5.1I) is not a good starting point for symmetry -broken phases with a fermionic 
gap, such as charge density wave phases or antiferromagnets. A fermionic gap leads 
to a suppression of the dynamical term (linear in frequency) in the action, since it 
suppresses low energy particle-hole excitations. If not treated by a suitable resum- 
mation in the beginning, this eff"ect is hidden in high orders of perturbation theory 
(Rosch 2001). We do not deal with this complication in the present chapter as our 
main interest lies in systems with Pomeranchuk deformed Fermi surfaces where the 
particle-hole continuum remains gapless in the phase with broken symmetry. How- 
ever, our results for the transition temperature should not be affected by a gap in the 
symmetry -broken phase, since it vanishes continuously at Tc. 



5.3 Method 

To analyze the quantum field theory defined by S [0] we compute the flow of the effec- 
tive action r'^{(p\ with approximate flow equations derived from an exact functional 
RG flow equation (Wetterich 1993, Berges 2002, Delamotte 2004, Gies 2006). The 
exact flow equation and the vertex expansion including symmetry-breaking was de- 
rived and discussed in detail in chapter [3l The effective action r^{(p\ is the generating 
functional for one-particle irreducible vertex functions in presence of an infrared cut- 
off A. The latter is implemented by adding a regulator term of the form / ^(f>R^(f) to the 
bare action. The effective action interpolates smoothly between the bare action S [0] 
for large A and the full effective action r{(p\ in the limit A ^ (cutoff removed). Its 
flow is given by the exact functional equation (Wetterich 1993) 

where = BaR^, and T^^^ = 6^r^{(p\l6(p^ . In momentum representation {(pp), 
the trace sums over momenta and frequencies: Tr = T / For the regulator 
function i?^(p) we choose the optimized Litim cutoff (Litim 2001) 



i?>)=Zp(yl2-p2)0(A2-p2) , 



(5.6) 
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where Zp is a renormalization factor (see below). This regulator function i? (p) re- 
places ZpP^ with ZpA^ for |p| < A. Following (Berges 2002) we neglect Zp in R^(p), 
such that i?'^(p) = 2ZpA0(A^ - p^) . The scale-derivative of the Z-factor leads to ad- 
ditional terms on the right-hand- side of the flow equations proportional to the anoma- 
lous dimension rj defined in below in Eq. (15.171) . These terms are neglected here as 
they are of higher order in the vertices. 



5.3.1 Truncation 

The RG flow of the local potential of the form Eq. (15.31) will be followed by cutofi"- 
dependent parameters u and 0o (or, alternatively, 6). 

The quadratic part of the efl^ective action for the Hertz action contains two renor- 
malization factors, one for the frequency and one for the momentum dependence: 

1 . /„ \Ci>„\ „ 2 



r<l>,p,H = J :^(f>~p y^^-^^ + ^pP +^j(f'p' (5-7) 

where J^^ comprises frequency integration (Matsubara summation) and momentum 
integration at zero (finite) temperature. As it will turn out below, the frequency renor- 
malization for the Hertz action does not play a role, in contrast to the Quantum Ising 
model, where we have 

= \'P'P i^^^l + ■ZpP^ + ^) ^P • (5-8) 

The inverse propagator exhibits relativistic invariance and Z^ is equally important to 
Zp at zero temperature. 

The Green function for the Hertz action endowed with the regulator reads accord- 
ingly, 

Gr (p) = -{(p-n(p„) = ^ ^ , (5.9) 

' ' Z^|J^+Zpp2+5+i?A(p) ' 

and similarly for the Quantum Ising model with the frequency term replacement 
Zt^l^p r Z^o;^. In all expressions with the frequency term of the Hertz action, Eq. 
(15.71) . the corresponding terms for the Quantum Ising model are obtained by the just 
described replacement of frequency terms. 

We emphasize that we have used a relatively simple parametrization of the effec- 
tive action r'^\_(p\. In particular, we kept only the dominant terms in the derivative 
expansion (Berges 2002) and neglected the field dependence of the Z-factors. Fur- 
thermore, the simple parametrization of the eff'ective potential Eq. (15.31) allowed us 
to substitute the partial diff'erential equation (15.101) governing the flow of U{(p\ by 
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the two ordinary differential equations (15.321) and (|5.33l) . The latter approximation is 
equivalent to neglecting all higher order vertices generated during the flow. More so- 
phisticated truncations have been applied, among others, in the context of the classical 
Ising universality class, where they lead to improved results for the critical exponents 
(Ballhausen 2004). 



5.3.2 Flow equations 



We derive our flow equations keeping the notation general so as to account for zero 
and finite temperature. Evaluating Eq. (|5.5I) for a momentum-independent field (p 
yields the flow of the efi"ective potential U[(p\ 



2 ^-z^|^+Zpp2+i?^(p) + [/"[0] 



(5.10) 



from which we derive the flows of the parameters and u, following the procedure 



in (Berges 2002). Viewing t/ as a function of p = and using U'[pq\ = 0, we can 



write 



= 4jU'[po] = d^U'lpo] + U"[po] d^po . 
dA 



(5.11) 



Inserting dy\U'[po] as obtained by difi'erentiating Eq. (15.101) with respect top atp = po, 
and using U"[po] = one obtains the flow equation for po 



Z^I^+Zpp2+i?'^(p) + fMpO 



(5.12) 



The flow of u is obtained by diff"erentiating Eq. (|5.10l) twice with respect to p: 



Oau = 3u Tr 



Z,jJj^+Zpp2+i?^(p) + |i.po 



(5.13) 



Inserting the above flow equations into the A-derivative of 6 = |mpo, we obtain the 
flow of 6. To complete the system of flow equations one still needs to derive the 
evolution of Zp and Z^, which parametrize the momentum and frequency dependence 
of the propagator. Taking the second functional derivative of Eq. (15.51 ), we obtain the 
flow equation for the propagator 



dAG'^\p) = 3u6Tr [R^ Gl(q) G^iq + p)] 



(5.14) 



where Gr(^) is given by Eq. (15.91) . Here we skipped the contribution from the tadpole 
diagram, since it involves no dependence on momentum and frequency, and therefore 
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does not contribute to the flow of the Z-factors. The momentum renormalization factor 
is then given by 

^ ^ ^ ■ ^ , (5.15) 



Zp= —dl[Gl\v,o^n = 0)\ 



where is the Laplace operator evaluated at constant cutoff" function, that is, does 



not act on i?^(p). We now turn 
We first consider the case z = 2, where 
at finite temperatures can be related 
to the propagator by 



to 



p=0 

cuto 

the renormalization 



of Z„. 



Z.,= 



1 



G«Hp = 0,27rr) 



InT 
-G^i(p = 0,0) 



(5.16) 



But we shall later show that Z^^ renor- 
malizes only very weakly at finite tem- 
peratures for z = 2 and therefore even 
more weakly for Pomeranchuk transi- 
tions where z = 3. At T = 0, the 
derivatives become continuous and for 
the Quantum Ising model z = 1, Z^ 
renormalizes as strongly as Zp. In the fol- 
lowing, we will often employ the anoma- 
lous momentum scaling exponent de- 
fined as: 



7] = 



d log Zp 
J log A 



(5.17) 



and similarly for t/^j. The contributions 
to the flow of 6, u, Z^ and Zp are illus- 
trated in terms of Feynman diagrams in 
Fig.EU 




u : 




7 7 



Fig. 5.1. Feynman diagrams representing 
the contributions to the flow equations 
BJ2\\5J3\\5l4\) . 



5.4 Zero-temperature solution at the quantum critical point 

In this section we present a solution of the flow equations (15.121 15.13[ 15.141 ) at zero 
temperature. For z > 2, we can linearize the flow equations around the Gaussian fixed 
in J = 2 and we provide an analytic expression for the value of the control parameter 
^0 corresponding to the quantum critical point. 

For z = 1, the two-dimensional QCP is described by a non-Gaussian fixed point 
and we resort to a numerical solution of our equations. 
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5.4.1 z > 2 

Convenient handling of the flow equations is achieved by introducing the variables 

A, 



6 = 



(5.18) 



For the Hertz action, the resulting flow equations are 



d6 



J log A 

dv 
J log A 



(77-2)5 + 



U6v 



1 



= (J + z - 4 + 2t])v + 



(l+6yd + z-2 



+ 



8v 



1 



l+6d+z-2' 
1 



(1 + S)2 J + z-2 ' 



6v 



d(l +6Y 



4 ~ (d- 2)(z -2) ~ , 4 
3 3{d + z-A) d + z 



(5.19) 



The flow equations (15.191) reveal that v is an irrelevant variable for J+z > 4 as pointed 
out earlier by Hertz (Hertz 1976). In this case Eqs. (15.191 ) have a stable Gaussian fixed 
point in V = 0, 5 = 0, 7/ = 0. Linearizing the flow equations around the fixed point, 
one obtains the solution 



5(A) = 



6 + 



8v 



(J + z-2)2 
r](A) = , 



{(A/Ao)'^'~^ - 1) 



(A/Ao)- 



(5.20) 



where 6 and v on the right hand sides are the initial values of the parameters at A = Aq. 
In the marginal case d + z = 4- one finds logarithmic convergence of v and 6 to zero. 
Expressing the order parameter (po in terms of 6 and v, substituting the above solution, 
and taking the limit A — > yields 



00 OC ^/S-So , 



where 



2A„ 



a: 



d+z-2 




(5.21) 



(5.22) 



n {d + z-2f 

is the quantum critical point's coordinate. From Eq. (15.211) we read ofl^ the value of the 
exponent /3 = ^, consistent with mean-field theory. From Eq. (15.201) one also straight- 
forwardly evaluates the correlation length ^, using ^"^ = lim^i^o 5{A) = lim^^o ZpA^d, 
which yields 
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^ = {5-6,)-'l\ (5.23) 

as expected within mean field theory. 

We have recovered the well-known fact that quantum phase transitions have prop- 
erties similar to their classical counterparts in effective dimensionality D = d + z 
(Sachdev 1999). In the case studied here D > A and one obtains mean-field behavior 
governed by a Gaussian fixed point. 



5.4.2 z = 1 

We now deal with the case where D = J -I- z is smaller than four and thus even at 
zero temperature the QCP falls into a non-Gaussian universality class for \ < d < 
3. Although this case is covered by the quantum-classical mapping and the zero- 
temperature quantum theory matches that of the classical theory in J -I- 1, it is worthfile 
to inspect the zero-temperature flow equations in hindsight to a subsequent compari- 
son with the finite-temperature equations in the following chapter. 

We put Z^^ = Zq as there is nothing that breaks the relativistic invariance except for 
our regulator function which cuts off" momenta but leaves the frequencies untouched. 
In addition to rj defined in Eq. (15.171) . we here employ the rescaled variables 

pdZn 

P = 



The frequency integrations can be performed analytically yielding the flow equations 

dp /i J N ~ 3 1 
-(l-d-r])p + 



dlogA 24(i + 2|p^3/2 



du T / 3 . 

= (J-3 + 2?7)m + 3mM ' (5-25) 



dlogA \^i6(l + 2M)5/2 

and for the anomalous dimension: 



.2.135V3(3 + 12J + 8 (J + 2) iip) 
128 (J + 2) (3 + 2upf 



We solve these equations numerically at the fixed point for 1 < d <3. The anomalous 
dimension comes out as 77 = 0.126 for d = 2. This can be compared with the best 
known estimates for the three-dimensional classical 0(\)-mode\, where 77 « 0.04 from 
seven-loop perturbation theory and e-expansions (Berges 2002). It is a general feature 
of functional RG computations in simple truncations that the anomalous dimension 
is overestimated. Improving the truncation by including a field-dependent eff"ective 
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potential and higher-order terms in the derivative expansion leads to accurate critical 
exponents (Berges 2002). 

A thorough discussion and flows versus the cutoff'-parameter follow in subsection 
[5331 

5.5 Finite temperatures 

Temperature provides for an additional scale in our problem that significantly enriches 
the analysis. There are now two relevant parameters, the first being the control param- 
eter and the second the temperature. In the {6, r)-plane there is now line of fixed points 
ending at the QCP (cJcnt, 0). Notably, for any, however small T , the critical behavior at 
A ^ is classical. Yet, we shall see below that the shape of the Tc-line is dominated 
by quantum fluctuations. 

At finite temperatures, the continuous frequency integrations are substituted by 
discrete Matsubara sums. These sums in the flow equations (15.121) and (15.131 ) can be 
expressed in terms of polygamma functions W„{x), defined recursively by Wn+\{x) = 
¥^{x) for n = 0, 1, 2, ... , and ¥o(x) = r'{x)/r{x), where r{x) is the gamma function. 
From the Weierstrass representation, r{xy^ = xe'^^ XX^=\^^ + where y is the 

Euler constant, one can derive the relation 



Taking derivatives with respect to x yields expressions for sums involving higher neg- 
ative powers of -I- x) in terms of polygamma functions of higher order. The d- 
dimensional momentum integrals on the right hand side of the flow equations can be 
reduced to one-dimensional integrals, since the integrands depend only on the modu- 
lus of p. 

Explicit dependences on A, Z-factors, and lengthy numerical prefactors in the flow 
equations can be eliminated by using the following rescaled variables: 




(5.27) 



P = Ipl/A , 
~ ^ 2;rZ^ , 

" ZpA^ 



(5.28) 



(5.29) 




(5.30) 



AdT 



(5.31) 



u , 



withA^ = {2n) ''Srf-i, where S^-i = 2;r'^^^/r(J/2) is the area of the ((i-l)-dimensional 
unit sphere. 
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5.5.1 z = 3 

For the case z = 3, the flow equations for 6 and u are obtained as 



d6 



J log A 



(?7 - 2) § + 4m 



1 1 



+ 



1 



+ \2u6 



1 1 



1 



d{\+6y Jo 



(5.32) 



dU 



J log A 



{d-4 + 2T])u+ \2u^ 



1 1 



1 , 

- / dpp 



d{l+6f Jo 



d+3z-l 



(5.33) 



where x = I + f-\\ + 6)p^-'^. Inserting Eqs. (I5TT41) and (l5T5l) into Eq. (ISTTl) . and 
performing the frequency sum, one obtains for the anomalous dimension 



u6 



d{\+6f 



2(1 + ~6)- 



8 



J + 2. 



T 



-5 



- 6p\z - 2){d + z- 4)f 2 W2{x) - 2p^^' (2(8 - d) 
{\ + ~5-p^) + [{d - 14)(1 + ~5) + 8p^]z + 3(1 + ~5)z^ 



f¥,(x) - p^'[2ip^ - 1) + ~6iz - 2) + z?¥,ix) 



(5.34) 



Taking the logarithmic derivative with respect to A, inserting Eq. (|5.14l) . and perform- 
ing the trace yields for rj^ for the case z = 2 



12 ~~ 
d 



T -{\+5) 
f (1 + ~6f 



+ T-'WA 1 + 



1+6 



(5.35) 



From the numerical solution of the flow equation we observe that rj^a is small at all 
scales and vanishes for A ^ 0. For example, in two dimensions /y^, varies between 
-0.033 and 0.005 for w = 1 at T = e and has practically no influence on the phase 
diagram. In three dimensions the values are at least one order of magnitude lower. 
For z > 2 one expects an even smaller 77^, since a larger z reduces the strength of 
fluctuations near the quantum critical point. Therefore we set t]^) = and Z^j = I from 
now on. The scaling variable T then obeys the flow equation 



dT 



^^logA 



= {ri-z)T . 



(5.36) 



In the flow equations (15.3215.3315.341) one identifies the classical mean-field (involv- 
ing only one power of 6 ox u), classical non-Gaussian and quantum (involving f) 
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terms. The quantum contributions vanish as the infrared cutoff tends to zero at con- 
stant non-zero temperature {t^^ «; 1). On the other hand the quantum terms dominate 
the high energy part of the flow, where » 1 . Our framework allows for a continu- 
ous connection of these two regimes of the flow. The scale Ad below which quantum 
fluctuations become irrelevant depends on temperature. It vanishes at T = 0. Li the 
absence of a sizable anomalous dimension rj in the quantum regime of the flow, one 
has 

Ad oc T'l' , (5.37) 

as follows directly the definition of T. It turns out that 77 is indeed negligible down to 
the scale Ad, except in the case z = 2 in two dimensions (see Sec. 15.51) . 

Below, we numerically solve the full RG flow equations (15.3215.3315.341) for fi- 
nite temperatures. As already announced, the analysis is focused on the region of the 
phase diagram where symmetry-breaking occurs. In our notation this corresponds to 
sufficiently large values of the control parameter 6. First we treat the case z = 3 (in 
d = 2,3) to which Eqs. (15.3215.3315.341) are applied directly. In all numerical results 
we choose an initial cutoff Aq = 1, and an initial coupling constant u = \. 

We first solve the coupled flow equations (|5. 3215. 3315. 341) with the aim of deter- 
mining the phase boundary Tc{5), or, equivalently 6c{T). To this end, for each given 
temperature we tune the initial value of 6 such that at the end of the flow (for A ^ 0) 
one obtains the critical state with 6(A) 0. The tuned initial value is then identified 
as 6ciT). Inverting the function SdT) yields T^S). In the variables 6, u this corre- 
sponds to seeking for such values of the initial 6, that both 6{A) and m(A) reach a fixed 
point as the cutoff is removed. 

The flow of rj and a as a function of the logarithmic scale variable s = - log(A/Ao) 
is shown in two exemplary plots in Fig. 15. 21 respectively. The flow is shown for various 
temperatures T , with 6 tuned to values very close to 6c{T). The plateaus in Fig. l5.2l cor- 
respond to non-Gaussian fixed point values of the flowing parameters. From Fig. 15.21 
where d = 2, one reads off the value of the anomalous dimension 77 ^ 0.5. The exact 
value from the Onsager solution to the Ising model is | . For the case J = 3 we find 
ri X 0.1 within our truncation, which is also about twice as large as the best estimates 
for the exact value (Goldenfeld 1992). To obtain scaling behavior in the range of a 
few orders of magnitude one needs to fine-tune the initial conditions with an accuracy 
of around 15 digits. The breakdown of scaling behavior observed for very small A is 
due to insuflicient accuracy of the initial value of 6 and numerical errors. The plateaus 
are more extended as we go on fine-tuning the initial condition. Only exactly at the 
critical point true scale invariance manifested by plateaus of infinite size is expected. 

The figures also reveal the Ginzburg scale Ac at which non-Gaussian fluctuations 
become dominant, such that the exponent 77 attains a non-zero value. By fitting a 
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2 4 6 S 10 12 14 16 5 10 15 20 



Fig. 5.2. Left: Anomalous dimension ri(A) plotted as function of s = -log(AMo) 
for different values of temperature in the case z = 3 and d = 2. The function 77(5) 
exhibits crossover from the mean-field value t] = 0, to the non-Gaussian result rj ^ 0.5. 
The crossover scale Ac is shifted towards smaller values (larger s) as temperature is 
reduced. Right: Quartic coupling u{A) plotted as function of s = - log(A/Ao) /or 
different values of temperature in the case z = 3 and d = 3. 



power-law we observe Ac oc for J = 3 and Ac oc ^/T^ for J = 2 with non-universal 
proportionality factors. As expected, Ac vanishes at the quantum critical point, be- 
cause at r = the effective dimensionality D = d -\- z is above the upper critical 
dimension dc = 4. At finite temperatures, Ac is the scale at which u is promoted from 
initially small values (of order T) to values of order one. From the linearized flow 
equations one obtains 

Ag oc r," (5.38) 

in agreement with the numerical results for d = 2 and d = 3. Note that Ac Ad, 
since Ad oc T^^^ for z = 3, see Eq. (15.371) . Hence anomalous scaling (finite 77) is 
indeed absent in the regime where quantum fluctuations contribute, and non-Gaussian 
fluctuations appear only in the classical regime. As already mentioned, the quantum 
contributions influence the flow only at relatively large A for T > 0. In particular, 
they do not alter any fixed point values. However, at the beginning of the flow they 
dominate over the classical part and therefore are crucial for a correct computation of 
the initial value of 6 leading to a scaling solution in the infrared limit. 

Results for the transition line T^d) for J = 2 and J = 3 are shown in Fig. 15. 3[ In 
both cases we recover the shape of the transition line as derived by Millis (1993), who 
used the Ginzburg temperature in the symmetric phase as an estimate for . Namely, 
we find 

{6-6o)<^TAogT, (5.39) 

for d = 2, and 
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Fig. 5.3. The transition lines Tc{5) obtained for z = 3, d = 2 (a) and z = 3, d = 3 
(b). The phase with broken symmetry is located below the line. The phase boundary 
Tc(6) obeys {6 - 6o) oc T^ log T^for d = 2 and T^ oc (S - ^o)"'^^ in agreement with the 
result by Millis (1993). The crosses indicate the Ginzburg temperature Tofor various 
choices of 6. 



T,<x(6- 6of (5.40) 

for d = 3. The exponent 3/4 in the three-dimensional case matches with the general 
formula for the shift exponent 

'^=JT^ (5.41) 
d + z-2 

for arbitrary z in dimensions J > 2 as long as d + z > 4. Generally, we see from 
this formula that increasing z has qualitatively (neglecting possible log-corrections) 
the same effect as decreasing d. Note that the phase boundary Tc{6) approaches the 
quantum critical point with vanishing first derivative for d = 2 and with singular first 
derivative in the case d = 3. 



An advantage of the present approach is that one can also follow the RG flow 
into the strong coupling regime, where non-Gaussian critical behavior occurs. This 
in turn allows an estimate of the critical region's size as a function of temperature or 
the control parameter 6. To evaluate the Ginzburg line in the symmetry-broken phase 
one solves the flow equations (15.3215.3315.341 ) for fixed T and at diff"erent values of 
6 > 5c{T), observing the behavior of fixed point values of the average order parameter 
(po (or, alternatively, the correlation length ^) as 6 approaches 5c. Typical results are 
plotted in Fig. 15. 41 from which we read off" the value of the exponent /3 describing the 
decay of the order parameter upon approaching the transition line (pQ oc {6 - 5cf. In 
the truly critical region (for 5 - 5c small enough) one obtains /3 « 0.11 for d = 2 and 
/3 « 0.30 for d = 3. These results come out close to the correct classical values 0.125 
and 0.31, respectively. This is unlike the other critical exponents {rj and the correlation 
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Fig. 5.4. Order parameter 0o a function of {6 — 5c) atT = e~^ for z = 3, d = 2 (a), 
and z = 3, d = 3 {b). The exponent /3 governing the decay ofcpQ upon approaching 6c 
exhibits a crossover from a mean-field value yS = 0.5 to a non-Gaussian ~ 0.1 1 for 
d = 2, and /3 ~ 0.30 for d = 3. 



length exponent v) which within our truncation differ by factors close to 2 from their 
correct values. 

Indeed, as discussed in (Ballhausen 2004), to obtain accurate values of the critical 
exponents in the Ising universality class, and in particular in J = 2, one not only needs 
to consider the full partial differential equation governing the RG flow of the effective 
potential U{(p\, but also the field dependence of the wave function renormalization 
and higher orders in the derivative expansion of the effective action (Berges 2002). 

From Fig. 15. 41 we can extract the Ginzburg value 6g below which true critical be- 
havior is found at the chosen temperature T . Around 6g, the exponent (3 exhibits a 
crossover from its mean-field value [3 = 0.5 to a non-Gaussian value. In other words, 
6g marks the boundary of the non-Gaussian critical region at a given T. At zero tem- 
perature, 6g coincides with the quantum phase transition point Sq, since there the fluc- 
tuations are effectively d-\-z > 4- dimensional, leading to mean-field behavior. Several 
Ginzburg points in the 6 - T plane are plotted as Tg{6) in Fig. 15. 3[ where they can 
be compared to the phase transition line. In three dimensions Tg and almost coin- 
cide, such that Tg provides an accurate estimate for T,.. In two dimensions a sizable 
region between Tg and Tc appears in the phase diagram. In that region non-Gaussian 
classical fluctuations are present. 

We stress that accounting for the anomalous exponent 77 is necessary to describe 
the classical scaling regime, that is to obtain the plateaus in Fig. 15.21 Upon putting 
77 = the scaling plateaus do not form. On the other hand, the shapes of the phase 
boundaries and the Ginzburg curves become very similar at T ^ in the present 
cases, where d + z>'4. 
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Fig. 5.5. The transition lines T(.{6) obtained for z = 2, d = 2 (a) and z = 2, d = 3 
(b). The phase with broken symmetry is located below the line. For d = 2 the phase 
boundary Tc{5) is consistent with the relation {6 - 6q) oc T^ log log Td log T^ derived 
by Millis (1993). For d = 3 it follows the expected power law Tc{6) oc (5 - Sq)"^ with 
ij/ = 213. The crosses indicate the Ginzburg temperature Tcfor various choices of 6. 



5.5.2 z = 2 

In the case z = 2 the flow equations (I5.32I5.33I5.34I ) are significantly simplified, as 
all the integrals can be evaluated analytically. One obtains 



d6 



J log A 

du 
J log A 



= (77 - 2) § + 2a 



2 14 

d{\ +6Y d 



+ 3u6 



4 14 

- — -^--T-'W^iy) 
d{l+6f d 



= (d-4 + 27])n + 3it 



1 



d{l+6f d 



f-'W2iy) 



and for the anomalous dimension, 



^7 = 3 



d(l +6y 



2(1+5) 



d + 2 



4u6. 



fW,{y) + 



d + 2 



(5.42) 



(5.43) 



where the argument of the poly gamma functions is given hy y = 1 + (1 +6)/T. 

The procedure to evaluate the phase diagram and the Ginzburg line is the same as 
in the previously discussed case z = 3. In Fig. 15.51 we show results for the transition 
line Tc{6) in two and three dimensions. We also show the Ginzburg temperature Tq 
for various choices of 6. 

In two dimensions, the transition line is consistent with the almost linear behav- 
ior (6 - 60) oc Tc log log TJ log Tc , derived previously for the Ginzburg temperature 
(Millis 1993). However, a sizable region with non-Gaussian fluctuations opens be- 
tween and Tq. In three dimensions, T^d) obeys the expected (Millis 1993) power 
law Tc(6) oc (S - with shift exponent i/r = 2/3, and the Ginzburg temperature is 
very close to Tc for any 6. 
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5.5.3 z = 1 



We now consider the Quantum Ising model at finite temperatures. Tliis case lias not 
been considered before in d < 3 and also not by Millis (1993). The Millis analysis 
relies on the fact that the QCP is described by a Gaussian fixed point which is not 
applicable here. The zero-temperature equations (15.251 [5?26l) are generalized to finite 





Adp 


p = 


IKdTA^-^ 




IKdTu 


u = 






InT 


T = 






AA ' 



(5.44) 



to write the flow equations as 



dp 3 
—lL_ {2-d-ii)p+- 
(2 log A 2 



_(i+f)^'''5((„f)^i+f)' 



du 



J log A 



= (J - 4 + 2?7) ii + 3m 



.-,2 



(1 + f) 



+ 



3 / 3 

n=i ( [nff + 1 + 



(5.45) 



The anomalous dimension is determined by 
1 2 



TJ = 2u p 



L(l + 2|e)' (j + 2)(l + 2M) 



+ 



1 



(d + 2) (^{nTy + I + ^ 



and the effective temperature scales as 

df 



dlogA 



(5.46) 



(5.47) 



The Matsubara summation can be performed analytically, which we have done for a 
computationally cheaper numerical solution, but they do not deliver analytic insights 
here. 
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We discuss the numerical solution 
of the finite temperature flow equations 
(|5.45[ 15.461) and we compare it to the 
zero-temperature flow equations of sub- 
section I5.4.2[ Note that now d + z = 
3 < 4, and this case has not been 
covered previously (Hertz 1976, Millis 
1993). Non-Gaussian fluctuations persist 
even at zero temperature which is illus- 
trated in a qualitative comparison of the 
Ginzburg line with the Tc-line in Fig. 

m 

In Fig. 15. 7[ the zero-temperature flow 
is juxtaposed with finite-temperature 
flows. Clearly, there is no continuous 
crossover from the finite-temperature 
behavior to the zero-temperature one. 
Hence, taking the limit T is a 
discontinuous process. Both, the finite- 
T and the T = find a description in 
terms of two distinct non-Gaussian fixed 




QCP 



Fig. 5.6. Schematic comparison of the 
Ginzburg line with the Tc-line for d = 2, 
z = 1. The critical region extends down to 
zero temperature. 



pomts. I From Fig. 0(b), we deduce the Ginzburg-scale to vary with temperature as 



An oc T. 



1/2 



(5.48) 



fitting the formula derived for <i -I- z > 4, Ac T^J''^~'^\ As expected, the QCP- value 
for 7] shown in Fig. 15.81 (a) is smaller than its classical counterpart as the QCP is 
eff"ectively d + 1 -dimensional thus taming critical fluctuations. The phase boundary in 
double-logarithmic plot is shown in Fig. 15.81 (b). The points do not lie on a straight 
line, indicating log-corrections correcting the power-law behavior, i/r attains effective 
values in the range: i/r x 0.5 - 0.7 over six (!) orders of magnitude in good agreement 
with 



9.28- 



|ln(5-5.; 



(5.49) 



which, as anticipated, does not fit the Millis formula, Eq. (15.411) . which was derived 
under the assumption d + z > 4. However, the existence of these log-corrections 
remains inconclusive. Interestingly, i/r is reduced compared to its two-dimensional 



^ As a counterexample, consider the merger of the non-Gaussian Wilson-Fisher-type 
fixed point with the Gaussian fixed point when taking the limit d d^, with d^ the 
upper critical dimension, in standard (p^ -theory (Goldenfeld 1992). 
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(a) (b) 




Fig. 5.7. Flows for the Quantum Ising model (z = I) in d = 2. (a): Flows of p as a 
function of logarithmic cutoff-scale s = - log [A/Aq] for various temperatures and at 
T = 0. The values of the classical fixed point (p = 0.385j and the QCP (p = 0.197j 
are marked on the vertical axis. We set Aq = I. The infrared (ultraviolet) is to the 
right (left) of the graphs, (b): Corresponding flows ofu. The values of the classical 
fixed point (u = 6.87j and the QCP (u = 3.11) are marked on the vertical axis. 

analogs when z > 1 although the strength of non-Gaussian classical fluctuations is 
increased when reducing z. This underpins the importance of Gaussian quantum fluc- 
tuations which on the other hand decrease when reducing z- 




Fig. 5.8. Results for the Quantum Ising model (z = \) in d = 2. (a): Cutoff-scale 
dependences of the anomalous dimension for the same temperatures as in Fig. 15.71 
The value of the anomalous dimension at the classical fixed point (rj = 0.51 1 ) and the 
QCP (t] = 0.126j are marked on the vertical axis, (b): Double logarithmic plot of the 
phase boundary fitted to Eg. d5.49P . 
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5.6 Conclusion 

We have analyzed classical and quantum fluctuations in the symmetry-broken phase 
near a quantum phase transition in an itinerant electron system and in the Quantum 
Ising model. The analysis is restricted to the case of discrete symmetry breaking, 
where no Goldstone modes appear. Following Hertz (1976) and Millis (1993), we use 
an eff'ective bosonic action for the order parameter fluctuations as a starting point. 
For the Quantum Ising model, the effective continuum action is also of the ^'^-type 
but with dynamical critical exponent z = 1. The renormalization of the Hertz action 
and the Quantum Ising model by fluctuations is obtained from a system of coupled 
flow equations, which are derived as an approximation to the exact flow equation for 
one-particle irreducible vertex functions in the functional RG framework. In addi- 
tion to the renormalization of the effective mass and the four-point coupling, we also 
take the anomalous dimension of the order parameter fields into account. In the 
symmetry-broken phase, contributions to 77 appear already at one-loop level. Quan- 
tum and thermal fluctuations are captured on equal footing. The flow equations are 
applicable also in the immediate vicinity of the transition line at finite temperature, 
where fluctuations deviate strongly from Gaussian behavior. 

We have computed the transition temperature Tc as a function of the control pa- 
rameter 6 near the quantum critical point, approaching the transition line from the 
symmetry-broken phase. Explicit results were presented for dynamical exponents 
z = 3 and z = 1 in two and three dimensions corresponding to Pomeranchuck tran- 
sitions and the Quantum Ising model, respectively. For z = 3, T^id) agrees with the 
behavior of the Ginzburg temperature above T^. derived previously by Millis (1993). 
For z = 1 in two dimensions, even the zero-temperature theory is characterized by 
a non-Gaussian fixed point. This case is not covered in the Hertz-Millis theory. We 
computed the phase boundary and found logarithmic corrections to power-law be- 
havior. The effective shift-exponent, (^(J = 2,z = 1) « 0.5 - 0.7 is smaller than the 
two-dimensional Hertz-Millis value, ilf(d = 2,z>2) = l. 

We have also computed the Ginzburg temperature Tg below T^, above which non- 
Gaussian fluctuations become important. Although and almost coincide in three 
dimensions, a sizable region between Tq and opens in two dimensions. While for 
z > 2 this critical region shrinks to zero for ^ 0, a finite difference between the 
Ginzburg and critical coordinates persists even at zero-temperature for z = 1. 

It will be interesting to extend the present approach to the case of continuous sym- 
metry breaking. Then, symmetry-breaking at finite temperature is suppressed com- 
pletely for J < 2 by Goldstone fluctuations. We consider the fluctuation effects of 
Goldstone modes on the quantum critical behavior at zero and at finite temperature in 
the next chapter. 
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6.1 Introduction 

A quantum critical point may separate a phase with broken continuous symmetry 
from a normal phase (see Fig. l6.1l) as is the case for example in quantum magnets, su- 
perfluids (Belitz 2005, Loehneysen 2007), or continuum systems with Pomeranchuk 
Fermi surface deformations (Oganesyan 
2001, Wolfle 2007, Quintanilla 2008). 
The distinct signature of phases with 
broken continuous symmetry is the oc- 
curence of Goldstone modes. In a re- 
cent experiment on quantum magnets 
with pressure as non-thermal control pa- 
rameter (Riiegg 2008), inelastic neutron 
scattering data reconfirmed the vanish- 
ing mass for the Goldstone mode in the 
magnetic phase and provided informa- 
tion about the magnitude of the longi- 
tudinal mass in the vicinity of the QCP. Fig. 6.1. Phase diagram of the QCP with 
Theoretically, it is largely unexplored Goldstone modes for d > 2. At d = 2 ther- 
how Goldstone modes affect the zero and mal fluctuations destroy the order and sup- 
finite temperature properties on the or- press the T^-line to zero temperature ac- 
dered side of the quantum phase transi- cording to the Mermin-Wagner theorem. 
tion (Wetterich 2008). 

For classical bosonic 0(N)-models with Goldstone and longitudinal fluctuations, 
the functional RG provided a unifying picture (Berges 2002). In the present chapter, 
we compute RG flows in the vicinity of a QCP accounting for transversal and longi- 
tudinal fluctuations of the order parameter restricting ourselves to Goldstone modes 
dispersing linearly in momentum {Q ~ |q|) as is the case in superfluids and antiferro- 
magnets (Belitz 2005). 



Symmetric 



Broken continuous 
symmetry 
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At finite temperatures, we extend our computation of the phase boundary for dis- 
crete symmetry -breaking (see chapter [51 Jakubczyk 2008) to account for Goldstone 
modes. We compute the shift-exponent ij/ which characterizes how the Tf-line varies 
as a function of control parameter (Loehneysen 2007) in three dimensions. We also 
clarify the singularity structure of the theory and establish the infrared asymptotics 
away from the phase boundary. 

We start by presenting the cr-II model for phases with broken continuous sym- 
metry in section 16.21 We clearly distinguish longitudinal fluctuations (cr) from the 
Goldstone modes (11) both dispersing relativistically with dynamical exponent z = 1. 
In section [631 we describe our functional RG set-up. In section [641 we compute the 
shift-exponent characterizing the shape of the phase boundary at finite temperatures 
for J = 3. In section [631 we extend our calculation to regions away from the phase 
boundary and we analytically establish the infrared properties of both, the Goldstone 
and longitudinal propagator, as well as the bosonic self-interaction. Summary and 
conclusions follow in section [6!6l 

6.2 cr-n Model for continuous symmetry-breaking 

We are interested in the effective action for A^-component real-valued bosonic fields O 
in phases with spontaneously broken continuous symmetry {N > 2). The most general 
action to fourth order in the fields that respects the 0(A^)-symmetry and contains no 
more than two derivatives reads (Wetterich 1991, Tetradis 1994): 



where Z and Y are renormalization factors, d denotes the standard derivative operator, 
and the integration symbol / denotes the standard space integration. Including the 
field component index a and the dimensionality index /i, the products in the derivative 
terms are defined as: {d<^f = df^^^'d^^a, and (050)^ = {'P"d^<Pa) {<P''&'<Pb) where 
a,b = 1, N and fi = I, d+l. Local interaction terms are retained in the effective 
potential: 



where the effective potential has a minimum at O = (a, 0, 0...) with a real and pos- 
itive. It is convenient to decompose O into a londitudinal component cr and (A^ - 1) 
transverse components 11 via 




(6.1) 




(6.2) 



O = ((r,0,0,...) + (0,n) . 



(6.3) 
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Substituting O ^ (a + cr, 11) in Eq. (|6.1|) yields various interaction terms (see below) 
and two distinct quadratic parts: 

r„, = ^ju^,{zq^) n, , (6.4) 

where ^ = (^o, q) collects the frequency and momenta and = J ^ J com- 
prises the corresponding integrations. The components of the Il-field are mass- 
less as these are the Goldstone modes. The propagators for the cr- and Il-field 
have the form 



GA^, q) = -{cr^cr ) = ——^ — - 

Z^q^ + ml 

-1 Ja-n^ = 

(6.5) 

with p = lap, Z^. = Z, and = 
Z^ = Z + Yp 

m^. = up . (6.6) 
For the interaction terms, we obtain: ^^^^ ~ 



< 



■ TaV = 

/ \ / 
/ \ / 



X 



= / ya-*0--q-pO'-q'+pO'tj'0-q , 

'qA'^p 



GAk) = 



'q^q'^p 



q 



q^q ,p 



Fig. 6.2. Interaction vertices and propa- 



/rig. o.^. inieracHon ve races ana propa- 
7cr3 cr -q-pcr pcr q , gators of the cr-II model for continuous 

symmetry-breaking. 



^0-^2 = / 7o-7T^ O-.^-pIIpII^ , 

hP 

(6.7) 



with = Tr* = (" + Yp^) /8, y^2„2 = (u + Yp^) /4, and 7^3 = 7, 
(m -I- 7/?^) a/2 as depicted in Fig. 16. 2[ In total, the cr-II model contains four in 
dependent parameters: a, u, Z, and Y. 
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Fig. 6.3. Linearly realized 0(N)-symmetry. Here, if/Q = a is the order-parameter and 
the trough is understood to be N - I dimensional. From Weichman (1988). 

A comment on the choice of field basis in Eq. (16.31) is in place. The rotational 0{N)- 
symmetry in internal space is realized only linearly. The periodicity of the Goldstone 
field is deformed to a linear trough of infinite extension, see. Fig. 16.31 Consequently, 
periodic field configurations such as vortices are hard to capture in this basis. 

Alternative choices are the phase-amplitude representation (schematic): <P ~ 
Po^'^®, with S the - 1 group generators, the conjugate Goldstone modes, and po 
the amplitude of the order-parameter. This phase-amplitude representation manifestly 
preserves the rotational symmetry and allows, for example, a convenient description 
of the non-linear sigma model when the mass of the radial boson is sufficiently large 
so that Goldstone bosons remain the only fluctuating fields (Amit 2005). An analytic 
treatment of some aspects of the Kosterlitz-Thouless phase for N = 2 is also easily 
accessible (Goldenfeld 1992). 

At criticality, however, when the effective potential is flat in the center and the 
order-parameter is still small, the phase- amplitude basis has a parametric singularity at 
the origin for po ^ 0. And at higher energies, in the symmetric phase, neither Eq. (16.31) 
nor the phase-amplitude representation is appropriate: then, the particle representation 
in terms of the 0-field is most transparent. 

We are therefore confronted with the situation that different choices of basis in field 
space are useful only in limited regions of the phase diagram. This problem has been 
the subject of an extensive body of literature (Popov 1987, Nepomnashchy 1992 and 
references therein, Pistolesi 1997, Diener 2008, Strack 2008). A derivation of Ward 
identities for the cr-II model can be found in Amit (2005). 

A promising RG approach employing the linear basis in Eq. (16.31) but addition- 
ally invoking constraints from Ward identities has recently been put forward in (Pis- 
tolesi 1997, 2004). How one might smoothly connect the linear basis with the phase- 
amplitude basis via scale-dependent transformations is shown in Diener (2008). 



6.3 Method 53 



6.3 Method 

In this section, we present the functional RG method and our truncation for the 
symmetry-broken phase. We have explained in subsection 13.2.21 in detail how to ac- 
count for spontaneous symmetry-breaking with the exact flow equation for the effec- 
tive action (Wetterich 1993, Berges 2002, Salmhofer 2001, Metzner 2005): 

where r'^^^'^ = d^r^[fP]/d<P^ is the second functional derivative with respect to the 
order parameter field and the trace (Tr) traces over all indices (see section 13.21 for a 
detailed derivation). When expanding this equation around a non- vanishing minimum 
or of the eff'ective action, additional terms will arise in the flow equations for the n- 
point vertex involving the scale-derivative of the order-parameter, a, and the n-l-1- 
point vertex, see Fig. 13.11 

The objective is now to compute the flow of the independent parameters a, u, Z, 
and Y of the cr -H model. A convenient cutoff" function that respects the 0{N)- 
symmetry is a self-energy independent sharp cutoff for space-like momenta: 

R\k) = [Go(k)]-' - [/(k) Go(k)r' (6.9) 

where ;^;''^(k) = 0(|k| - A). This term replaces the bare propagator Go by Gg = X^Gq. 

For a sharp momentum cutoff the momentum variable running around the loop is 
pinned by GJj(k) to |k| = A as the socalled single-scale propagator G'j^ has support 
only for momenta at the cutoff, that is, for |k| = A. Hence the momentum integral can 
be performed analytically. The problem that the integrand contains also step functions 
;t"^(k) = 0(|k| - A) can be treated by using the identity 



/ 



1 

dx 6(x - xo) f[x, 0(x - xo)] = I duf(xo,u), (6.10) 





which is valid for any continuous function /. More specifically, in the present case the 
one-loop diagrams are evaluated for vanishing external momenta and depend only on 
the modulus |k|, such that all internal propagators carry the same momentum and one 
can use the identity 

" / -^G'R(k)A[Gii{k)Ar-' =A'-'KAG(k)A]\\=^, (6.11) 
J (2;r)" 

valid for any matrix A and with being defined by / = J d\k\ \k\''~^. The 
factor n corresponds to the n possible choices of positioning G^ in a loop with n lines. 
For loop integrals in the flow equations, we use the short-hand notation 
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= hr-J T^T^W-A). (6.12) 
Jk\A J (27r) J (InY 

At finite temperatures, the continuous frequency integration becomes a discrete Mat- 
subara sum: J ^ ^ T Xlf3„' with Q„ = InnT . 



6.4 Finite temperature phase boundary in three dimensions 

The primary objective of this section is to extend the finite temperature work presented 
in chapter [5] to systems with Goldstone modes. Most importantly, we compute the 
shift-exponent ij/ of the finite temperature phase boundary in J = 3. 

On the phase boundary at T^, the effective potential is flat and there exists only 
one bosonic excitation. Following Tetradis and Wetterich (1994), we set 7 = in this 
section and therefore: Za- = Z„ = Z. This truncation gives the best results for critical 
exponents and correctly reflects that the phase transition becomes mean-field like with 
vanishing anomalous dimension in the large limit. For an extended discussion we 
refer to section 9 of Tetradis (1994). 



6.4.1 Flow equations 

As demanded by the condition of a vanishing bosonic one-point vertex for all 
scales, see Eq. (13.261) . the square of the order parameter obeys the flow equation 



dAP = - 3GAq) + iN-l) GM , 

\ / 



(6.13) 

The bosonic self-interaction is driven by 

dAU = u^ [ 9Gl(q) + (N-l)Gl{q). 

(6.14) 



V 



+ 




1\A \ , \ 



Lastly, the bosonic field renormalization is 
determined by the expression 



f d 

d^Z = u^p ^ [G,{q + k)GAq)] , 

JcjlA ^ 

(6.15) 



where dl = ^ Y.U ^l,- The flow equations pig. 6.4. Feynman diagrams for the flow 
(I6.13H6.15D are shown in terms of Feynman equations KT3WT5\) . 
diagrams in Fig. [ 
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We split the Matsubara summations into their classical and quantum parts 

(quant) 

T fi^n) = Tf{Qn = 0) + r J] f{Qn) , 



(6.16) 



On 



and utilize the following variables. 



P = 



Zp 



KdTu 



as well as a rescaled temperature and anomalous dimension 

~ InT 
T = 



(6.17) 



77 = 



A 

dXogZ 



J log A 

We therewith obtain for the squared order-parameter and the quartic coupling: 

3 



(6.18) 



J log A 



n=i {[nT)^ + 1 + up^ 
N-l 



+ 



du 



dlogA 



= (d-4 + 2T])u + u 



(1 + upy 



+ (N-1) + 2J2 



+ 



n=i \^{nf) +l + up^ 
N-l 

(6.19) 

Terms proportional to - 1 outside the summation stem from classical Goldstone 
fluctuations. But the Goldstone mode also yields new quantum terms proportional to 
A'^ - 1 inside the summation. The anomalous dimension is determined by 



J] = u p 



1 ^ 1 1 



(6.20) 



The Matsubara summations can be performed analytically yielding hyperbolic trigono- 
metric functions. The rather complex expressions do not deUver any additional in- 
sights. 
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6.4.2 Classical fixed point 

In the classical limit, considering only the zeroth Matsubara frequency/neglecting the 
summation terms, our equations (|6.19[[6^20l) describe the correct classical critical be- 
havior of the Heisenberg universality class for N = 3. Inspection of the first term of 
^ in Eq. (16.191) yields the lower critical dimension d~ = 2 for > 1 in agreement 
with the Mermin-Wagner theorem. The flow equations therefore correctly account for 
the fact that Goldstone fluctuations suppress to zero in two dimensions. The upper 
critical dimension is = 4 from the first term in as expected. The fixed point 
value for the anomalous exponent in three dimensions comes out as: 77 = 0.058 to be 
compared with the known value 77 = 0.036 from other methods (Berges 2002). 

6.4.3 Shift exponent i/f 

The shift exponent determines how the critical temperature varies with control param- 
eter in the immediate vicinity of the QCP: 

Tc<^(6-6,,,,f . (6.21) 

We compute (A employing the same procedure as explained below Eq. (15.371 ) in chap- 
ter [51 The numerical flows look identical to those of the previous section in Fig. 15.71 
except that the scaling plateaus form at the fixed point values of the Heisenberg uni- 
versality class for N = 3. We find that the Ginzburg scale decreases with decreasing 
T, as: 

Ag^T,, (6.22) 

fitting the formula Ac oc j'^li'^-d) derived in the previous chapter, c.f. Eq. (15.381) for 
d = 3. For the shift exponent, we find 

(A « 1/2 , (6.23) 

which matches the formula derived by Millis (1993): iff = j^p^ for J = 3 and z = 1. 
We have verified that for any value of A'^, the shift-exponent attains the Millis value 
1 /2 derived in the symmetric phase without accounting for Goldstone bosons. 

Therefore, the shift-exponent for three dimensional QCPs where a bosonic descrip- 
tion is appropriate does not depend on the number of Goldstone bosons. This is the 
main result of this section. 
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6.5 Infrared asymptotics in the symmetry-broken phase 

Away from the QCP and the phase boundary, the system is characterized by a finite 
minimum in the efi'ective potential. In this section, we clarify the singularity structure 
of the cr - n model in the limit of vanishing cutoff A ^ in the presence of a finite 
order parameter p. We will compute at zero temperature but the obtained results for 
the infrared properties may be transferred to finite temperature by increasing spatial 
dimensionality by one. We establish the scaling properties of the various parameters 
and show that the Goldstone propagator and longitudinal propagator exhibit very dif- 
ferent momentum scaling, in stark contrast to the phase boundary where both modes 
become degenerate. 

Technically, it is crucial to consider the F-term which distinguishes the Goldstone 
field renormalization Z,^ from that of the longitudinal mode Z^-, see Eq. (16.61) . 

6.5.1 Flow equations 

The flow equations for p and u remain unchanged from Eqs. (I6.13l6.14h . The flow 
equation for G^^{k) contains terms from the four diagrams shown in Fig. 16.51 and has 
the form: 



The flow equation for G„^{k) corresponding to the three diagrams in the second 
line of Fig. [63] reads: 




q\A 




q\A 



(6.24) 



'cr 






G~\k\ Z, 



'n 




Fig. 6.5. Feynman diagrams for Eqs. A6.24\6.25\) . 
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d^G-\k) = j^^^ {^-^^^ 3G^(^) + 3^ j G„{q + k) + j^^^ + k) , 

(6.25) 

where the factor of three in front of Ga-{q) in the first term arises from three different 
contractions of the first diagram in the second line of Fig. |63] as drawn in Eq. (4.36) of 
Pistolesi, et al. (2004). One obtains the flow for Z^- and Z,^ by applying the Laplacian 
d\l2 evaluated at |k| = to both sides of Eqs. (16.241 [6^25] ). respectively. The flow of 
Y follows then from the scale derivative of the relation: 

Y = ^^ZZl^ . (6.26) 
P 

6.5.2 Analytical results 

We now identify the leading singularities in the flow equations for w, Zo-, and Z-^ to 
determine their scaling behavior in the infrared, that is, in the limit of vanishing cutoff 
A ^ 0. The order parameter p is finite and does not scale to zero as we are now in the 
symmetry-broken phase away from the phase boundary and the QCP. 

For the quartic coupling m, the diagram with two Goldstone propagators in Fig. 16.41 
is the most singular, leading to the expression: 

d^u = u^ [ (N-\)Gl{q). (6.27) 

Jq\A 

This equation can be integrated analytically and yields in J = 2 at zero temperature: 

u^Z^A, (6.28) 

provided > 1. As the longitudinal mass is proportional to the quartic coupling, 
see Eq. (16.61) . it also vanishes linearly as a function of scale. In three dimensions, the 
linear behavior becomes logarithmic. 

For Zo-, the first diagram in the first line of Fig. [63] contains the leading singularity: 

SaZ^ = "^{N - 1) / ^ [ G^iq + k)GM] ,_o , (6-29) 

^ Jq\A ^ 

Upon inserting Eq. (16.281) . we obtain in two dimensions 




(6.30) 



In three dimensions Zq- diverges only logarithmically. As we will demonstrate below, 
Z„ remains finite and therefore by virtue of Eq. (16.261 ). we obtain for the F-term: 



6.5 Infrared asymptotics in the symmetry-broken phase 
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7 7^ 

y--^^, (6.31) 

where we have used Eq. (16.301) . What remains to be shown is that all singularities for 
the Goldstone propagator cancel such that Z-^ is rendered finite. The leading term for 
the Goldstone propagator from Eq. (16.251) can be written: 

d^G-'(k)= [ 7'J(q)GAq + k), (6.32) 

Jq\A 

where the effective momentum-dependent self -interaction among the Goldstone bosons 
contains two terms: 

7f*(q) = fAq) + 7L2GAq), (6.33) 

with y,fi(q) = \Yq^ and y^„2 = |p ^^^-^j . Since Go-iq) is negative, the interaction 
among Goldstone bosons from exchange of longitudinal fluctuations has the opposite 
sign of the direct interaction y^^ (Pistolesi 2004). By inserting the longitudinal prop- 
agator from Eq. (16.51) and the obtained scale dependences of Eqs. (16.281 16.301 16.311) 
into Eq. (16.331) . we obtain as the central result of this section: 

rjiq) ^ 0, (6.34) 

therefore liberating the corrections to the Goldstone propagator from singularities and 
ensuring the finiteness of Z„ for J > 1 (J > 2) at zero (finite) temperature. In the 
important work by Pistolesi, Castellani, et al. (2004), this proof has been extended to 
arbitrary loop order for the case of an interacting Bose gas. 

We have established the following infrared (A 0) behavior for the longitudinal 
propagator at zero temperature: 

GM) T for ^ = 2 

GM) °^ log A for J = 3 , (6.35) 
while the Goldstone propagator remains quadratic: 

GM) oc ^ for J = 2, 3 , (6.36) 
and the bosonic self-interaction scales as: 

u(A) oc A for J = 2 

"(^)'^ -1^7 forJ = 3. (6.37) 



These results may be transferred to finite temperatures by incrementing spatial dimen- 
sionality d by one. 
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6.6 Conclusion 

This chapter was devoted to an analysis of the effects of Goldstone modes in the 
vicinity of quantum critical points where a description in terms of bosonic fields is 
apppropriate. Within the functional RG framework, we derived flow equations for the 
relativistic cr - 11 model with dynamical exponent z = 1 suitable for the ordered phase 
at zero and finite temperature. 

First, we computed the shift-exponent ijj characterizing the shape of the phase 
boundary at finite temperature in three dimensions and found that ij/ is independent of 
the number of Goldstone bosons in three dimensions. 

We then clarified the infrared properties of the Goldstone and longitudinal prop- 
agator as well as the bosonic self-interaction away from the phase boundary. We 
showed that the longitudinal propagator and the self-interaction are strongly renor- 
malized. The Goldstone propagator, on the other hand, remains quadratic in momenta 
as the effective interaction among Goldstone bosons flows to zero due to symmetries. 

In the future, it would be worthwile to apply an extended version of the flow equa- 
tions devised in this chapter to quantum criticality in the Kosterlitz-Thouless univer- 
sality class (see also Outlook r9.3.2l) and to compute, for example, the shift-exponent 
for the Kosterlitz-Thouless temperature. 
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Fermi-Bose renormalization group for quantum 
critical fermion systems 



7.1 Introduction 

As noticed already more than 30 years ago by Hertz (1976), quantum phase transi- 
tions in correlated fermion systems may be described in terms of the order parameter 
field alone when one integrates out fermions from the path integral in one stroke and 
subsequently deals with an exclusively bosonic theory. The resultant effective action 
is expanded in powers of the bosonic field <p and often truncated after the ^^'^-term. Al- 
though an analysis in terms of ^'^-type theories seems mundane, the presence of two 
relevant energy scales, one given by temperature and the other given by a non-thermal 
control parameter acting as a mass term in the bosonic propagator, gives rise to a rich 
finite temperature phase diagram (Millis 1993). 

The Hertz-Millis approach relies on integrating out fermions first. In general, how- 
ever, the zero temperature fermion propagator is gapless and consequently may lead to 
singular coeflacients in the efi'ective bosonic action (Belitz 2005, Loehneysen 2007). 
In such cases, it is obligatory to keep the gapless fermions in the theory and consider 
them on equal footing to the bosons as performed in several works with (resummed) 
perturbation theory (Altshuler 1995, Abanov 2000, Abanov 2003, Rech 2006). 

Recently, various quantum critical exponents for U(l) gauge theories with Dirac 
fermions and complex-valued bosonic fields in two spatial dimensions were computed 
in the limit of a large number of fermion and boson species Nf and Ni, (Kaul 2008). 
However, the complicated interplay of two singular propagators promotes a controlled 
perturbative treatment to a formidable task (Rech 2006) and represents a clear calling 
for the renormalization group (RG). 

When computing RG flows for coupled fermion-boson theories, where the propa- 
gator of the zero temperature fermions becomes unbounded for momenta on the Fermi 
surface and the propagator for massless bosons is singular at the origin of momentum 
space, it is crucial to synchronize the evolution of correlation functions so that both 
singularities are reached simultaneously at the end of the flow. 
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Within the functional RG framework 
formulated for fermionic and bosonic 
fields (Berges 2002, Baier 2004, Schiitz 
2005, Strack 2008), this poses no partic- 
ular problem as the regulator functions 
for fermions and bosons can be cho- 
sen to synchronize both types of fluctu- 
ations. The mutual feedback of gapless 
fermions coupled to massless bosons 
has been studied already with func- 
tional flow equations in Quantum Elec- 
trodynamics (Gies 2004), non-abelian 
gauge theories (Pawlowski 2004), and 
the Gross-Neveu model (Rosa 2001). 

In the present chapter, we study a 
model of attractive Dirac fermions rel- 
evant for neutral Graphene and cold 
atoms in the half- filled honeycomb lat- 
tice (Zhao 2006, Castro Neto 2007, Cas- 
tro Neto 2008) exhibiting a quantum 

phase transition from a semimetal to a superfluid as shown in Fig. 17. 1[ The flow 
equations are analytically transparent and a simple truncation yields the complete set 
of quantum critical exponents. As a central result, the fermion and order parameter 
two-point correlation functions develop non-analytic dependences on frequency and 
momentum at the QCR 



Fig. 7.1. Schematic phase diagram of at- 
tractive Dirac fermions. The QCP at a 
critical interaction strength Uc separates 
the semimetal from the superfluid. In the 
present chapter, we approach the QCP at 
T = from the semimetallic phase indi- 
cated with an arrow. 



In section 17.21 we introduce the Dirac cone model and show that the mean-field 
theory of this model leads to a semimetal-to-superfluid quantum phase transition at 
a critical interaction strength 11^. In section 1731 the RG method, truncation and flow 
equations are presented. Results for the quantum critical behavior follow in section 
17.41 We finally summarize and conclude in section 1731 



7.2 Dirac cone model 



We consider an attractively interacting Dirac fermion system with the bare action 



roil//, ilf] = - ^kAih - ^k) ^ka- + 



ka- 



k,k'. 



(7.1) 



7.2 Dirac cone model 
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where the variables k = (^o^k) and q = (^o.q) collect the zero temperature Mat- 
subara energies and momenta, and we use the short-hand notation = J^^ J^^ = 

I-oo ^ / momentum and energy integrals, and /^^ includes also a spin sum. 

The dispersion for Dirac fermions with 
the chemical potential directly at the 
Dirac point (see Fig. 17.21) is given by: 



= ±v/|k| 



(7.2) 



with +Vf for the empty cone and -v/ 
for the filled cone. Energy integra- 
tions over are cut off in the ul- 
traviolet and restricted to the band 
{min^k = -^0, max^k = For at- 

tractive interactions the coupling con- 
stant U is negative and drives sponta- 
neous breaking of the global U(l) gauge 
symmetry. Therefore, we decouple the 
Hubbard interaction in the s-wave spin- 
singlet pairing channel by introducing a 
complex bosonic Hubbard-Stratonovich 
field (p^ conjugate to the bilinear com- 
posite of fermionic fields (Popov 1987) 



(7.3) 




Fig. 7.2. Energy spectrum 
dimensional Dirac fermions. 



of two- 



This yields a functional integral over if/, ij/ and (p with the new bare action 
ro[i/f, i/r, (/>] = - / ifJko-{ih - ^k) ^ko- - I fg-rA 



k<T 



' k,q 



(7.4) 



where (p* is the complex conjugate of (p, while ip and if/ are algebraically independent 
Grassmann variables. 



7.2.1 Mean-field theory 

When neglecting bosonic fluctuations by replacing (pq with its expectation value (pq=Q, 
the saddle-point approximation exactly solves the functional integral of Eq. (17.41) and 
leads to the standard BCS gap equation (Popov 1987): 
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^o = -U [ ^' . , (7.5) 

with U and (po independent of momenta. To compute the critical interaction strength 
Uc at which the gap equation has a solution with non-zero 0o, we rearrange Eq. (17.51) : 

-1 _ fdko 
U^.~ J 2^ 

where the density of states has the form: 



A^(a = ^^, (7.7) 

with Kd being defined by / = J d\k\ \k\^~^. The expression in Eq. (|7.6I) is 
equivalent to the Thouless criterion for superconductivity which involves summing 
the particle-particle ladder in the normal phase to infinite order and determining the 
divergence of the effective interaction. Since the density of states vanishes -linearly in 
two dimensions- at the Dirac point, the zero temperature system is stable against for- 
mation of a superfluid for any weak attraction. Instead, a finite attraction larger than a 
certain threshold is necessary to cause superfluidity at zero temperature. Performing 
the integrations in Eq. (17.61) . we obtain the mean-field position of the quantum phase 
transition in Fig. l7.1[ The control parameter provides a mass term for the boson prop- 
agator in Eq. (17.41) and is therefore inversely related to the four-fermion attraction. For 
the physical case of two dimensions, we have: 



1 2nv 



e - ~f^c,MFT - — ; — , (7.8) 

where here Vf has units of energy. For electrons in the honeycomb lattice, the Fermi 
velocity in the vicinity of the Dirac points is proportional to the nearest-neighbor hop- 
ping matrix element t (Castro Neto 2007). Upon setting Vf = Aq = I, the numerical 
value for the mean-field control parameter value is Smft = 0.159. Note that the posi- 
tion of the QCP is non-universal and depends on microscopic parameters such as the 
ultraviolet band cutoff Aq and the Fermi velocity. 

In the following, we will conduct a renormalization group study which enables 
us not only to compute the non-universal renormalized position of the QCP but also 
yields the complete set of universal quantum critical exponents at and in the vicinity 
of the QCP 



7.3 Method 

We derive flow equations for the scale-dependent effective action F/, [i//, if/, 0] within 
the functional RG framework for fermionic and bosonic degrees of freedom in its 
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one-particle irreducible representation (Berges 2002, Baier 2004, Schiitz 2005, Strack 
2008). Starting from the bare fermion-boson action /"a=Ao ['^^'k,^] in Eq. (17.41) . 
fermionic and bosonic fluctuations are integrated along the continuous flow param- 
eter A simultaneously. In the infrared limit A ^ 0, the renormalized, effective action 
T/i^o [lA' 0] is obtained from which physical properties can be extracted. Its RG 
flow is governed by the exact functional flow equation 

-^r'^iS, S] = Str ^^^^ , (7.9) 

where we have collected fermionic and bosonic fields in superfields S, S defined in 
s ub section 13 . 1 . 1 1 and r^^^^ = d^r^[S,S]/dSdS denotes the second functional deriva- 
tive with respect to the superfields and the supertrace (Str) traces over all indices 
with an additional minus sign for fermionic contractions. We refer to Section [3]2] for 
a detailed derivation and discussion of the IPI functional RG framework for both 
fermionic and bosonic fields. 



7.3.1 Truncation 

When evolving towards A ^ 0, infinitely many terms involving fermionic and/or 
bosonic fields with possibly complicated dependences on frequency and momenta are 
generated necessitating a truncation of the effective action. The purpose of this sub- 
section is to explain how the effective action for the model Eq. (17.41) is truncated with 
the objective to capture the most relevant quantum critical renormalization effects. 



Fermion propagator 

To account for a renormalization of the single-particle, fermionic properties by order 
parameter fluctuations, the quadratic fermionic term in the action is modified by a 
field renormalization factor, 

= - i/fka-Zfiiko - ^k) lA/tcr , (V.IO) 

J /to- 
yielding the fermion propagator 

Gf{k) = -{ilfkh) = ^r^ ■ (V.ll) 

A diverging Zf suppresses the quasi-particle weight to zero. If it diverges as a power 
law, Zf ~ IjA^f, the fermion self-energy becomes a non-analytic function of fre- 
quency {ka = io) with 



(7.12) 
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upon identifying the cutoff scale A with oj. As already mentioned above, the Fermi 
velocity is not renormalized separately but kept fixed at unity for the following rea- 
sons. The flow equations for the frequency- and momentum renormalization factors 
are obtained by derivatives of expressions corresponding to IPI-diagrams. The flow 
equation for an independently parametrized Vf would come from one derivative of the 
diagram involving G/ and Gb, shown in the first line of Fig. 17.41 with respect to the 
deviation of momenta from the Dirac point. In parallel, the flow equation for Zy also 
comes from one derivative with respect to frequency of the same expression. Hence, 
both expressions are almost identical for a linearized fermion dispersion and keeping 
Vf does not yield additional information. The initial condition for Zy is Zy = 1. 



Boson propagator 

The bosonic quadratic part of the bare action, Eq. (|7.4I) . consists only of a local mass 
term. Upon integrating out fluctuations, the efl"ective four- fermion interaction medi- 
ated by boson exchange will become momentum- and frequency-dependent. To cap- 
ture this propagating order parameter field, we deploy a renormalization factor Zb 
multiplying the lowest order frequency and momentum terms in a derivative expan- 
sion of the fermionic particle-particle bubble. The real part of the particle-particle 
bubble spanned by fermion propagators endowed with the regulator is a quadratic 
function of external frequencies and momenta leading to the bosonic quadratic part of 
the action: 

r^.^ = [ 0; {Zb {ql + q') +6)4>,. (7.13) 

Jq 

Note that there is no complex linear term in frequency here, as we consider the half- 
filled band and the imaginary linear frequency part of the particle-particle bubble and 
hence the boson propagator vanishes exactly. The control parameter term 6 controls 
the distance to the continuous quantum phase transition and is also renormalized by 
fluctuations. If the initial value of 6 is fine-tuned so that 6^0 for vanishing cutoff" 



boson propagator, parametrized by two RG parameters, reads 



A ^ 0, we are in the quantum critical state with infinite susceptibility, ~ ^U->o- The 



Gb(q) = -(Ml) = ^ , ^ \. , . (7.14) 
^ Zb [ql + q2) + 5 

Recall that this initial ansatz for the boson propagator is motivated by the intermediate 
energy behavior of the particle-particle bubble as the low energy regime is excluded 
from the spectrum by the regulator R'^. The low energy regime is controUably ac- 
cessed subsequently in the coupled Fermi-Bose RG flow when R'^, A 0. We will 
see below that, although we have started with a difl"erentiable function of frequen- 
cies and momenta for the boson propagator at intermediate energies, at criticality in 
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the low energy regime a diverging Z^-factor induces strong analyticities in the bo- 
son propagator -the order parameter field attains an anomalous dimension. The initial 
condition for is Z^ = and for the control parameter 6 = l/U. 



Order parameter self interactions 



We now demonstrate how to account for direct self-interactions of the order parameter 
field within our coupled Fermi-Bose scheme. Recollect that in the conventional Hertz- 
Millis theory of quantum criticality (Hertz 1976, Millis 1993, Loehneysen 2006) one 
expands the purely bosonic action in (even) powers of the bosonic fields: 



u„ 



n>4, even 



2(1)! (I)! 



n/2 



(7.15) 



where the coefficients m„ are generated by fermion loops with n fermion propa- 
gators and n external bosonic legs as shown in Fig. 17.31 At finite temperatures 
and when the fermions are gapped, such expressions are possibly large but finite. 
With the bare, gapless fermion propaga- 
tor at zero temperature defined in Eq. 
(17.1 II) . however, these coefficients are 

singular for small momenta: ^ ^ 



1 



1 



M„ : 



k" k''-^d+l) ' 



(7.16) 




Fig. 7.3. Fermion rings which generate the 
bosonic n-point vertex (here for n = 10). 



leading for instance to a power law sin- 
gularity for the (f)'^ coefficient M4 = w ~ 
l/k in d = 2. The presence of these 
singularities completely invalidates the 
Hertz-Millis approach of expanding the 

eff'ective action in powers of the ordering field alone. Our coupled Fermi-Bose scheme 
is capable of renormalizing such power law singularities, as during the RG flow, 
fermionic fluctuations will generate order parameter self-interactions even if there is 
no ^'^-term present in the bare action. To capture this, we keep the local fourth-order 
term in the eff'ective action: 



riM4 — 



8 



(7.17) 



The initial condition for w is w = 0. Note that the fermion-boson vertex in Eq. (17.41) is 
not renormalized within our truncation. The standard vertex correction one-loop dia- 
gram ~ g^ vanishes by particle conservation. In the following, we keep the fermion- 
boson vertex at its bare value g = I, which can be read off from Eq. (17.41) . 
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7.3.2 Flow equations 

In this subsection, we derive analytic expressions for our flow equations of the trun- 
cated effective action. Both Green functions, Eqs. (17.1 1[ ITTT?!) . display singularities 
for certain choices of momenta, the fermion propagator everywhere in the phase dia- 
gram, and the boson propagator at the QCP when the bosonic mass vanishes. These 
potential infrared singularities for |momenta| < A are regularized by adding optimized 
momentum cutoff's (Litim 2001) for fermions (subscript f) and bosons (subscript b), 

i?/^(k) = Zf {-A sgn [^k] + ^k) e [A - l^kl] 

i?M(q) = (A^ - q2) [A^ - q2] (7.18) 

to the inverse of the propagators in Eqs. (|7.1 1 [177141) . The cutoff-endowed propagators 
are denoted with by GfR and G},r in the following. The scale-derivatives of the cutoffs 
read, 

d^RfA =RfA = -Zf sgn [^k] e[A- l^kl] 

Rt^ = 2Z,Ae[A^-q^] , (7.19) 

where terms proportional to rjf and rjt defined below in Eq. (17.261) are neglected here. 
These additional terms are of higher order in the vertices and subleading. Further 
arguments buttressing this commonly employed procedure are given in (Berges 2002). 

A comment on the choice of relative cutoff scale between fermions and bosons 
is in order. In principle, fermion fluctuations can be cut off at Aj and their bosonic 
counterparts at A/, with both cutoffs being independent functions of A. Without ap- 
proximations, the results do not depend on the concrete choices of Ay(A) and Ab(A). 
In Eq. (|7.18l) . we have chosen Af = Ah = A; the standard choice for critical Fermi- 
Bose theories (Litim 2001, Gies 2004). 

The recipe to obtain the flow equations is now the following: one executes a cutoff- 
derivative acting onRf hA in the analytic expressions corresponding to all 1-loop one- 
particle-irreducible Feynman diagrams for the parameters Zj, Zy, 6, u, which can be 
generated with GfR and GbR as shown in Fig. 17.41 The subsequent trace operation is 
abbreviated by 




(7.20) 



The flow equation for the fermion self-energy is obtained from the Fock-type dia- 
gram in the first line of Fig. 17.41 involving GbR and GfR. This feedback of bosonic 
fluctuations on the fermionic propagator is captured via the flow equation: 



d^Zf = / / dika GfRiq- k) GbR iq) \k=o 



(7.21) 
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This expression vanishes for a local boson (Z^ = 0) from multiple poles in the 
upper complex frequency half-plane. This correctly reflects that for a momentum- 
independent four-fermion interaction, the fermionic Zj-factor is renormalized only at 
the two-loop level. When reinserting the particle-particle bubble which generates Zh 
in the third line of Fig. 17. 41 into the bosonic propagator in the first line of Fig. 17. 4[ we 
indeed observe that effectively our flow equations capture two-loop effects by keeping 
the boson propagator as momentum- and frequency dependent (Z^ 0). 

For the control parameter, we evalu- 
ate the IPI-diagrams with two external 
bosonic legs, and we obtain the two con- jp:!P^ 
tributions: 2/ : , A ^ A, 



d^6 = g^ [ GfAk)GfA-k) 

Jk,Rf 

+ ^ / GtRiq) , (7.22) 

The fermionic contribution on the right- 
hand-side is positive leading to a re- 
duction of 6 for decreasing A whereas 
the bosonic contribution counteracts the 
fermions and tends to increase 6. This is 
the generic behavior of the bosonic fluc- 
tuations as they always tend to restore 
the symmetry. 

The flow of the bosonic frequency 
renormalization is obtained as the sec- 
ond frequency derivative of the particle- 
particle bubble: 




Fig. 7.4. Feynman diagrams representing 
the flow equations. 



d^Z, = / /" l-dlGfR (k + q) GfR i-k) |,=o . (7.23) 

Jk,Rf ^ 

The bosonic tadpole diagram does not contribute here, as the ^'^-vertex u is taken as 
momentum- and frequency-independent. 

Finally, the bosonic self-interaction flows according to: 

Oau = - 4/ / [GfR(-k)nGfR(k)]^ + ^u^ [ [GtRiq)]^ , (7.24) 

Jk,Rf ^ Jq,Rt 

where the first terms generates u and the second, bosonic term tends to reduce u 
in the course of the flow. All frequency and momentum integrations in the above 
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flow equations can be performed analyticaUy. We now elegantly reformulate the flow 
equations by employing the following scaling variables: 

5 ' 



^ A^ZfyfZyfd 
u Ka 

where Kj is defined by J = J d\k\ \k\^''\ By defining the anomalous exponents 
for the fermionic and bosonic Z-factor, respectively, 

d log Zf 
" ~ JlogA 
dlogZh 

rib = --7] T , (V.26) 

J log A 

the explicit dependence on the Z-factors disappears from the expressions. The flow 
equations for the control parameter and the bosonic self-interaction are obtained as: 

dS , r. u 

= (r]b-2)6 + f 



dlogA 4(1+5) 



3/2 



du A 15 

id-3 + 2r],)u-6f + — - -vo ■ (V.27) 



JlogA " ■ "° ■ 16 

The rescaled fermion-boson vertex g obeys the equation: 



J log A 



1 3) 

= [rif + :;Vh-^\~8 . (7.28) 



The equations (17. 27117. 281) have to be considered in conjunction with the fermion and 
boson anomalous exponents: 



3-2 



2 



m = f\ 1 ^ + 77^ I ^ ; ^ • (7-29) 



We will now investigate the analytic properties of these equations and then solve them 
numerically. 



7.4 Solution at the quantum critical point 
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At the quantum critical point, the system becomes scale-invariant, that is, the param- 
eters show no dependence on the scale A. Translated into an algebraic condition, this 
becomes 



d6 



dg 



du 



dlogA dlogA dlogA 



= 



(7.30) 



and one has to solve Eqs. (|7.27[r7?28]) for fixed points together with the anomalous ex- 
ponents of Eq. (I7.29I ). It is immediately apparent that both, the fermion-boson vertex g 
and the boson self-interaction S are relevant couplings below three dimensions. For the 
d = 2 case possibly relevant for graphene and cold atoms in the honeycomb lattice, the 
RG equations for g and u have stable non-Gaussian {g i^O, U 0) solutions with finite 
The anomalous scaling exponents are in- 
terrelated at the QCP: 



3-d 1 



(7.31) 



anomalous exponents rjf and rit, as fol- 
lows directly from Eq. (17.281 ). The values 
of the anomalous exponents from a nu- 
merical solution of Eqs. (I7.27H7.291 ) un- 
der the constraint at the QCP, Eq. (I7.30L 
are shown in Fig. 17.51 

At quantum criticality, the fermionic 
properties of the system cannot be de- 
scribed in terms of conventional Fermi 
liquid theory (Nozieres 1964). A finite 
fermion anomalous dimension entails a 
fermion propagator of the form. 




Fig. 7.5. Fermion and boson anomalous 
exponents at the QCP for \ < d <3. 



Gf(ko,k) oc — 



1 



oc 



with a non-analytic frequency-dependence of the fermion self-energy (^o = 



(7.32) 



a> 



l-r]f 



(7.33) 



where I - rjf = 0.68 for J = 2 as can be read off" from Fig. 17.51 (a). The quasi-particle 
picture breaks down. A similar quantum criticality induced breakdown of the Fermi 
liquid has been noted, among others, in the context of Fermi surface fluctuations near 
a Pomeranchuk instability (Metzner 2003, Dell'Anna 2006 and references therein). 
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There, the frequency exponent of the fermion self-energy at the QCP is ~ or' . 
The fermion propagator also develops a small Luttinger-like anomalous dimension 
in QED3 (Franz 2002) and in the Gross-Neveu model (Rosa 2001, Herbut 2006). 

Concerning the collective properties of the system at quantum criticality, a finite 
boson anomalous dimension entails an order parameter propagator of the form, 

with a non-analytic boson self-energy 

^hiq) - q^-"' , (7.35) 

where 2 - rjh = 1.63 in two dimensions, see Fig. 17.51 (a). In a completely diff"erent 
context, at the antiferromagnetic QCP of the spin-fermion model in two dimensions, 
Chubukov et al. (Abanov 2000, Abanov 2003, Abanov 2004) estimated for the spin 
susceptibility the momentum scaling ~ |q|^^^ within a perturbative 1/A^ calculation 
where is the number of hot spots on the Fermi surface. The physical origin of 
the anomalous momentum exponent for the spin-fermion model is diff"erent from that 
of the present work. There, peculiarities from antiferromagnetic scattering processes 
cause all (/""-vertices to be marginal in J = 2. Summing these logarithms, Chubukov 
et al. then obtain the power law for the quantum critical bosonic self-energy alluded 
to above. This non-analyticity thus vanishes for > 2 as the logarithms are special 
to J = 2. The present work most closely resembles the Gross-Neveu model with 
two fermion flavors where the boson anomalous dimension is however expected to be 
rather large /7z,,Gross-Neveu ~ 0-7 (Rosa 2001, Herbut 2006). 

By comparing the exponents for the quantum critical fermion and order parame- 
ter self-energies of attractive Dirac fermions to other physical contexts, we offer that 
some of the overarching qualitative features generic to various QCPs found with vari- 
ous methods are elegantly accessible within our coupled Fermi-Bose RG framework. 

7.4.1 Quantum critical flows in two dimensions 

We now establish a direct continuous link between the microscopic bare action in Eq. 
(|7.4I) and the infrared properties of the effective action at the QCP. For this purpose, 
we solve the flow equations (|7.27H7.29l) numerically as function of the flow parameter 
A in two dimensions. The initial conditions of our parameters are chosen to precisely 
match the bare action in Eq. (17.41) : a = 0, Z/, = 0, Zj = 1, and g = 1. As in the mean- 
field calculation of subsection 17.2. 1[ the Fermi velocity and the ultraviolet cutoff are 
set to unity Ao = v/ = 1. To reach the quantum critical state, the initial value of 6 is 
fine-tuned so that at the end of the flow, for A ^ 0, 5 ^ 0. 
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Fig. 7.6. (a) Flows of the rescaled control parameter and fermion-boson vertex versus 
s = -log[AMo]. The ultraviolet (infrared) regime is on the left (right) side of the 
plots, (b): Flow of the control parameter 6 at the QCP. 



In Fig. |7.6l (b). we exhibit such a characteristic flow of the control parameter at the 
QCP versus cutoff scale in a double logarithmic plot. We observe scaling behavior 
over more than 20 (!) orders of magnitude limited only by numerical accuracy (note 
the small values of 6 on the vertical axis). The slope is the value of the exponent with 
which the control parameter decreases: 

6 ~ A^""* , (7.36) 

illustrated by the fixed point plateaus of the rescaled control parameter 6 and also 
the rescaled fermion-boson vertex shown in Fig. 17.61 (a). The point we make here is, 
that we can verify whether the effective action is really attracted toward the quantum 
critical fixed point when one starts from the microscopic model and whether this 
behavior is numerically stable. Beyond doubt, this is indeed the case here. 

In Fig. 17. 7[ we present flows for the fermionic and bosonic frequency- and mo- 
mentum renormalization factors Zy and Z/, as well as the associated anomalous di- 
mensions T]f and 7]!,. In Fig. 17.71 (b), we observe that although Z/, is zero initially, it is 
generated for small s (large A) and subsequently diverges as a power law with slope 
rjb'. Zh ~ 1/A''*, which is underlined by fixed point plateaus of rji, shown in Fig. 17.71 
(a). The fermionic Zf is initially equal to unity and then also diverges as a power law 

Zf — , (7.37) 

with a slightly smaller slope than Zb to be read off from the fixed point plateau of ?7/ in 
Fig. 17.71 (a). The numerical solution for both Z- factors exactly fulfills the interrelation 
condition Eq. (|7.31l) and the fixed point values precisely match those of Fig. 17.51 at 
d = 2. We show flows of the ^'^-vertex in Fig. 17.81 Starting from zero initial value, 
fermion fluctuations quickly generate u as is observed in the peak for small s in Fig. 
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Fig. 7.7. (a): Scaling plateaus of the fermion and boson anomalous dimensions at 
quantum criticality versus s = - log [A/Aq]. The fermion anomalous dimension starts 
off at zero ( inset) and becomes finite as soon as Zi, 4^ 0, see discussion below Eq. 
t\7.21\) . (b): Flow of the bosonic and fermionic frequency- and momentum factors Zf 
and Zh- 



17.81 (b). and then the interplay of fermionic and bosonic fluctuations leads to the power 
law scaling behavior 



(7.38) 



again accompanied by fixed point plateaus of u depicted in Fig. 17.81 (a). 



We emphasize that our RG equations not only yield various scaling exponents at 
the QCP, but are also useful in determining non-universal properties of the system. 
The characteristic scale Aqc at which the quantum critical asymptotics sets in, for 
example, can be determined. From Fig. 17.71 we find that the anomalous dimensions 
attain their fixed point values at 5 « 5 and therefore Aqc ~ AqC'^ . 
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Fig. 7.8. (a): Fixed point plateau of the rescaled cp"^ -coupling versus s = - log [A/Aq]. 
(b): Flow of the quartic coupling u. Inset: Small-s behavior when fermionic fluctua- 
tions generate u. 
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Additionally, we can compute the renormalized interaction strength Uc at which 
the system becomes quantum critical and compare it to the mean-field value derived 
in Eq. (17.81) since we have chosen the same values for the upper band cutoff (Aq = 1) 
and the Fermi velocity (v/ = 1) as in the mean-field calculation. Uc is the inverse of the 
fine-tuned initial value of the control parameter {6a=Ao = 0.009032269344423279716) 
such that 6 vanishes at the end of the flow. We find for the ratio between mean-field 
and RG interaction strengths: 

« 17.6 , (7.39) 

therewith renormalizing the mean-field position of the QCP in Fig. 17.11 to the left. 
Fluctuations drastically reduce the size of the superfluid phase. 



7.4.2 Quantum critical exponents 



When approaching a critical point along the relevant parameter axis, the susceptibil- 
ity and correlation length diverge as a power law. The presence of gapless fermions 
places the quantum phase transition considered in this chapter outside the usual 
Ginzburg-Landau- Wilson paradigm, as distinctively signalled by the fermion anoma- 
lous dimension rjf (see Figs. 17.51 and 17.71) which governs the disappearance of quasi- 
particles when approaching the QCP from the semimetallic side, rjf influences all 
other critical exponents similar to the 
universality class of the Gross-Neveu 
model to which the semimetal to anti- 
ferromagnetic insulator QCP for repul- 
sively interacting electrons in the 2d- 
honeycomb lattice has been proposed to 
belong to (Herbut 2006). 

We now compute the susceptibility 
exponent for ^^ = 2 in the vicinity of the 
QCP, 



le-09 - 



le-12 



X 



1 



le-15 



(^-^crit)' 



(7.40) 



le-18 




le-18 



where the inverse susceptibility is iden- 
tified with the (non-rescaled) mass term 
6 of the bosonic propagator, see Eq. 
(|7.14l) . In Fig. 17.91 we present a double- 
logarithmic plot of the susceptibility at 
the end of the flow = ^U^o versus 

the difference of initial values of the control parameter, where (Jcdt corresponds to the 
location of the QCP: A'^i^oC^crit) ^ 0. y can be read off from the slope in Fig. |7.9t 



Fig. 7.9. Logarithmic plot of the inverse 
susceptibility versus 6 in d = 2. 
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r = 1.3 . 



(7.41) 



The correlation length exponent v now follows from the scaling relation (Goldenfeld 



1992, Kaul 2008): y = v(2 - 77^). With 77^ = 0.37 in J = 2 from Fig.llTl one obtains: 



The critical exponents of the structurally similar Gross-Neveu model were computed 
with functional RG methods by Rosa et al. (2001). In a recent paper on a field theory 
containing bosons, fermions, and a gauge field, various quantum critical exponents 
were estimated in the limit of large number of fermion and boson species (Kaul 2008). 

7.5 Conclusion 

In this chapter, we derived coupled flow equations for fermionic and bosonic degrees 
of freedom in quantum critical fermion systems. These equations capture the mutual 
feedback of two distinct types of gapless fluctuations: the first associated with zero 
temperature fermions and the second with massless order parameter fluctuations. As 
a first application, we computed various scaling exponents for the quantum critical 
point between the semimetal and superfluid phase of attractive Dirac fermions. The 
fermion and order parameter propagators are non-analytic functions of frequency and 
momentum. In two dimensions, the fermionic self-energy as function of frequency 
scales ~ o)"^^ leading to a complete breakdown of the Fermi liquid. We demonstrated 
how to compute the susceptibility and correlation length exponents when approaching 
the QCP along the control parameter axis and presented first estimates thereof. 

Extending the work of this chapter to finite temperatures is interesting. For exam- 
ple the correlation length in the quantum critical regime as a function of temperature, 
^{T, ^crit), when approaching the QCP vertically from the top (see Fig. 17. II) . would be 
worthwile to investigate within our coupled Fermi-Bose RG. 

In the next chapter, we combine the bosonic truncations devised for phases with 
broken symmetry in chapters [5] and [6] with the normal phase analysis in terms of 
fermionic and bosonic fields in chapter|7]and we compute various physical properties 
of fermionic superfluids at T = with a coupled fermion-boson action. 



y = 0.8 . 



(7.42) 
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Fermionic superfluids at zero temperature 



8.1 Introduction 

In this chapter, we analyze the attractive Hubbard model as a prototype of a Fermi 
system with a superfluid low temperature phase (Micnas 1990). An experimental re- 
alization of the attractive Hubbard model is conceivable by trapping fermionic atoms 
in an optical lattice and tuning the interaction close to a Feshbach resonance (Hofstet- 
ter 2002, Jaksch 2005, Chin 2006). 

We here theoretically focus on the superfluid ground state for which the impor- 
tance of quantum fluctuations has been emphasized recently in the context of the 
BCS-BEC crossover (Diener 2008). Although the long-range order is not destroyed 
by fluctuations in dimensions d > I, the order parameter correlations are nevertheless 
non-trivial in J < 3. The Goldstone mode leads to severe infrared divergences in per- 
turbation theory. A detailed analysis of the infrared behavior of fermionic superfluids 
has appeared earlier in the mathematical literature (Feldman 1993), where the pertur- 
bative renormalizability of the singularities associated with the Goldstone mode was 
established rigorously. To a large extent divergences of Feynman diagrams cancel due 
to Ward identities, while the remaining singularities require a renormalization group 
treatment. Since the fermions are gapped at low energy scales, the infrared behavior of 
the collective, bosonic sector in fermionic superfluids is equivalent to the one of an in- 
teracting Bose gas, where the Goldstone mode of the condensed state strongly aff"ects 
the longitudinal correlations, leading to drastic deviations from mean-field theory in 
dimensions d < 3 (Nepomnyashchy 1992, Pistolesi 2004) as already demonstrated 
within the functional RG setting in chapter |6l 

Technically, we implement the order parameter via a Hubbard-Stratonovich field 
and we compute the renormalized eff"ective action for the coupled theory of bosons 
and fermions by truncating the exact hierarchy of flow equations for the one-particle 
irreducible vertex functions (Baier 2004, Schiitz 2005, Schiitz 2006). A truncation of 
this hierarchy has been applied a few years ago to the antiferromagnetic state of the 
two-dimensional repulsive Hubbard model (Baier 2004). Important features of the 
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quantum antiferromagnet at low temperatures were captured by the flow. More re- 
cently, various aspects of superfluidity in attractively interacting Fermi systems have 
been studied in the fRG framework with a Hubbard-Stratonovich field for the su- 
perfluid order parameter. Approximate flow equations were discussed previously for 
the superfluid ground state (Krippa 2005), for the Kosterlitz-Thouless transition in 
two-dimensional superfluids (Krahl 2007) and for the BCS-BEC crossover in three- 
dimensional cold atomic Fermi gases (Diehl 2007). 

The purpose of this chapter is to construct a relatively simple truncation of the 
exact fRG flow which is able to describe the correct infrared asymptotic behavior, 
and which yields reasonable estimates for the magnitude of the order parameter at 
least for weak and moderate interaction strength. From a numerical solution of the 
flow equations, which we perform in two dimensions, we obtain information on the 
importance of Goldstone modes and other fluctuation efi'ects. 

In Sec. 18. 2[ we introduce the bare fermion-boson action obtained from the attrac- 
tive Hubbard model by a Hubbard-Stratonovich transformation. Neglecting bosonic 
fluctuations, one recovers the standard mean-field theory for fermionic superfluids, as 
recapitulated in Sec. 18. 3[ By truncating the exact fRG hierarchy, we derive approxi- 
mate flow equations involving fermionic and bosonic fluctuations in Sec. 18.41 At the 
end of that section we reconsider mean-field theory from a flow equation perspective. 
Sec. l8.5l is dedicated to a discussion of results obtained by solving the flow equations. 
We discuss the asymptotic behavior in the infrared limit in two and three dimensions 
and then present numerical results for the flow in two dimensions, where fluctuation 
efi'ects are most pronounced. Finally, we summarize our results in Sec. l8.6[ 

8.2 Bare action 

As a prototype model for the formation of a superfluid ground state in an interacting 
Fermi system we consider the attractive Hubbard model represented by the Hamilto- 
nian 

^ = E E ^'j ^^-^j- + ^ E "'T«a . (8.1) 

ij o- i 

J- 1 

where c-^ and Cj^ are creation and annihilation operators for spin-^ fermions with spin 
orientation cr on a lattice site i. For the hopping matrix we employ t^^ = -f if i and j 
are nearest neighbors on the lattice, and fij = otherwise. On a J-dimensional simple 
cubic lattice, this leads to a dispersion relation = -2? X]f=i '^^^ ^i- the attractive 
Hubbard model the coupling constant U is negative. 

The attractive Hubbard model has a superfluid ground state for any particle den- 
sity n m d > 2 dimensions (Micnas 1990), provided the lattice is not completely 
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filled (n = 2) or empty {n = 0). The Fermi surface for the fiducial case of a quar- 
ter filled band is shown in Fig. I8.1[ At half filling (n = 1) the usual U(l) global 
gauge symmetry becomes a subgroup of a larger S0(3) symmetry group, and the 
order parameter for superfluidity mixes with charge density order (Micnas 1990). 



Our analysis is based on a functional 
integral representation of the eff'ective 
action, that is, the generating functional 
of one-particle irreducible correlation 
functions. For the Hubbard model, the 
starting point is a functional integral over 
fermionic fields ifr and ifr with the bare 
action 



7r/2- 



k, 



-n/2- 




tl^kAiko - ^k) (Aad- 



kcr 



' k,k',q 

(8.2) 



-71/2 7r/2 

kr 



Fig. 8.1. Exemplary Fermi surface for a 
quarter-filled band. 



where = fk - is the single-particle 

energy relative to the chemical potential. The variables k = (ho, k) and q = {qo, q) col- 
lect Matsubara energies and momenta. We use the short-hand notation = /^^ = 

I-oo ^ I-7T (i^ momentum and energy integrals, and /^^ includes also a spin sum. 
We consider only ground state properties, so that the energy variables are continuous. 

The attractive interaction drives spin-singlet pairing with s-wave symmetry and a 
spontaneous breaking of the global U(l) gauge symmetry. Therefore, we decouple 
the Hubbard interaction in the s-wave spin-singlet pairing channel by introducing a 
complex bosonic Hubbard-Stratonovich field 0^ conjugate to the bilinear composite 
of fermionic fields (Popov 1987) 



(8.3) 



^q = U jj/k^1;ll/-k^l ■ 

This yields a functional integral over ij/, ij/ and with the new bare action 

ro[l/^, (A, 0] = - / tffkAik) - ^k) ^ko- - I 4>*qJj(Pq 



kcr 



k,q 



(8.4) 



where 0* is the complex conjugate of (p, while ifj and ij/ are algebraically indepen- 
dent Grassmann variables. Our aim is to compute fermionic and bosonic correlation 
functions with focus on the correct description of the low energy (infrared) behavior. 
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8.3 Mean-field theory 

As a warm-up for the renormalization group treatment it is instructive to recapitulate 
the mean-field theory for the superfluid phase in the functional integral formalism 
(Popov 1987). In mean-field approximation bosonic fluctuations are neglected, that is, 
the bosonic field is fixed instead of being integrated over all possible configurations. 
The fermion fields can be integrated exactly. The (fixed) bosonic field is determined 
by minimizing the eff'ective action as a functional of (p. For a homogeneous system, 
the minimizing (pq can be be non-zero only for ^ = 0. We denote the minimum by a. 
Substituting cpo ^ a + yields 

ro[iA, lA, a + 0] = - / iffka-(iko - ^k) ipka- - a* + [ {iJf-ki^k-[ oc -I- ^k\^-ki ot) 

Jkcr ^ Jk 

-— {a*(t>o + a(pl) + / (iA-,t+fi'Ayt+fT (f>q + ^Ai+f T'A-/i+fi (f>\ 

'-' J k,q 

A necessary condition for a minimum of the effective action is that its first derivative 
with respect to cf) (or 0*), that is, the bosonic 1-point function r'"^\q), vanishes. In other 
words, terms linear in (p (or (p*) have to vanish in the effective action. For q 0, r''^\q) 
vanishes for any choice of or in a homogeneous system. For q = and in mean-field 
approximation, the 1-point function is given by 

r^^\0) = ~a + jiiPu^iP-ui) , (8.6) 

where (. . .) denotes expectation values. The first term on the right hand side corre- 
sponds to the contribution -jjacpQ to Fq in the second line of Eq. (18.51) . while the 
second term is generated by contracting the fermions in the contribution proportional 
to 0* in the second line of Eq. (18.51) . In the absence of bosonic fluctuations there is no 
other contribution to F^p. From the condition F'jj\0) = one obtains 

a = U J{iPkrip-ki) , (8.7) 

which relates or to a fermionic expectation value. We now turn to the fermionic 2-point 
functions. The normal fermionic propagator G/o-ik) = -{if/ko-^Pkcr) and the anomalous 
propagators Ff{k) = -(ip^ip^ki) and Ff{k) = -{ip^kiipkf) can be conveniently col- 
lected in a Nambu matrix propagator 
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The anomalous propagators satisfy the relations Ff(k) = F*j(k) and Ff{-k) = Ff{k). 
In (our) case of spin rotation invariance the normal propagator does not depend on cr, 
and therefore, Gf^{k) = Gfi(k) = Gfik). 

In mean-field theory, the fermionic 2-point vertex function Vj^ = -G'^/ can be 
read off directly from the bare action in the form Eq. (|8.5I) : 



The off-diagonal elements are due to the last bracket in the first line of Eq. (18.51) . 
Tadpole contributions which are generated from the terms in the second and third 
line of Eq. (18.51) cancel exactly by virtue of the condition F^^^ = 0. In the absence of 
bosonic fluctuations there are no other contributions to r^^\ Inverting F^' and using 
^-k = from reflection symmetry yields 



/Cq -I- 

^#) = 7I^, (8.11) 

where = (^^ -I- \Af-y^^ and A = a. We observe that in mean-field theory the bosonic 
order parameter a is equivalent to the gap A in the fermionic excitation spectrum. 
Eq. (|8.6I) corresponds to the BCS gap equation 

A = -U FAk) . (8.12) 



We finally compute the bosonic 2-point functions in mean-field theory (see Keller 
(1999) for a more detailed explanation of the diagrammatic ingredients and the 
notation in the superfluid phase). The bosonic propagators Gb{q) = -{(pq(fQ and 
Fbiq) = -{(pq(p-q) = -{f-q(p*q)* form thc matrix propagator 

G.W=f!:*L^*0. (8.13) 



FHq) Gt{-q) 

Note that Fhi-q) = Fh(q). The bosonic 2-point function Vf^ is equal to -G^^ We 
define a bosonic self-energy E/, via the Dyson equation (GhY^ = (G^)"' - ^fc, where 
the bare propagator corresponding to the bare action Fq is given by 

G,o(^)= (qI)- (8.14) 



In mean-field theory, only fermionic bubble diagrams contribute to the bosonic self- 
energy: 
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Fig. 8.2. Momentum (left) and frequency (right) dependence of K(q) in the normal 
phase. The small linear frequency dependence in the imaginary part vanishes at half- 
filling and is discussed in Sec. \8.4\ 



f K(q) L(q) 
\L*iq)K{-q) 



where 



K(q) = - Gf{k + q)Gf{-k) 



L(q)= J^Ff(k + q)Ff(-k), 



(8.15) 

(8.16) 
(8.17) 



see Fig. 18.21 for an exemplary plot of K(q) in the normal phase. In the absence of 
bosonic fluctuations, there are no other contributions to E/,. Tadpole diagrams cancel 
due to rl^^ = 0. Note that K{-q) = K*{q) while L{-q) = L{q). Inverting the matrix 
Gio - one obtains the bosonic propagator in mean-field approximation 



Gb{q) = 



( - K{-q) L(q) 
d(q)\ L\q) U-'-Kiq) 



(8.18) 



with the determinant d(q) = \U ^ - K(q)\^ - \L{q)\^. Using the explicit expressions 
(18.101) and (18.111) for G/ and Ff, respectively, one can see that 



- K(0) + |L(0)| = U-'+^ jPfik) , 



(8.19) 
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which vanishes if A is non-zero and satisfies the gap equation. Hence d(q) has a zero 
and Ghiq) a pole in ^ = 0. This pole corresponds to the Goldstone mode associated 
with the sponaneous breaking of the U(l) symmetry of the model. For small finite q 
and qo, the leading ^-dependences of \U~^ - K{q)\ and L{q) are of order |qp and q^. 
Hence the divergence of Gb(q) for ^ ^ is quadratic in q and qo. Continuing qo to 
real frequencies one obtains a propagating mode with a linear dispersion relation. The 
second pole of d(q) is gapped and features a quadratic momentum dispersion. 

By appropriately tailoring the fRG-trunction in Sec. l8.4l to the results of this mean- 
field calculation in the superfluid phase, we incorporate the effects of transversal Gold- 
stone fluctuations as well as longitudinal fluctuations into our computation. 



8.4 Truncation 

At the core of our analysis is the exact functional flow equation for the scale- 
dependent effective action F^, which generates all IPI correlation functions (Wet- 
terich 1993, Berges 2002, Salmhofer 2001, Metzner 2005): 

—r'^lS, S] = Str 5t , (8.20) 

where r^^^^ = d^r^[S,S]/dSdS denotes the second functional derivative with re- 
spect to the superfields and the supertrace (str) traces over all indices with an addi- 
tional minus sign for fermionic contractions. We refer to Section [X2l for a detailed 
exposition of the functional RG formalism including symmetry -breaking. 

Untruncated, the exact effective action contains an infinite number of terms of ar- 
bitrary order in fermionic and/or bosonic fields. We now describe which terms are 
kept and how they are parametrized. We keep all terms which are crucial for a quali- 
tatively correct description of the low-energy behavior of the system. We distinguish 
the symmetric regime, where a = 0, from the symmetry broken regime, where a 0. 
The former applies to large A, the latter to small A. 



8.4.1 Symmetric regime 

The bare action Eq. (18.41 ) contains quadratic terms for fermions and bosons, and an 
interaction term where bosons couple linearly to a fermion bilinear. In the effective 
action we keep these terms with generalized cutoff-dependent parameters and add 
a bosonic self-interaction of order The latter is generated by the flow and be- 
comes crucial when the quadratic part of the bosonic potential changes sign. Other 
interactions generated by the flow are neglected. For our choice of parameters (rela- 
tively small U), the fermionic propagator receives only Fermi liquid renormalizations, 
leading to a slightly reduced quasi particle weight and a weakly renormalized disper- 
sion relation. We neglect these quantitatively small effects and leave the quadratic 
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fermionic term in the action unrenormal- 
ized, that is, 



Gf{k) 



^ko-iiko - ^ko- , (8.21) 



/to- 



corresponding to an 
fermionic propagator 

Gf{k) = Gfoik) 



unrenormalized 



(8.22) 



Gbiq) 



A = 




iko - 



Fig. 8.3. Diagrammatic constituents of our 
truncation in the symmetric regime as de- 
scribed in Sec. \8.4.1\ 



In the bare action the term quadratic in 
bosons contains only a mass term. In the 

effective action this mass decreases with decreasing cutoff until it vanishes at a criti- 
cal scale Ac, which marks the transition to the symmetry-broken regime. As the mass 
decreases, the momentum and frequency dependence of the bosonic 2-point func- 
tion becomes important. The latter is generated in particular by fermionic fluctua- 
tions. For small q, the leading q-dependence is of order |qp. The leading frequency- 
dependent contribution to the real part of the bosonic 2-point function is of order 
qly The frequency-dependence of the imaginary part is generally of order but the 
prefactor is very small which is related to the fact that it vanishes completely in case 
of particle-hole symmetry. Furthermore this small imaginary part does not have any 
qualitative impact on the quantities we compute in the following. We therefore neglect 
this term and make the ansatz 



f^{ml + Zbql + Abiol) (pg , 



(8.23) 



where oj^ - 



= 2J2i=i (1 ~ cos qi) is fixed, while ml, Zb, and Ab are cutoff-dependent 
numbers. The function co^ has been chosen such that the quadratic momentum depen- 
dence for small q is continued to a periodic function defined on the entire Brillouin 
zone. The initial conditions for the parameters in the bosonic 2-point function. 



Gbiq) = - 



Zbql + ^bcol + ml 



(8.24) 



can be read off from the bare action as ml = \2/U\ and Zb = Ab = 0. 

The interaction between fermions and bosons remains regular and finite near Ac. It 
can therefore be parametrized as 



fq] , (8.25) 

J k,q 



where the coupling constant g depends on the cutoff, but not on momentum and fre- 
quency. The initial condition for ^ is ^ = 1. 
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The flow generates a bosonic self-interaction which plays a crucial role near and 
in the symmetry-broken regime, that is, when the bosonic mass term becomes small 
and finally changes sign. The most relevant term is a local 



il'^-interaction 



rui4 = - I ^^^piij^,_piij,fipci 

O ■lq,q',p 



(8.26) 



where the coupling constant X depends on the cutoff but not on momentum and fre- 
quency. The initial condition for Ah A = 0. 

The propagators and vertices are represented diagrammatically in Fig. l8.3[ 

8.4.2 Symmetry-broken regime 

For A < Ac the efl'ective action develops a minimum at (pq=Q = a 0. Due to the 
U(l) symmetry associated with charge conservation the minimum is degenerate with 
respect to the phase of a. We here employ the cr-II model presented previously in 
section \62\ consisting of the quartic coupling A and the minimum a, see Eqs. (16. 2[ 
I6.3[|6.7I) . In this chapter we focus on the case N = 2 degrading the vector field 11 to a 
one-component field n representing the Goldstone mode. The leading momentum and 
frequency dependence of the bosonic 2-point function is quadratic in q and qo, both 
for the cr- and ;7r-component. Hence 
we make the following ansatz for the 
quadratic bosonic contributions to the ef- 
fective action 



Gf(k) = 



Ffik) = 




F<rcr = ^ J cr-cj{ml + Z^ql + A^oj\) cr^ G^k) - - - 

If 2 2 

(8.27) 

where mo-, Zo-, Ao-, Z„, and A„ are cut- 
off" dependent real numbers. The propa- 8.4. Propagators and fermion-boson 

gators for the cr and n fields thus have vertex for our truncation in the symmetry- 
broken regime as specified in Sec. \8.4.2\ 

1 The bosonic self-interactions are exhibited 

ml + Z^ql+A^co\ P^S- ^ 

'^-(^^ = - 7 il, 2 ■ (8.28) 
The longitudinal mass is determined by the 101^^ coupling A and the minimum a as 



the form 

GAl) = 



ml = A lap 



(8.29) 
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The small imaginary contribution to T^.^ generated by fermionic fluctuations in the 
absence of particle-hole symmetry, mentioned already above, gives rise to an off- 
diagonal quadratic term Fa-n with a contribution linear in q^. Since this term would 
complicate the analysis without having any significant effect, we will neglect it. 

We now discuss terms involving fermions in the case of symmetry breaking. In 
addition to the normal quadratic fermionic term T^^^ defined as before, see Eq. (I8.21|) . 
the anomalous term 

F^^ = jy {p.uih^ + zIV.T'A-^i) (8.30) 

is generated in the symmetry -broken regime, where is a cutoff-dependent energy 
gap. The phase of A is inherited from the phase of a while its modulus is gener- 
ally different, due to fluctuations. Since we have chosen a real and positive, A is real 
and positive, too. The normal and anomalous fermionic propagators Gf and F f cor- 
responding to and F^tjj have the standard mean-field form as in Eqs. (|8.10l) and 
(|8.11l) . with £k = (^k + Ml^)'^^' but now A is not equal to a. 

In addition to the interaction between fermions and bosons of the form Eq. (I8.25|) . 
an anomalous term of the form 

J k,q 

is generated in the symmetry-broken regime. Inserting the decomposition of (p in lon- 
gitudinal and transverse fields into the normal and anomalous interaction terms, we 
obtain 

/ ('A-zt+uiA/t+fTcr? + ^A/t+fTiA-zt+fiCT-g) . (8.32) 

J k,q 

r^2„ = ig„ Vp-k^i^k+l^^q-^k^^^-k+li^-q) ^ (8.33) 

J k.q 

where go- = g + 8 and g„ = g - g. Fermions couple with different strength to the cr- 
and TT- field, respectively. 

A diagrammatic represention of the various progagators and interaction vertices in 
the symmetry-broken regime is shown in Fig. 18.41 

We finally note that in a previously reported truncation (Birse 2005) of the fRG 
flow in a fermionic superfluid, no distinction between longitudinal and transverse 
fields was made for the bosonic Z-factors in the symmetry-broken regime. 



8.4.3 Flow equations 

Inserting the above ansatz for the truncated effective action into the exact flow equa- 
tion and comparing coefficients yields a set of coupled flow equations for the cutoff 
dependent parameters. The various contributions can be conveniently represented by 
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Feynman diagrams. The prefactors and signs in the flow equations could be extracted 
from the expansion of the exact flow equation, Eq. (|3.20l) . However, in practice we 
determine them by comparison to a conventional perturbation expansion. 

All contributions to our flow equations correspond to one-loop diagrams with only 
one momentum and frequency integration, as dictated by the structure of the exact 
flow equation in the form (13.201) . One of the propagators in the loop is a bosonic 
or fermionic component of the single-scale propagator G^, the others (if any) are 
components of Gr. 

In this chapter, we use sharp frequency cutoffs which exclude bosonic fields with 
Ifrequencyl < and fermionic fields with |frequency| < Af from the functional 
integral. Thereby both fermionic and bosonic infrared divergences are regularized. 

Both cutoffs are monotonic functions of the flow parameter, Afo(A) and A/(A), 
which vanish for A ^ and tend to infinity for A ^ oo. Concretely, we employ 

Rl*(fc) = [G,o(fc)]-^ - [;rl*(fco) G,o(fc)]-^ (8.34) 
for s = b,f, and;^^(A:o) = ©{\h)\ - As). This term replaces the bare propagators G.vo by 

For a sharp frequency cutoff the frequency variable running around the loop is 
pinned by G^(fco) to = ±^fo or ko = +Af as the socalled single-scale propaga- 
tor G^'^ has support only for frequencies at the cutoffs, that is, for |fcol = As. Hence 
the frequency integral can be performed analytically. The problem that the integrand 
contains also step functions ;^ff(^o) = 6>(l^ol - As) can be treated by using the identity 

j dx6{x - xo) f{x,0{x - xo)\ = j duf(xo,u), (8.35) 

which is valid for any continuous function /. 

More specifically, in the present case the one-loop diagrams are evaluated for 
vanishing external frequencies, such that all internal propagators carry the same fre- 
quency. In loops involving only either only bosonic or only fermionic propagators, 
one can use the identity 

n I dh G',«(fco) A [GsRiko) A]"-' = a; [Gs(ko) Af , (8.36) 

valid for any matrix A, to replace the frequency integration by a frequency sum over 
+As while replacing all the propagators in the loop by Gv. The factor n corresponds to 
the n possible choices of positioning G',^ in a loop with n lines, and A'^, = dAs/dA. For 
Ah = Af the above formula holds also for the superpropagator G^, such that is applies 
also to loops with mixed products of bosonic and fermionic propagators. For At Af, 
mixed loops contribute only if the single-scale propagator is associated with the larger 
cutoff. For example, for At > Af, the single-scale propagator has to be bosonic, since 
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Gi vanishes at |fcol = On the other hand, for |fcol = one has G/rC^q) = Gfiko), 
and for the integration of the bosonic factors in the loop one can again use Eq. (18.361) . 

For loop integrals involving the frequency sum over +As and the momentum inte- 
gral over the Brillouin zone we use the short-hand notation 




(8.37) 



where A'^^ = dAJdA. 

The frequency cutoff is convenient for fermion-boson theories at zero temperature 
as the singularity for both particle species is situated at the origin \qQ\ = on the fre- 
quency axis. This makes mixed diagrams particularly easy to treat as alluded to above. 
The downside, on the other hand, is that cutting off frequencies ruins the analyticity 
properties of Green functions in the complex frequency plane. In the present chapter, 
this issue is not important but in other parts of this thesis, when for example calculat- 
ing a Z- factor for the fermion self-energy as performed in chapter |71 we resorted to 
other cutoffs. 

Symmetric regime 

Here we choose Ah = Af = A. One is in principle free to choose the fermionic and 
bosonic cutoff independently. We have checked that the concrete choice of Ah and 
Af in the symmetric regime does not change the final results for A ^ much. The 
diagrams contributing to the flow in the symmetric regime are shown in Fig. 18.51 

The flow of the bosonic mass is given by the bosonic self-energy at vanishing 
external momentum and frequency, that is, 

5a^=// Gf{k)Gf{-k) + ^ f Ghiq). (8.38) 

Jk\A ^ Jq\A 

The fermionic contribution to d/^rr^ is positive, leading to a reduction of upon 
decreasing A, while the bosonic fluctuation term is negative (since Gh{q) < 0). The 
flow of Zh is obtained from the second frequency derivative of the bosonic self-energy 
as 

d^Zh =g^ [ dlPfik + q) Gfi-k)\ . (8.39) 

Jk\A ^ 

Similarly, the flow of Ah is obtained from a second momentum derivative of the 
bosonic self-energy: 

d^Ah = / dlGfik + q) Gf{-k)\^^^^ , (8.40) 
where 5q = 2Yli=[d]ir Since the bosonic self-energy is isotropic in q to leading 
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(quadratic) order in q, the results do not 
depend on the direction in which the mo- 
mentum derivative is taken. The bosonic 
tadpole diagram in Fig. 18.51 contributes 
only to nil,, not to and A/,, since it 
yields a momentum and frequency inde- 
pendent contribution to the self-energy. 

Finally, the flow of the l^l"^ coupling 
is given by 

Sa^ = -^g' [ [Gf(k)f[Gf{-k)f 

Jk\A 

+ ^A^ I VGi,{q)f . (8.41) 




Within the truncation of the effective ac- 
tion described in Sec. [8AT] there is no Fig. 8.5. Feynman diagrams representing 
contribution to the flow of the interac- the flow equations in the symmetric regime. 
tion between fermions and bosons in 
the symmetric phase. The coupling g re- 
mains therefore invariant. 



Symmetry broken regime 

In the limit A ^ we are forced to choose Af «; to avoid an artificial strong 
coupling problem, as will become clear below. We therefore choose Af < A/, in the 
entire symmetry broken regime, which implies that the frequencies in mixed loops 
with bosonic and fermionic propagators are pinned at the bosonic cutoff". The precise 
choice of the cutoffs will be specified later. 

We first derive the ffow equation for 
the minimum of the bosonic potential a, 
which is derived from the condition that 
the bosonic 1 -point vertex 7"*^'^ be zero 
for all A. The ffow equation for F^^^ reads 



9 V Q 



dAFy = midAa + 2g^ Ff(k) 

/ [3G^{q) + G„{q)] ■ 



2 



Fig. 8.6. Contributions to the bosonic 
(8.42) p^Ij^j^ vertex and fermion gap below Ac. 



The various contributions are repre- 
sented diagrammatically in Fig. 18.61 
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The first term is due to the cutoff-dependence of the expansion point around which 
the effective action is expanded in powers of the fields. The condition dy^F^^^ = 
yields 



Ff(k) 



k\Af 



1 

2^ 



[3GAq) + G,iq)] . 



(8.43) 



We have used Eq. (18.291 ) to simplify the last term. The fermionic contribution to d^a is 
negative, leading to an increase of a upon decreasing A, while the bosonic fluctuation 
term is positive and therefore reduces a. The behavior of Eq. (18.431) in the vicinity of 
the critical scale, A < Ac, when a and are small, is shown below in Sec. 18.5.11 

The flow of A is obtained from the flow of the anomalous component of the 
fermionic self-energy as 



= gcr - / Ff{q - k) [glG^(q) - glG„(q)] |i=(o,k^) 

JqlAb 



(8.44) 



The first term, due to the cutoff-dependence of the expansion point for the effective 
action, links the flow of the fermionic gap to the flow of the bosonic order parameter. 
The second term captures a correction to the relation between a and A due to bosonic 
fluctuations as illustrated in Fig. 18.61 

The flow of the mass of the longi- 
tudinal order parameter fluctuations (cf. 



Fig. 18.71) is obtained from the self-energy 
of the cr fields at zero momentum and 
frequency. The flow of A can then be 
computed from the flow of and a via 
the relation Eq. (18.291 ). For m^, we have 
the flow equation 



--0-- 



■O 



- -I- 



dA^- 



-I- 




— -I- 



Gf(k)Gf(-k) 



k\Af 

- Ff(k)Ffi-k) 
+ ^ [ [3GAq) + GM] 



+ 



N ^ ^ ' 



■o- 



+ 



' q\Ab 
{Aaf 

Aa 



q\Ab 



{9Gl{q) + Gliq)] 



+ 3— OaO: 



(8.45) 



Fig. 8.7. Diagrammatic representation of 
the contributions to the bosonic mass in the 
symmetry -broken regime. 



The second term in this equation is due to a product of the 3-point vertex yg-a and 5/iQr 
arising from the cutoff dependence of the expansion point for the effective action. The 
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7 A 



Fig. 8.8. Diagrams renormalizing the bosonic Z- and A-factors for A < Ac. 



flow of Zo- is obtained from the second frequency derivative of the self-energy of the 
cr fields, which yields 



d^Z^ = gl[ dllGfik + q) Gf(-k) - Ff(k + q) 

Jk\At 



lk\Af 

+^ / {9G^{k + q) G^k) + GAk + q) G,{k)] 

^ Jk\A,, 



q=0 



q=0 



(8.46) 



The flow of Aa- is given by the same equation with d^^ replaced by d^. 

In the flow of Z„ there are strong cancellations of different terms originating from 
bosonic fluctuations. We have shown with a slightly extended truncation in chapter 
[6] that the efi'ective self-interaction among Goldstone bosons flows to zero such that 
Z;r remains finite and the Goldstone mode is not renormalized substantially (Pistolesi 
2004). We therefore keep only the fermionic fluctuations, that is, 

d^Z, = gl[ dl[Gfik + q)Gfi-k) + Ffik + q)Ff(-k)]\ . (8.47) 

Jk\Af ^ 

For the flow of A„ we obtain the same equation with d^^^ replaced by d^. The terms 
contributing to the flow of the Z- and A-factors are illustrated in Fig. 18. 8[ 

In the symmetry broken regime, there are also contributions to the flow of the 
interaction between fermions and bosons due to vertex corrections with bosonic fluc- 
tuations (see Fig. 18.91) , yielding 



dAgo- = go 



dAgn = gn 



q\Ab 



F]{k-q)-\Gfik-q) 



Fl{k-q)+\Gf{k-qf 



J *:=(0,kf ) 



<:=(0,kf) 



[glGAq)-glGAq)\ 

[glGM)-glGM] . 



(8.48) 



The right hand sides are dominated by the contribution from the n propagator, which 
tends to reduce go-, gn for decreasing A. 
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gn ■■ 




Fig. 8.9. Fermion-boson vertex corrections below Ac. 



8.4.4 Relation to mean-field theory 

Before solving the flow equations derived above, we first analyze what happens when 
contributions due to bosonic fluctuations are neglected, and relate the reduced set of 
equations to the usual mean-field theory (Sec. 18.31) . 

In the absence of bosonic fluctuations, g = go- = gn = 1- Furthermore, the bosonic 
order parameter a and the fermionic gap A are identical: a = A. Since the bosonic 
cutofl" is irrelevant here, we can choose Ay = A. The flow equation for the bosonic 
mass in the symmetric regime, Eq. (18.381) . simplifies to 



^ Jk\A 

where Gf{k) = Gfo(k) = (iko - ^k)~^- This equation can be easily integrated, yielding 



nih vanishes at a critical scale Ac > 0. The flow equation for A (= a) in the symmetry 
broken regime A < A^. , Eq. (18.441 ). is reduced to 



It is complemented by the flow equation for the mass of the cr field, Eq. (|8.45l) . which 
becomes 




(8.49) 




(8.50) 




(8.51) 




(8.52) 
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in the absence of bosonic fluctuations, with y^-s = AA/2 = m^/(2zl). The propagators 
Gf and Ff have the usual BCS form, as in Eqs. (18.101) and (18.111) . 

A numerical solution of the coupled flow equations (18.511) and (18.521) yields a gap 
A which is a bit smaller than the BCS result obtained from the gap equation (18.121) . 
The reason for this discrepancy is the relatively simple quartic ansatz (16.21) for the 
bosonic potential. The complete bosonic potential is non-polynomial in |0p even in 
mean-field theory. Restricted to the zero momentum and frequency component of it 
has the form (Popov 1987) 

t,MP(^)=l^-/-,„M^filW!. (8.53) 



1^1 Jt'" *ii + fk 

The kernel of the 3 -point vertex y^i obtained from an expansion of this mean-field 
potential around a finite order parameter A reads 



7^3 = -2 



^-F}ik)-Ff(k)\Gf(kf 



(8.54) 



at zero frequencies and momenta. Inserting this into (18.521) . the flow of can be 
written as a total derivative 

Oa"^ = -dA [ [ [\Gf{kf - Fj{k)] , (8.55) 

^ J\ka\>A Jk 

where the A-derivative on the right hand side acts also on A, generating the term 
proportional to y^3. Integrating this equation with the initial condition nia- = at 
A = Ac, yields 

2 

= K{0) + L(0) - , (8.56) 

which is the correct mean-field result. With nia- given by (|8.56h . the flow equation 
(18.511) yields the correct mean-field gap. The easiest way to see this, is to write the 
BCS gap equation in the presence of a cutoff" in the form 1 = -U ij^p|>^ /k^ ' Ffik), 
and take a derivative with respect to A. 

It is instructive to relate the above flow equations for A and mo- to the flow equations 
for the BCS mean-field model obtained in a purely fermionic RG (Salmhofer 2004). 
For a sharp frequency cutoff, those flow equations have the form 



d^A = -(V + W) Ffik) , (8.57) 
d^V+W) = {V+Wf dJ [ [\Gf(k)\^ - F}(k)] , (8.58) 

J\kn\>A Jk 



where V is a normal two-fermion interaction in the Cooper channel, while W is an 
anomalous interaction corresponding to annihilation (or creation) of four particles. 
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With the identification = V + V7 these equations are obviously equivalent to 

(18.511) and (18.551) . The above flow equation iox V + W is obtained from a one-loop 
truncation complemented by additional self-energy insertions drawn from higher or- 
der diagrams with tadpoles (Salmhofer 2004, Katanin 2004). These additional terms 
correspond to the contractions with yo-3 in the present bosonized RG. 



8.5 Results 

In subsections l8.5.1l and l8.5.2[ analytic properties of our flow equations are discussed. 
Numerical results for two-dimensional systems are exhibited and interpreted in sub- 
section [8331 



8.5.1 Flow for A < A, 

For A slightly below Ac, the flow equations can be expanded in the order parameter. 
To leading order, the order parameter a and the gap A are identical, a = A . The 
fluctuation term in the flow equation (18.441) for A is quadratic in a, and also the flow 
of go- yields only corrections beyond linear order to the relation between a andzl. 
Near A^, the flow equation (18.431) for a can be written as 



d^a^ = -1 / - / \3GM) + GM^ , (8.59) 

lq\At 



A 



where / = /^|^^(^o + Ik) ^ evaluated for A = Ac- Note that we have replaced the ratio 
ml./a^ by A in the first term on the right hand side of the flow equation. Integrating 
the flow equation one obtains 



q\Ab 



A 



{Ac -A) (8.60) 

A=Ac 



for A < Ac. The order parameter a and the fermionic gap are thus proportional to 
{Ac- Afl^ for A < Ac. 

Inserting Eq. (|8.60l) into the flow equation (18.451) for and neglecting the last 
fluctuation term, which is of higher order in a, one obtains 



d^mi = -4I-A [3GAq) + G„{q)] ■ (8-61) 

Jq\Ab 

This shows that the flow of a and is indeed consistent with the relation = Aa^ 
following from the ansatz for the bosonic potential. 
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8.5.2 Infrared asymptotics 

In the infrared limit (A 0), the key properties of the flow can be extracted from 
the flow equations analytically as already demonstrated in chapter [6] for the general 
cr-n model with A^- 1 Goldstone bosons. The behavior of the bosonic sector depends 
strongly on the dimensionality of the system. We consider dimensions d > 2, focus- 
ing in particular on the two- and three-dimensional case. The bosonic order parameter 
and the fermionic gap saturate at finite values in the limit A ^ 0. The fluctuation 
corrections to d^a and involve the singular Goldstone propagator Gj^ only lin- 
early and are therefore integrable in d > 1. The fermion-boson interactions go- and 
g„ also saturate. The finiteness of Zj^ and A„ is guaranteed by the correct implemen- 
tation of symmetries (see chapter [6l Pistolesi 2004). We choose At, = A in the fol- 
lowing. The choice of Af (as a function of A) will be discussed and specified below. 



The flows of m^, A, Z^-, and Aa- are 



dominated by terms quadratic in G„ for 

A ^ 0, see Fig. 18.101 Using = Aa^, "\ 

we obtain the asymptotic flow equation 

for A from Eg. (18.451) in the simple form r^. o m * ^ ♦ 

^ ^ ' ^ tig. o.lO. Goldstone fluctuations deter- 



d^A = A' [ Gliq) . (8.62) 

Jg\A 



mining the infrared asymptotics for A 0. 



The above integral over is proportional to A"^ for small A in dimensions d < 4, 
implying that A scales to zero in J < 3. In two dimensions one obtains /^|^ Gl(q) = 
47T^A z small A, such that the rescaled variable ~A = AIA obeys the flow equation 

=-1 + -^, (8.63) 



d log A Ati^A^Z^ 

which has a stable fixed point at 1* = Ai{^A,^Zt^. Hence, the bosonic self-interaction 
vanishes as 

A^An^A^Z^A for A^O (8.64) 

in two dimensions in agreement with Eq. (16.371) . Consequently, also the radial mass 
of the Bose fields vanishes linearly in A. In three dimensions one has /^|^ G\{q) oc 

A"' for small A such that A and scale to zero logarithmically for A ^ 0. Since 
is the dominant contribution to the denominator of G^- at small momenta and 

frequencies, the scaling of to zero as a function of A implies that G^- (at A = 0) 

diverges as 



Ga{sq)oz 5-1 forJ = 2 (8.65) 
Gcr(sq) oc log s for d = 3 (8.66) 
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in the limit s ^ matching Eq. (16.351) . Although derived from an approximate trun- 
cation of the functional flow equation, this result is exact even in two dimensions, 
where the renormalization of is very strong. This is due to the fact that the scaling 
dimension of is fully determined by the scaling dimension of the Goldstone prop- 
agator and the existence of a fixed point for A, but does not depend on the position of 
the fixed point (Pistolesi 2004). 
The flow of Zo- is given by 



for small A. The integral over the second derivative of G„Gn is of order A^'^. In two 
dimensions the coupling A vanishes linearly in A, such that OaZo- oc A~^, implying that 
Zg- diverges as A"' for A ^ 0. Hence, the term Zq-^q with l^d = A in the denomina- 
tor of Go- scales linearly in A, as m^. In three dimensions J^|^ d'^pjj„{p + q)G^(q)\p=o 
diverges as A"^, while A vanishes only logarithmically. Hence ^^Zq- oc (log A)~^A~^, 
which is larger than in two dimensions. Integrating over A one finds Z^- oc (A log A)"^, 
which means that Z^-ql vanishes as (log A)~^ in the infrared limit. This yields a sub- 
leading logarithmic correction to the mass term in the denominator of Go-. An 
analogous analysis with the same results as just obtained for Zo- also holds for the 
momentum renormalization factor Ag-. A strong renormalization of longitudinal cor- 
relation functions due to Goldstone fluctuations appears in various physical contexts 
(Weichman 1988, Zwerger 2004). 

Recently, a singular effect of Goldstone fluctuations on the fermionic excitations in 
a superfluid was found in Gaussian approximation (Lerch 2008). This singularity ap- 
pears only after analytic continuation to real frequencies, and its fate beyond Gaussian 
approximation remains to be clarified. 

Since and A scale to zero in the infrared limit in J < 3, all purely bosonic 
contributions to the effective action scale to zero. On the other hand, the fermion- 
boson coupling remains finite. One is thus running into a strong coupling problem, 
indicating a failure of our truncation, if fermionic fields are integrated too slowly, 
compared to the bosons. The problem manifests itself particularly strikingly in the 
flow equation for the order parameter, Eq. (18.431 ). in two dimensions. Since oc Ab 
for small At, the fermionic contribution to dy^a is of order A"' if one chooses A/ = 
Ah = A, leading to a spurious divergence of a for A ^ 0. 

The problem can be easily avoided by integrating the fermions fast enough, choos- 
ing A/ «; Ah in the infrared limit. In our numerical solution of the flow equations 
in the following section we will choose Ah = A and A/ = A^/Ac for A < A^, which 
matches continuously with the equal choice of cutoffs for A > A^. The fermionic con- 
tribution to d^a in Eq. (18.431) is then finite for small A, since the factor Ay = 2A in 
compensates the divergence of m^^ in front of the integral. Since the fermions 




(8.67) 
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are gapped below one could also integrate them completely (set Af to zero), and 
then compute the flow driven by At only. 

The freedom to choose fermionic and bosonic cutoff's independently was exploited 
also in a recent fRG-based computation of the fermion-dimer scattering amplitude in 
vacuum (Diehl 2007). 

8.5.3 Numerical results in two dimensions 

In this section, we present a numerical solution of our flow equations from Sec. 18.4.3] 
in two dimensions. Technically, we employ a fifth-order Runge-Kutta integration rou- 
tine to solve coupled, ordinary diff"erential equations. At each increment of the Runge- 
Kutta routine, two-dimensional integrations over the whole Brillouin zone have to be 
executed. For this purpose, we employ an integrator for singular functions with rel- 
ative error of less then 1%. In particular, for A ~ Ac it is imperative to operate with 
sufficiently accurate routines as the integrands are large and small deviations result in 
a significant spread in the final values for A ^ 0. For further details on the numerical 
procedure, we refer to Appendix lAl 

We fix our energy units by setting the hopping amplitude t = I. We choose a chem- 
ical potential fi = -1.44 corresponding to an average electron density of 1 /2 (quarter- 
filled band). This choice represents the generic case of a convex Fermi surface remote 
from van Hove singularities. The only varying input parameter is the Hubbard U, 
which determines the initial value of the bosonic mass via ml = \2/U\. Initially, the 
flow starts in the symmetric regime with A = Ao = 100, where Eqs. (18.381 - [84TI) 
determine the evolution. The critical scale is determined by the condition m|(Ac) = 0. 
In the symmetry-broken regime (A < Ac), Eqs. (18.431 - 1848]) determine the evolution. 




Fig. 8.11. (a): Flows of the bosonic mass, ml for A > Ac and ml- for A < Ac. (b): 
Quartic coupling A. 
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(a) (b) 




U A 

Fig. 8.12. (a): Fermion gap A, order parameter a, and the critical scale Ac versus 
Hubbard U. (b): Exemplary flows for A (triangles) and a (lines) each corresponding 
to one point in (a). 



In Fig. 18 ■ 1 1 1 characteristic flows of the bosonic mass and the quartic coupling are 
shown in double-logarithmic plots for diff"erent choices of the Hubbard U. The sharp 
de- and increase of the bosonic mass marks the region around Ac. For small A, the flow 
reaches the infrared asymptotic regime (see Sec. 18.5.21) . The scale Air at which this 
scaling sets in decreases for decreasing U. The numerically obtained scaling m^, A oc 
A is consistent with the analytical result of Sec. 18.5.21 

In Fig. 18.121 (a) we compare the fermion gap with the bosonic order parameter 
(final values at A = 0) and the critical scale as a function of U. In Fig. 18.121 (b) the 
flow of A and or as a function of A is shown for various choices of U. We observe 
A, a oc (Ac - Ay^^ for A < Ac as derived below Eq. (18.601) . The ratio A /Ac, where A 
is the final gap for A ^ 0, is approximately 1 .4 for the values of U studied here. As 
a result of fluctuations, the gap is reduced considerably compared to the mean-field 
result 

~ 0.25 (8.68) 

^BCS 

for 2 < i7 < 4. The main reduction here stems from the bosonic self-interactions 
in the symmetric regime leading to a substantial decrease of Ac via the second term 
of Eq. (18.381 ). A reduction of the gap compared to the mean-field value is generally 
present even in the weak coupling limit U ^ 0. In fermionic perturbation theory 
second order corrections reduce the prefactor of the BCS gap formula even forU 
(Martin 1992). The reduction obtained here is slightly stronger than what is expected 
from a fermionic renormalization group calculation (Gersch 2008). 

For A < Ac, Goldstone fluctuations slightly reduce a via the term involving G„(q) 
in Eq. (18.431) . On the other hand, the term due to the Goldstone mode in Eq. (|8.44l) 
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enhances the fermionic gap A relative to a, such that A is generally slightly larger 
than a. This difference will become larger upon increasing the interaction strength as 
one enters the regime of a Bose gas made of tightly bound fermions. Here, however, 
for relatively weak interactions, the impact of Goldstone fluctuations on both a and 
A is very modest. By contrast, the impact of Goldstone fluctuations is known to be 
dramatic at finite temperatures (not treated here) in two dimensions, since they drive 
the order parameter to zero (Goldenfeld 1992). 

In the right plot of Fig. I8.13[ we show flows of the Z- factors of the cr- and 7r-field, 
respectively. In the symmetric regime, the evolution is independent of U. At Ac, the (f>- 
field splits into the cr- and ;r- modes with Za- diverging in the limit A ^ as Z^- oc 
(cf. Sec. 18.5.21) . The Z-factor of the Goldstone field saturates for A <^ A. The flows 
for the A-factors (not shown) parametrizing the momentum dependence of the cr- and 
;r-propagators exhibit very similar behavior. 

Finally, in the left plot of Fig. l8.13l we show flows of the fermion-boson vertices 
gu- and g„ for A < Ac. Their relative changes are only of the order of a few percent 
(note the scale of the vertical axis) with g^. being a bit larger than g„. 



8.6 Conclusion 

Truncating the exact fRG flow, we have derived approximate flow equations which 
capture the non-trivial order parameter fluctuations in the superfluid ground state 
of the attractive Hubbard model, which has been chosen as a prototype model for 
attractively interacting fermions. The superfluid order parameter is associated with 
a bosonic field which is introduced via a Hubbard-Stratonovich decoupling of the 
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fermionic interaction. Below a critical scale Ac, the bosonic effective potential as- 
sumes a mexican hat shape leading to spontaneous symmetry breaking and a Gold- 
stone mode. The bosonic order parameter is linked to but not equivalent to a fermionic 
gap. The fermionic gap is significantly smaller than the mean-field gap, mostly due to 
fluctuations above the scale Ac. Transverse order parameter fluctuations (Goldstone 
mode) below Ac lead to a strong renormalization of radial fluctuations. The radial 
mass and the bosonic self-interaction vanish linearly as a function of the scale in two 
dimensions, and logarithmically in three dimensions, in agreement with the exact be- 
havior of an interacting Bose gas (Pistolesi 2004). On the other hand, the average 
order parameter, the fermionic gap, and the interaction between fermions and bosons 
are affected only very weakly by the Goldstone mode. 

Supplementing the flow equations derived above by a shift of the chemical po- 
tential, to keep the density fixed, one may also try to deal with larger values of U 
(see also Outlook [9. 3. 3 1) . Eagles (1969) and Leggett (1980) have shown that already 
the BCS mean-field theory captures many features of the condensed Bose gas ground 
state made from strongly bound fermion pairs in the limit of strong attraction. Beyond 
mean-field theory, the difference between the fermionic gap and the order parameter 
a increases at larger U. 

It will also be interesting to extend the present analysis to T > 0, in particular in 
view of the possibility of a finite fermionic gap in the absence of long-range order in 
a Kosterlitz-Thouless phase at low finite temperatures (see also Outlook |9.3.2|) . 
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Conclusions 



The topic of this thesis has been the computation of the effects of order parameter 
fluctuations in two- and three-dimensional interacting Fermi systems with empha- 
sis on symmetry-broken phases and quantum critical behavior. For this purpose, we 
have solved flow equations within the functional RG framework for fermionic (single- 
particle) and bosonic (collective) degrees of freedom at zero and finite temperature. 

9.1 Key results 

Three important methodological advancements were achieved during this work. 
9.1.1 Non-Gaussian fixed points 

We demonstrated in chapters [5l[6l and|7]how to compute RG flows of quantum critical 
fermion systems at zero and finite temperature where the interaction vertex is attracted 
towards a non-Gaussian fixed point with nontrivial critical exponents. 

In chapters [5] and [6l we extended the Hertz-Millis theory formulated entirely in 
terms of a bosonic order parameter field to phases with discrete and continuous 
symmetry-breaking. Extending the previous work by Millis (1993), we computed the 
shape of the Tf-line in the vicinity of the QCP for the case where even the zero tem- 
perature theory is described by a non-Gaussian fixed point. We fully captured the 
interplay and relative importance of quantum and thermal fluctuations in the vicinity 
of the QCR 

In chapter |7l we set up a renormalization group framework for coupled theories of 
fermions interacting with their own collective order parameter fluctuations at quantum 
criticality. As a first application, we considered the semi-metal to superfluid QCP of 
attractive Dirac fermions in two dimensions possibly relevant for cold atoms on the 
honeycomb lattice and graphene. At the QCP, the fermion-boson coupling becomes 
relevant and induces anomalous exponents for both the fermion and boson propagator, 
respectively. 
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9.1.2 Goldstone modes 

Systems whose spectrum contains massless Goldstone bosons as a consequence of 
spontaneously breaking a continuous symmetry were analyzed in chapters [6] and [81 

In chapter [6l we devised a truncation of the effective action for zero tempera- 
ture phases with continuous symmetry-breaking (cr-II model) that yields the exact 
momentum behavior of the longitudinal and transversal (Goldstone) mode, respec- 
tively. The Goldstone mode strongly renormalizes the longitudinal propagator while 
all singularities in the Goldstone channel cancel as a consequence of symmetries - 
both modes therefore exhibit strongly different dependences on momenta. An even 
simpler truncation also yields the expected critical exponents directly at the critical 
point where both modes become degenerate. Additionally, the cr-II model permits 
investigation of the finite temperature effects in particular correctly reproducing the 
Mermin-Wagner theorem. 

In chapter[8l we thoroughly studied the impact of the Goldstone mode in fermionic 
superfluids at zero temperature. The Goldstone mode only weakly affects the fermion- 
boson vertex, the fermionic gap, and the order parameter. In the infrared regime, the 
fermions decouple from the flow and the collective bosonic sector described by the 
0--II model yields the exact momentum scaling of the transversal and longitudinal 
mode. 

9.1.3 Non-universal and universal quantities 

In chapters |7] and El the functional RG was implemented as a versatile framework for 
coupled fermion-boson systems to compute universal as well as non-universal quan- 
tities simultaneously. The emergence of -bosonic and universal- effective theories 
from microscopic -fermionic and non-universal- models upon lowering the energy 
scale is a strength of the Fermi-Bose RG approach. In addition to the singular terms 
also collected in other RG methods, our flow equations also capture finite but possi- 
bly large contributions which affect non-universal quantities such as the position of 
a QCP, or the size of a fermionic gap. Furthermore, the energy scale at which the 
universal asymptotics sets in can also be computed conveniently within this frame- 
work enabling a controlled comparison of the size of the various couplings with their 
residual low energy phase-space volume. 

9.2 Criticism 

Despite the methodological advancements of this thesis, some sore spots, subject to 
valid criticism, remain. Three central shortcomings of this work are as follows. 

As outlined in section I27T1 in many microscopic models such as the repulsive Hub- 
bard model the situation of competing ordering tendencies arises. For an accurate 
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computation of non-universal properties and to obtain certain qualitative phenomena 
such as the generation of an attractive d-wave coupling from antiferromagnetic spin 
fluctuations, it is not sufficient to take into account only one ordering channel. 

In chapter m the attractive Hubbard model was considered in a parameter regime 
where the Cooper channel is undoubtedly the most dominant channel and it was there- 
fore justified to utilize only one Hubbard-Stratonovich field for superfluid order. In 
the future, one would want to either (i) decouple the microscopic four-fermion inter- 
action into multiple bosonic channels and a posteriori check how much the results 
depend on the precise decoupling procedure (Baier 2004), or (ii) perform a combined 
fermionic plus bosonic RG study, that is, first integrate out the high energy modes 
with a fermionic RG and subsequently perform the Hubbard-Stratonovich decoupling 
at a scale Ai, when only a few dominant bosonic susceptibilities are left over. Then, the 
results have to be checked on their dependence on the switching scale At- A similar 
strategy has been adopted within a fermionic RG plus mean-field theory approach in 
Reiss (2007). 

The number of couplings employed in the truncations of this thesis has been kept to 
a minimum and often involved merely multiplicative Z-factors for the frequency and 
momentum dependencies of the various propagators. Although it is impressive how 
much interesting physics one can extract already on this level, to ascertain the results 
and improve on quantitative accuracy, it is necessary to parameterize the self-energies 
and vertices on a frequency and/or momentum grid resulting in a more elaborate nu- 
merical eff"ort. In some circumstances, as for example in the superfluid Kosterlitz- 
Thouless phase (see subsection 19.3.21 ). this is necessary to obtain even the correct 
qualitative behavior. 

It is under debate how well the functional RG copes with strongly coupled systems 
(Salmhofer 2007). On the one hand, non-Gaussian fixed points in the infrared regime 
are captured at least qualitatively. The system is strongly coupled in only a limited re- 
gion of phase- space and the small phase- space volume when the cutoff" goes to zero is 
helpful in controlling the vertex expansion when truncating the exact effective action. 
On the other hand, truly strongly coupled systems when the local interaction exceeds 
the bandwidth are strongly coupled everywhere in phase-space and cannot be treated 
in manageable truncations yet (Metzner 2005). 

9.3 Outlook 



There exist numerous interesting directions for future research based on the work 
presented in this thesis. We here list three of them and describe the strategy for an RG 
treatment as well as the operational risks associated with each. 
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9.3.1 QED3 (extending chapter |7l) 

We have shown in chapter |7] how to compute RG flows for coupled fermion-boson 
systems at quantum criticality. An interesting extension concerns Quantum Electro- 
dynamics in 2 + 1 dimensions which has been advocated as a model for the zero tem- 
perature behavior in underdoped high superconductors (Franz 2002, Herbut 2002). 



Massless QED3 can be tuned to a 
QCP by varying the number of fermion 
flavors to a critical A^/,crit- For fewer fla- 
vors than Nf^crit the ground state spon- 
taneously breaks the chiral symmetry 
while for more flavors than A'y;crit the 
gauge interaction is weakened and the 
ground state is chirally invariant. The 
properties of QED3, in particular the 
momentum dependence of the fermion 
and gauge boson propagator {universal) 
at the QCP and A^/,crit {non-universal), 
are presently only poorly known (Kaveh 
2005, Strouthos 2008 and references 
therein). To investigate this in a func- 
tional RG setting, one would replace the 
eff"ective action written in Eq. (17.41) with 
the fermion-boson action of QED3. The 
minimal truncation includes a fermion 
and (gauge) boson self-energy, possibly 
in the form of Z-factors. The flow of the 
renormalized gauge coupling {g) gener- 
ates four-fermion interactions {X f) in the 
chiral sector necessitating two interac- 
tion channels for the truncation. Diverg- 
ing four-fermion interactions signal the 
onset of chiral symmetry-breaking and 
thus permit the determination of A^/,crit 
(Gies 2006). The set of flow equations 
is depicted graphically in Fig. 19. II 





Fig. 9.1. Diagrammatic representation of 
the set of flow equations for QED^. Wiggly 
lines stand for gauge bosons and straight 
lines for fermions. 



Operational risks here involve the fulfillment of gauge symmetry and the depen- 
dence on the choice of gauge-fixing, typically Landau gauge. To get warmed up, 
one may want to start with an eff"ective model for QED3 without gauge bosons -the 
Thirring model in 2 -h 1 dimension (Hands 1995, Christofi 2007, Gies 2008). 
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9.3.2 Superfluid Kosterlitz-Thouless phase (ext. chapter H]) 

The correct finite temperature extension of the truncation employed in the ground- 
state study for two-dimensional fermionic superfluids in chapter[8]must yield the qual- 
itative features peculiar to the Kosterlitz-Thouless phase (Kosterlitz 1973, Kosterlitz 
1974, Goldenfeld 1992). Although there is no spontaneous symmetry-breaking, this 
phase is characterized by a line of fixed points and a temperature dependent anoma- 
lous dimension for the order parameter field. The order parameter correlation function 
decays algebraically due to the presence of a massless Goldstone mode. At the tran- 
sition temperature, the superfluid density jumps to a universal value due to a conspir- 
acy of vortex and Goldstone excitations (Nelson 1977). The single-particle spectrum 
displays a gap around the Fermi level reminiscent of pseudo-gap phenomena in the 
cuprates which also fall into the Kosterlitz-Thouless universality class (Rohe 2001, 
Eckl 2002). 

To date, no satisfactory study of the Kosterlitz-Thouless phase starting from a mi- 
croscopic fermionic model and capturing the above features has been performed and 
would clearly be desirable. First attempts within the functional RG framework for 
bosonic systems have been put forward 
in (Grater 1995, v. Gersdorff 2001). 
However, in both studies the flow runs 
back into the symmetric phase, though 
very slowly, therefore not capturing the -» 
vanishing of the yS-function of the super- 
fluid density in the Kosterlitz-Thouless i^- n i * •; ^ i 
■' rig. 9.2. Contribution to the normal 

phase. Notwithstanding the notorious r ■ ir ■ ui -j 

^ ° jermion self-energy in the superfluid 

difficulties of the bosonic sector, captur- , ^, ■ i v ♦ j ; 

' ^ phase. I he wiggly line stands for the bo- 

ing the gap in the single-particle spec- propagator and the straight line for 
trum without symmetry-breaking seems r 
in reach. Since the anomalous part of 
the fermionic self-energy, that normally 

gaps the spectrum, vanishes, one has to parametrize the frequency- and momentum 
dependence of the normal self-energy of the fermions rather accurately. 

The crucial diagram is shown in Fig. 19. 21 where the order parameter field renormal- 
izes the fermion self-energy. The gap should be visible in the temperature dependence 
of the specific heat deviating strongly from the linear temperature behavior expected 
in a (gapless) Fermi liquid. The scale-dependent specific heat can be computed via 
temperature derivatives of the thermodynamic potential in Eq. (13.221) : 

C\T) = -T . (9.1) 

Operational challenges lie in the numerical implementation of the normal fermion 
self-energy on a momentum- and frequency grid. 
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9.3.3 BCS-Bose crossover (ext. chapter |8]) 

Upon increasing the attractive interaction in a superfluid Fermi gas, the nature of the 
fermion pairs evolves from BCS-type Cooper pairs to tightly bound bosons which 
undergo Bose-Einstein condensation (Leggett 1980). Nozieres and Schmitt-Rink first 
considered this BCS-Bose crossover for lattice fermions and they extended Leggett's 
approach to finite temperatures (Nozieres 1985). 

An interesting extension of chapter[8]concems the properties of the attractive Hub- 
bard model when increasing the attraction U to larger values (Singer 1996, Keller 
1999, Keller 2001, Burovski 2006). The main addition to the RG truncation already 
included in chapter [8] is a constraint to keep the density fixed during the RG flow. 
Implementing a flow equation for the density can be achieved by differentiating the 
flow equation of the thermodynamic potential, Eq. (|3.22l) . with respect to the chemical 
potential: 

dAn^ = -^d^U\ (9.2) 

Ofi 

Including the fermionic, bosonic and fluctuation contributions (Ohashi 2002, Diehl 
2007) to the density is necessary to self-consistently adjust the chemical potential 
to keep the density constant during the RG flow. As a first step, one would try to 
understand and reproduce the Hartree mean-field theory in the zero sound channel 
within the RG equations analogously to the mean-field BCS model (Salmhofer 2004). 
Subsequently, a Hartree plus BCS-model with fixed density would have to be, at least 
qualitatively, reproduced by the RG flow. Finally, the order parameter fluctuations are 
added in and the effects thereof studied. 
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Deutsche Zusammenfassung 



Die vorliegende Dissertation entwickelt Renormierungsgruppenstrategien zur Berech- 
nung makroskopischer Eigenschaften von wechselwirkenden Fermi-Systemen. Beson- 
ders das Verhalten bei niederen Energien ist interessant, da hier wechselwirkungsin- 
duzierte, kollektive Phanomene auftreten. Durch Herleitung und Losung von Re- 
normierungsgruppengleichungen analysieren wir in dieser Arbeit den Einfluss von 
Ordnungsparameterfluktuationen auf thermodynamische Observable und Korrela- 
tionsfunktionen in Fermi-Systemen mit spontaner Symmetriebrechung und quan- 
tenkritischem Verhalten. 

Die Dissertation ist in zwei Hauptteile gegliedert. Teil IJ stellt das theoretische 
Riistzeug der Arbeit bereit; in Teil |II] wird dieses dann auf verschiedene Problem- 
stellungen angewendet. 

In Kapitel[2]werden fundamental Konzepte der Physik korrelierter Fermionen und 
der Renormierungsgruppe, die zum Verstandnis der in Teil |ll] folgenden Anwendun- 
gen notwendig sind, vorgestellt. Das Herzstiick des ersten Teils stellt Kapitel [3] dar, 
in dem die einteilchenirreduzible Darstellung der funktionalen Renormierungsgruppe 
hergeleitet wird. 

In Abschnitt 13 .1.11 wird erlautert wie man die in Superf elder gebiindelten fermion- 
ischen und bosonischen Freiheitsgrade in den Formalismus einbringt, so dass man 
simultan zu den fermionischen auch die koUektiven, bosonischen Fluktuationen aus- 
integrieren kann; dies wird in den Kapiteln|7]und[8]in die Praxis umgesetzt. 

Im Hinblick auf Anwendung in den Kapiteln[5l[6]und[8l erweitem wir in Abschnitt 
|3.2.2| den Formalismus auf symmetriegebrochene Systeme. In Abb. |3.1| zeigen wir die 
durch einen nichtverschwindenden Erwartungswerts einer bosonischen Feldkompo- 
nente im Minimum des effektiven Potentials modifizierte Hierarchic der Vertexfunk- 
tionen. 

Der zweite Teil der Dissertation ist vier verschiedenen Anwendungen der in TeilU 
vorgestellten Konzepte und Methoden gewidmet. 
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In Kapitel [5] verallgemeinern wir die urspriinglich von Hertz und Millis entwick- 
elte Renormierungsgruppentheorie, die Fermionen auf einen Schlag ausintegriert und 
die daraus resultierende efFektive, bosonische Theorie mit Hilfe von Flussgleichun- 
gen analysiert, auf Phasen mit diskreter Symmetriebrechung. Im Gegensatz zu Millis 
konnen wir durch das Losen unserer Flussgleichungen bei endlichen Temperaturen 
auch stark wechselwirkende, nicht-GauBsche Fluktuationen behandeln, die das Ver- 
halten entlang der im quantenkritischen Punkt miindenden Phasengrenze dominieren. 

Im ersten Schritt studieren wir bosonische Propagatoren mit dynamischen Expo- 
nenten z > 2 in zwei und drei Raumdimensionen. Obwohl das Infrarotverhalten bei 
endlichen Temperaturen durch einen nicht-GauBschen Fixpunkt beschrieben wird, 
bestatigen wir das von Millis in GauBscher Approximation von der symmetrischen 
Phase kommend abgeleitete Potenzgesetz, das die Phasengrenze als Funktion des 
KontroUparameters in der Nahe des quantenkritischen Punktes beschreibt. Dies ist 
eine Konsequenz des geringen Phasenraumvolumens, das den EfFekt der starken Kop- 
plung im Infraroten unterdriickt. Da wir die Phasengrenze durch Ausloten der Position 
des Fixpunktes bestimmen, sind wir in der Lage, die kritische Region zwischen der 
von Millis als Phasengrenze angenommenen Ginzburg-Linie und der wahren Phasen- 
grenze numerisch auszumessen: Wahrend beide Linien in drei Raumdimensionen 
kaum voneinander abweichen, ofFnet sich in zwei Raumdimensionen eine ausgepragte 
Liicke zwischen Ginzburg-Linie und Phasengrenze bei endlichen Temperaturen. 

Im zweiten Schritt nehmen wir uns im Abschnitt 15.5.31 das bis dato nicht unter- 
suchte Quanten-Ising-Modell in zwei Raumdimensionen mit dynamischem Exponent 
z = I vor. Nun liegt auch die Grundzustandstheorie unterhalb der oberen kritischen 
Dimension, was die von Millis gemachte Naherung -um den nicht-wechselwirkenden 
GauBschen Fixpunkt zu entwickeln- voUstandig invalidiert. Wie in Abb. 15.61 illus- 
triert, weicht in diesem Fall die Ginzburg-Linie auch bei Temperatur null von der 
Phasengrenze ab. Resultate der numerischen Losung der Flussgleichungen sind in 
Abb. l5.7l und l5.8l abgebildet. Fiir die Phasengrenze als Funktion des KontroUparame- 
ters linden wir Wurzelverhalten mit einer logarithmischen Korrektur, siehe Gl. (15.491) . 

In Kapitel [6] analysieren wir Falle, wo der quantenkritische Punkt einer Univer- 
salitatsklasse zugehorig ist, in der eine kontinuierliche Symmetriegruppe gebrochen 
wird. Dann treten masselose Goldstone-Moden auf, die das kritische Verhalten radikal 
beeinflussen -besonders in niedrigen Dimensionen. 

Basierend auf dem in Abschnitt 16 . 21 definierten cr - 11 Modell, dass den Goldstone- 
Anregungen mit einem separaten Propagator Rechnung tragt, fiihren wir analog zu 
Kapitel [5] eine Flussgleichungsstudie in der Nahe eines quantenkritischen Punktes 
durch. Wir fokussieren uns dabei auf Systeme mit quadratischer Frequenzabhangigkeit 
-also dynamischer Exponent z= 1- in den Propagatoren beider, radialer und transver- 
saler, Anregungen des Systems. Die Flussgleichungen spiegeln das Mermin- Wagner- 
Theorem richtig wider: Der Ordnungsparameter, und damit die kritische Temperatur, 
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wird in Raumdimensionen kleiner gleich zwei von thermisch angeregten Goldstone- 
Moden zu null gedriickt. 

Als ein unerwartetes Hauptresultat linden wir in Abschnitt l6.4.3l dass die Form der 
Phasengrenze bei endlicher Temperatur in der Nahe des quantenkritischen Punktes in 
drei Raumdimensionen nicht von der Anzahl der Goldstone-Moden abhangt. 

Abseits der Phasengrenze, in der symmetriegebrochenen Phase, extrahieren wir 
in Abschnitt 16.51 das exakte Infrarotverhalten beider Progagatoren und der Wechsel- 
wirkungen analytisch aus den Flussgleichungen. Die effektive Selbstwechselwirkung 
der Goldstone-Moden llieBt aus Symmetriegriinden zu Null (siehe Gl. (16.341 )). und 
der Goldstone-Propagator behalt seine unrenormierte, quadratische Frequenz- und 
Impulsabhangigkeit bei. Die radiale Mode und deren Wechselwirkungen hingegen 
werden, wie in Gl. (16. 35116.371) zusammengefasst, stark durch die Goldstone-Moden 
renormiert. 

In Kapitel|7]berechnen wir Renormierungsgruppenfliisse gekoppelter Fermi-Bose 
Theorien, die sich durch ihr kompliziertes, und ob der beiden masselosen Propa- 
gatoren am quantenkritischen Punkt daher mit anderen Methoden bisher nur unzu- 
reichend erforschbares, Infrarotverhalten auszeichnen. 

Wir analysieren attraktiv wechselwirkende Dirac-Fermionen (siehe Abb. 17. 21) : hier 
treibt die Wechselwirkungstarke bereits in Molekularfeldtheorie einen kontinuier- 
lichen Quantenphaseniibergang von einem Halbmetall zu einer Superflussigkeit (siehe 
Abb. 17. II und Abschnitt 17.2.11) . Physikalisch treten Dirac-Fermionen in Kohlenstoff- 
schichten und in optischen Honigwabengittern in kalten, atomaren Quantengasen auf. 

Unsere Trunkierung und Flussgleichungen der Fermion-Boson Wirkung, in Abb. 
17.41 diagrammatisch dargestellt, ermoglichen die Berechnung aller quantenkritischen 
Exponenten in unmittelbarer Nahe des quantenkritischen Punktes. Abb. 17.51 zeigt 
die Werte der anomalen Dimensionen der Fermionen und Boson direkt am quan- 
tenkritschen Punkt fiir verschieden Raumdimensionen. Wie unter Gl. (17.331 ) erlautert, 
impliziert eine anomale Dimension des Fermionfeldes den Zusammenbruch des Fer- 
mifliissigkeitsverhalten: Die lineare Frequenzabhangigkeit der Selbstenergie wird zu 
einem fraktionellen Potenzgesetz, das die hohe Zerfallsrate der niederenergetischen 
Quasiteilchen in der Nahe des Dirac-Punktes beschreibt. 

In Abschnitt 17.4.11 werden numerische Fliisse der in der Trunkierung vorkom- 
menden Parameter prasentiert. Durch kontroUiertes Verkniipfen der mikroskopischen 
Anfangsbedingungen, die auch der Molekularfeldrechung von Abschnitt 17.2.11 zu- 
grundeliegen, ist es uns moglich, zusatzlich zu den universellen Exponenten auch 
die nicht-universelle, renormierte Position des quantenkritischen Punktes entlang 
der KontroUparameterachse mit dem Wert der Molekularfeldrechnung zu vergle- 
ichen; Gl. (17.391 ) veranschaulicht, dass Quantenfluktuationen die AusmaBe des Su- 
perflussigkeitsbereichs im Phasendiagram Abb. l7.1l drastisch reduzieren. 
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Die Berechnung des Satzes quantenkritischer Exponenten wird durch die Bestim- 
mung der Potenzgesetze fiir die Korrelationslange und die Suszeptibilitat als Funktion 
des Kontrollparameters in Ab schnitt 17 . 4 . 21 ab gerundet . 

In Kapitel [8] wird ein Gros der in den Kapiteln [5] bis |7] aufgebauten Neuerun- 
gen synthetisiert, und wir prasentieren eine umfassende Analyse des attraktiven 
Hubbard-Modells bei Viertelbandfiillung, als Prototyp eines fermionischen Systems 
mit suprafluidem Grundzustand. In einem gekoppelten Fermion-Boson Fluss, der 
sowohl die Selbstwechselwirkungen des Ordnungsparameters, zwei verschiedene 
Fermion-Boson Kopplungen, als auch die anomalen Komponenten der fermionischen 
Selbstenergie beriicksichtigt, werden simultan zu diversen nicht-universellen GroBen, 
wie der Energieliicke des fermionischen Energiespektrums, auch universelle Infrarot- 
eigenschaften des koUektiven, bosonischen Sektors bestimmt. 

Strukturell zeichnet sich unsere Analyse durch Unterscheidung des Minimums des 
bosonischen, effektiven Potentials und der fermionisches Energieliicke -beide fallen 
nur in Molekularfeldtheorie zusammen, siehe Abschnitt l8.4.3l - und Trennung der ra- 
dialen und transversalen Moden (siehe Abschnitt l8.4.3l) . aus. Letzteres ermoglicht in 
Anlehnung an Abschnitt 16.51 die exakte Behandlung des Infrarotsektors (siehe Ab- 
schnitt[832l)- 

In Ab schnitt 18.5.31 werden numerische Losungen der gekoppelten Renormierungs- 
gruppengleichungen in zwei Raumdimensionen prasentiert. Ein wesentliches Resultat 
ist der in Abb. 18.1 2l dargestellte Vergleich zwischen fermionischer Energieliicke, Min- 
imum des effektiven Potentials und kritischer Skala. Die Werte fiir die fermionische 
Energieliicke sind um etwa einen Faktor vier gegeniiber der Molekularfeldtheorie re- 
duziert; hauptsachlich verursacht durch Ordnungsparameterfluktuationen, die bei rel- 
ativ hohen Energieskalen bereits im symmetrischen Teil des Flusses der Formierung 
suprafluider Ordnung entgegenwirken. 

Die linke Graphik in Abb. l8.13l zeigt. dass beide Fermion-Boson Vertizes, obschon 
diese durch Goldstone-Fluktuationen renormiert werden (siehe Abb. 18.91) . quantitativ 
nur um wenige Prozentpunkte von ihren nicht-renormierten Ausgangswerten abwe- 
ichen. Dies ist eine Konsequenz des geringen den Goldstone-Moden zur Verfiigung 
stehenden Phasenraumvolumens. 

Der unterschiedliche Verlauf der bosonischen Z-Faktoren ist in der rechten Abbil- 
dung von Abb. 18.131 sichtbar: Wahrend der Z-Faktor der radialen Mode als Funktion 
der Skala divergiert, saturiert der Z-Faktor der Goldstone-Mode, so dass diese ihre 
quadratische Frequenz- und Impulsabhangigkeit behalt -wie bereits in Abschnitt 16.51 
analytisch demonstriert. 

Die Dissertation schliesst mit einer kritischen Diskussion in Kapitel |9l die sowohl 
methodische Errungenschaften als auch Optimierungspotenziale der Arbeit und in- 
teressante zukiinftige Forschungsvorhaben umreiBt. Im Anhang |A] wird die zur nu- 
merischen Losung der Flussgleichungen entwickelte numerische Prozedur erklart. 



Part IV 



Appendices 



A 



Numerical procedure 



In this Appendix, we present the programs which were assembled and subsequently 
used to solve the RG flow equations in this thesis. The exact flow equation, Eq. (13.181) , 
projected on a finite number of parameters represents a system of coupled ordinary 
difi'erential equations. The right-hand- sides of these equations involve numerical inte- 
grations. In chapter[5J the right-hand-sides involve one-dimensional integrations and 
in chapter [8] two-dimensional momentum integrations over the whole Brillouin zone 
have to be executed. 

Hence, a routine for coupled difi'erential equations was combined with a numerical 
integration routine which provides the values of the right-hand-sides at each increment 
of the difi'erential equation solver. A data flow picture is exhibited in Fig. lA.l[ 

The routines employed were taken from the GSL-GNU scientific library Q and then 
assembled and modified in a C/C-1--1- code. Specifically, we used the Runge-Kutta- 
Fehlberg (4, 5) stepping function, gsl_odeiv_step_rkf45, for the difi'erential equation 
solver. This routine adaptively adjusts the step-size for given error margins and proved 
to be computationally efficient as well as accurate even at very low temperatures. We 
used relative error margins of le-4 to le-S. 

For the numerical integrations, we employed the routine gsl Jntegration_qags. This 
routine performs well for singular integrands as we encounter at the critical scale or 
toward the end of the flow. It adaptively adjusts the stepsize within each integration 
interval to satisfy the error margins. We employed relative error margins of le-2 to 
le-4. Two-dimensional integrations were executed by linearly combining two one- 
dimensional routines. 



^ http://www.gnu.org/software/gsl/ 
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Runge-Kutta Solver 

rhs[1]=f_1 [y[1 ],..., y[n]] 
rhs[2]=... 



rhs[n]=f_n[y[1 ],..., y[n]] 



Infrared Loop 

y_final[1 ],..., y_final[n] 



Integration Routine 


lnt_1{y[1],. 




Int_2{ y[1],. 




Int_n{ y[1],. 





Output-file1: values 
vs. scale 



Output-file2: final 
Infrared values 



Main Loop 

Initial conditions: 
y_lnl[1],...,yjnl[n] 
and physical parameters 
e.g.: T, density 



Fig. A.I. Schematic dataflow of the program employed to solve the flow equations 
in this thesis. The maximum number of coupled equations solved simultaneously in 
chapter\8\was n < 15. Typical computation times for forty consecutive flows as shown 
in Fig. \8.12\ were of the order of several days up to one week on a single modern 
processor 
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